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SECTION  1 
INTRODUCTION 


BACKGROUND 

The  United  States  Air  Force  has  launched  a program  directed  towat^l  the  devel- 
op(nent  of  a fourth-generation  Intercontinental  Ballistic  Missile  (ICBM)  desig- 
nated I'/.vi.  This  develofwient  program  is  aimed  at  producing  the 

technology  and  launch  concepts  for  a large  throw  weight,  survivable,  highly 

* 

accurate  MIRVed  ICBM  (ref.  1).  This  improvement  in  missile  technology,  spe- 
cifically that  related  to  accuracy,  has  led  to  the  consideration  of  different 
missile-basing  concepts.  Present  basing  concepts  for  Minuteiiwn  rely  upon 
hardened  silos.  The  perceived  possible  improvements  in  missile  accuracy  will 
radically  increase  the  hardness  requirentent  for  these  silos  or  the  alternative 
fixed-point  basing  concepts.  This  increase  can  be  both  impractical  and  ex- 
tremely costly.  Thus  the  Air  Force  is  considering  what  is  popularly  called 
the  nhi  /./  i^hv'u!  concepts.  This  simply  involves  an  areal  dispersion  of  poten- 
tial targets  not  all  of  which  contain  a missile;  e.g.,  two  of  the  more  popular 
concepts  are  the  shelter  and  the  trench.  The  shelter  concept  involves  provid- 
ing a number  of  semihardened,  relatively  inexpensive  shelters  for  each  missile. 
These  shelters  are  dispersed  over  a given  area  and  the  missile  is  occasionally 
moved  among  them.  Thus  the  probability  of  a direct  or  near  direct  hit  is  sim- 
ply the  inverse  of  the  number  of  structures.  In  the  trench  concept,  a long, 
shallow,  underground  tunnel  is  provided  for  the  launch  vehicle  and  the  missi'e 
is  randojtily  moved  to  different  locations  within  the  tunnel.  In  this  case  the 
probability  of  a near  miss  is  inversely  proportional  to  the  length  of  the 
tunnel . 

The  concept  of  areal  dispersion  of  potential  tat'gets,  as  opposed  to  superhard 
fixed  targets,  changes  the  mechanism  governing  the  vulnerability.  With  super- 
hardened,  fixed-point  structures  the  prime  concern  is  directed  toward  cratering 

★ 

Multiple  Independent  Reentry  Vehicle. 

1.  Slay,  Alton  D.,  Lt.  Gen.,  USAF,  "MX,  A New  Dimension,"  .l;)* 

September  1976. 
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and  high-stress,  ground-shock  propagation.  HoMever,  with  areal  dispersion  the 
emphasis  shifts  to  the  airblast  since  the  missile  Is  most  likely  to  be  some 
distance  from  the  crater  region.  It  Is  for  this  reason  that  the  Air  Force  has 
been  developing  a Dynamic  Airblast  Simulator  (DABS).  This  simulator  Is  being 
used  to  test  the  effects  of  airblast  loading  on  some  of  the  structures  to  be 
used  In  different  basing  concepts  for  MX. 

DABS  Is  sometimes  referred  to  as  a disposable  shook  tube.  It  simply  consists 
of  a roofed  tunnel  covered  with  an  overburden  which  Is  driven  by  high  explo- 
sives at  one  end.  It  Is  made  as  short  as  possible  for  cost  reasons.  Model 
launch  structures  are  located  on  the  floor  of  the  facility  and  their  response 
to  the  blast  environment  Is  measured. 

Since  the  basic  purpose  of  DABS  Is  to  generate  a nonnuclear  airblast,  simula- 
ting the  blast  from  a nuclear  device  at  a given  range,  the  blast  wave  produced 
must  adequately  represent  that  wave  In  both  magnitude  and  waveform.  This  Is 
somewhat  difficult  with  DABS  because  of  the  short  tunnel.  It  Is  necessary  that 
the  expanding  detonation  products,  driving  a clean  air  shock,  be  used  to  pro- 
vide a significant  portion  of  the  blast  loading  on  the  structure  being  tested. 
Herein  lies  the  basic  difficulty  for  both  the  theoretician  and  the  experimen- 
talist. The  thermodynamic  equations  for  shocked  air  are  quite  well  known.  This 
Is  not  the  case  for  the  detonation  products;  no  really  adequate  equatlon-of- 
state  exists.  With  no  adequate  model,  the  burden  of  parameterizing  this  por- 
tion of  the  flow  falls  upon  the  experimentalist.  However,  the  difficulty  of 
making  measurements  In  the  detonation  products  is  quite  great,  since  the  at- 
mosphere Is  composed  of  unknown  constituents  at  unknown  temperatures  and  Is 
loaded  with  debris.  As  far  as  the  theoretician  Is  concerned,  the  prime  vari- 
ables that  need  to  be  measured  In  reference  to  the  shock  wave  are  the  side-on 
or  Incident  pressure,  the  shock  velocity,  and  the  dynamic  pressure  In  the  shock 
wave.  The  first  two  are  quite  easily  and  routinely  measured;  the  third  Is  not 
so  easily  obtained.  If  one  defines  p as  the  density  In  the  flow,  V as  the  flow 
velocity,  and  Q as  the  dynamic  pressure,  Q = pV^/2. 
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OBJECTIVE 


The  prospects  for  measuring  the  dynamic  pressure  directly  are  quite  slim.  The 
flOM  velocity  can  be  measured  but  there  Is  no  reasonable  prospect  for  measuring 
the  density  of  this  dirty  flow.  Thus*  this  Instrumentation  development  was  di- 
rected toward  the  creation  of  measuring  techniques  that  could  parameterize  the 
flow  In  such  a manner  as  to  obtain  the  variables  from  which  the  dynamic  pres- 
sure could  be  Inferred.  These  variables  Include  the  flow  velocity,  the  velo- 
city of  the  Interface  between  the  air  and  the  detonation  products,  stagnation 
pressure,  temperature,  mach  number,  and  driver-chamber  pressure.  Methods  to 
accomplish  this  objective  were  thus  Investigated. 

This  Investigation  was  carried  out  on  a series  of  test  events,  termed  the 
DABS~i  Series . The  purpose  of  this  series  was  twofold- -to  develop  the  tech- 
niques necessary  to  generate  an  appropriate  blast  waveform  and  to  create  an 
environmerrt  In  which  the  developmental  Instrumentation  could  be  tested. 
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SECTION  2 
MICROWAVE  SYSTEMS 


DOPPLER  RADAR 

Introduction 

The  use  of  relatively  inexpensive  tra<  ic-control  doppler  radar  systems  to  make 
velocity  measurements  in  shock  waves  is  clearly  one  of  the  more  novel  and  sig- 
nificant advances  in  the  instrumentation  development  program  for  DABS.  A mass- 
produced,  commercially  available  item,  with  minor  modification,  has  been  adapted 
to  high  velocity  experimental  work.  Thus,  a relatively  sophisticated  microwave 
system  can  be  obtained  at  a cost  less  than  that  of  a normal  blast-pressure  gage. 

Radar  experiments  were  fielded  on  DABS  lA,  IC,  ID,  and  IE.  The  first  of  these 
experiments  on  DABS  lA  was  previously  reported  at  the  Defense  Nuclear  Agency 
Blast  and  Shock  Instrumentation  Conference  at  the  Waterways  Experiment  Station 
in  October  1975.  Briefly,  this  paper  (appendix  A)  describes  how  a radar  unit 
was  used  to  track  an  aluminum-foil  target  after  it  was  impacted  by  the  shock 
wave.  The  resultant  velocity  was  4489  ft/sec.  This  velocity  corresponds  with 
the  expected  particle  velocity  in  air  shocked  to  327  psi , a value  which  agreed 
with  the  pressure  readings  and  with  the  measured  shock  velocity.  This  Lagran- 
gian  particle  velocity  measurement,  although  quite  useful  in  checking  out  the 
system,  only  proved  that  the  ratio  of  the  isobaric  to  isovolumetric  specific 
heats  was  1.4,  This  was  already  known  for  clean  air  and  thus  the  experiment 
simply  confirmed  this  but  served  more  as  a check  on  the  measurement  system. 

One  of  the  prime  problems  associated  with  DABS  is  that  the  explosive  products 
are  used  as  a significant  portion  of  the  loading.  The  equation-of-state  for 
these  gaseous  products  combined  with  debris  is  unknown.  Thus,  making  meas- 
urements that  parameterize  the  explosive  products  is  extremely  important.  It 
is  to  this  end  that  the  additional  radar  experiments  were  directed. 

Possibly,  the  interface  between  the  explosive  products  and  the  driven  clean  air 
shock  can  be  tracked  by  a single  radar  unit  looking  down  the  DABS  tunnel.  The 
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primary  question  to  be  resolved  was  whether  this  interface  would  be  reflective 
or  transparent.  As  shown  in  appendix  A,  the  prime  variable  affecting  the  re- 
flectivity is  the  free-electron  concentration  in  the  hot  gases.  When  the  inci- 
dent microwave  radiation  is  10.525  GHz,  the  interface  will  be  reflected  when 
the  free-electron  concentration  is  high  enough  to  cause  a plasma  frequency  of 
10.525  GHz.  In  the  clean  air  portion  of  the  shock  wave  the  electron  concentra- 
tion at  which  the  peak  pressures  are  nominally  600  psi,  the  free-electron 
concentration  does  not  achieve  this  value.  However,  at  the  interface  between 
the  explosive  products  and  the  clean  air  shock  the  situation  is  quite  different. 
The  explosive  used  as  the  driver  (PETN)  is  oxygen  deficient,  and  as  a result,  a 
state  of  combustion  exists  at  and  near  the  interface.  This  combustion  is  con- 
tinually supported  (during  the  time  frame  of  interest)  primarily  by  the  oxygen 
in  the  air  mixing  with  the  detonation  products  as  a result  of  turbulence.  The 
characteristics  of  this  interface  are  controlled  by  the  rate  of  reaction  of  the 
unknown  nonequilibrium  detonation  products  with  the  air.  Thus,  the  plasma 
characteristics  of  the  interface  are  determined  by  the  rate  of  mixing  along  the 
interface  and  what  the  nonequilibrium  products  are.  The  present  data  preclude 
the  possibility  of  determining  these  quantities.  The  only  facts  known  are  that 
the  nonequilibrium  products  burn  (photographically  observed)  and  that  this 
burning  probably  occurs  in  a manner  consistent  with  the  formation  of  a large 
number  of  free  electrons.  These  electron  concentrations  could  easily  approach 
or  even  exceed  the  3.5  x 10‘®  electrons/cm®  needed  to  totally  reflect  the  radar 
beam.  Obviously,  this  hypothesis  needed  experimental  investigation.  Thus  a 
radar  experiment  was  fielded  on  DABS  IC  to  test  this  hypothesis.  The  results 
of  this  experiment  were  used  in  the  subsequent  experiments  to  direct  modifica- 
tions. 

Radar  Unit 

The  radar  unit  used  in  all  the  experiments  was  a Model  JFIOO  Speed  Gun  manufac- 
tured by  CMI,  Inc.,  Minturn,  Colorado.  In  its  normal  configuration  it  is  a 
hand-held,  traffic-control,  doppler  radar  unit  which  digitally  displays  veloc- 
ity in  miles  per  hour.  The  pertinent  specifications  as  given  by  the  manufac- 
turer are  as  follows: 

(1)  Frequency  — 10,525  ± 25  MHz 

(2)  Polarization  — Circular 
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(3)  Beam  Width  --  8°  or  Less  (Side  lobes  suppressed  greater 

than  24  dB  down) 

(4)  Antenna  — 4-Inch  Aperture;  Circular  Horn  Type 

(5)  Microwave  Oscillator  — Solid-State  Gunn  Diode 

(6)  RF  Power  — 20  mW  Minimum;  100  mW  Maximum 

(7)  Receiver  Diode  --  Shottky  Barrier  Type  Rated  for  100-mW  Burnout 

The  unit  was  ordered  without  the  digital  display  or  Its  associated  electronics 
since  they  were  not  pertinent  to  the  experiments.  Thus,  the  total  unit  consis- 
ted of  an  antenna,  a preamplifier,  an  ampllfler/regulator,  a 9-V  regulator,  and 
the  housing.  The  system  was  powered  by  a 12-V  wet  cell. 

Experimental  Procedures 

In  DABS  lA  the  signal  was  taken  from  the  output  of  the  preamplifier  and  ac- 
coupled  (high  pass  30  Hz,  3-dB  point)  because  of  a dc-offset.  The  only  other 
modification  made  was  the  removal  of  all  low-pass  filtering  from  the  preampli- 
fier. The  output  Impedance  of  the  system  was  1000  n. 

The  relatively  high  output  Impedance  In  DABS  lA  prompted  a change  to  the  radar 
system  for  DABS  IC.  To  Improve  the  frequency  response,  a buffer  amplifier  pow- 
ered by  a 45-V  battery  (fig.  1)  was  used  as  a line  driver.  This  amplifier  was 
an  Inverting  unity-gain  device  with  a 100D-f2  Input  Impedance  and  a 75-^}  output. 
Impedance.  This  modification  negated  the  need  for  ac-coupling  on  the  output. 
After  this  modification  was  made,  the  system  was  checked  out  with  a rifle  bul- 
let as  a high-velocity  target.  (The  procedure  described  In  appendix  A was 
used.)  The  results  (fig.  2)  agreed  quite  well  with  the  previous  rifle  data 
and  the  expected  bullet  velocity.  Also  the  previous  offset  problem  described 
In  appendl;  A was  eliminated  by  the  buffer  and  a terminated  line. 

The  output  of  the  radar  antenna  was  ac-coupled  (high  pass)  to  the  Input  of  the 
preamplifier  for  all  experiments  after  those  on  DABS  IC.  This  was  accomplished 
by  a simple  one-pole  RC  filter  which  was  formed  by  connecting  a 0.005  uF  capaci 
tor  In  series  between  the  antenna  and  the  preamplifier.  Since  the  preamplifier 
was  1000  ri,  this  formed  a high-pass  filter  with  a 3-dB  point  of  32  kHz.  The  re 
corded  signal  was  the  amplified  doppler  signal,  and  Its  relationship  to  the 
velocity  of  the  target  was  f^  * 2Vf^/C,  where  f^  Is  the  doppler  frequency. 
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Figure  1.  Radar  Buffer  Amplifier 


V is  the  velocity  of  the  target,  f^  is  the  radar  carrier  frequency,  and  C is 
the  speed  of  light.  If  the  unit  of  velocity  is  feet  per  second  and  the  fre- 
quency is  in  Hertz,  this  equation  reduces  to  fj  » 21.40V. 

For  each  of  the  experiments  the  data  were  recorded  on  a Bell  & Howell  3700  B 
Tape  Recorder.  With  the  exception  of  the  DABS  lA  data,  all  data  were  recorded 
with  Wideband  II  FM  (500-kHz  frequency  response).  On  DABS  lA  the  signal  was 
transmitted  by  RG213  cable  which  was  not  terminated.  This  was  necessary  be- 
cause of  the  high  output  impedance  of  the  preamplifier;  however,  this  was  not 
necessary  on  the  subsequent  experiments,  since  the  buffer  amplifier  was  used. 

A terminated  RG58  line  was  used  for  these  tests.  In  all  the  experiments  the 
line  was  approximately  1000  ft  long. 

In  all  cases,  the  radar  unit  was  aimed  down  the  center  of  the  tunnel  and  placed 
so  that  debris  created  after  impact  by  the  shock  wave  would  not  Interfere  with 
other  experiments  in  the  DABS  facility.  No  attempt  was  made  to  harden  the 
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(a)  Bullet  Fired  Toward  Radar  Unit 
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(b)  Bullet  Fired  Past  and  Away  From  Radar  Unit 


Figure  2.  Doppler  Radar  Signal  from  7-mm  Magnum  Bullet 


IS 


Figure  3.  Radar  Unit 


radar  unit  to  the  effects  of  the  shock  wave.  This  would  have  Involved  build- 
ing a bunker  with  a thick  quartz  or  Teflon  window  for  the  radar  unit  to  look 
through;  the  cost  would  far  exceed  the  cost  of  the  radar  unit.  Thus,  the  ra- 
dar was  presumed  to  be  disposable  and  was  simply  mounted  on  a wooden  stake 
(fig.  3),  usually  Inside  the  DABS  tunnel. 

With  the  exception  of  DABS  lA,  there  was  no  metal  target  placed  In  the  tunnel, 
since  the  objective  was  to  track  the  Interface  between  the  detonation  products 
and  the  clean  air  shock.  Thus,  the  DABS  IC  radar  test  was  performed  to  deter- 
mine whether  or  not  the  Interface  between  the  detonation  products  and  the  clean 
air  shock  would  reflect  the  microwave  radiation.  The  possibility  of  another  type 
of  microwave  system  (discussed  later)  prompted  Interest  In  the  propagation  of  a 
microwave  beam  of  10.525  GHz  frequency  through  the  detonation  products.  This 
additional  Information  made  the  experiment  quite  unique  from  the  experimental- 
ist's point  of  view.  The  fact  that  the  wave  could  propagate  through  the  shock 
front  had  been  established  In  DABS  lA.  Thus,  there  were  only  three  possibil- 
ities: (1)  the  beam  would  be  reflected  by  the  Interface  between  the  detonation 
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products  and  the  clean  air  shock;  (2)  the  beam  would  be  reflected  only  by  the 
back  walls  and  thus  be  transmitted  through  the  detonation  products;  or  (3}  the 
beam  would  be  totally  absorbed  by  the  detonation  products  and  thus  there  would 
be  no  return  signal.  Obviously,  the  first  two  possibilities  are  the  most  de- 
sirable, but  all  three  cases  represent  meaningful  data.  Thus,  excluding  the 
possibility  of  a severed  cable,  even  no  return  signal  would  be  meaningful. 

Results 

The  results  of  the  DABS  IC  radar  test  were  quite  Interesting  In  that  they 
showed  that  the  radar  unit  did  track  the  Interface  between  the  detonation  prod- 
ucts and  the  clean  air  shock  and  also  that  the  beam  was  transmitted  through  the 
Interface  and  reflected  off  the  rear  wall.  The  raw  data  are  shown  In  figure  4. 
These  data  show  that  there  were  high-frequency  data  superimposed  on  low- 
frequency  data.  The  reason  for  this  can  be  seen  by  means  of  the  overlayed  draw- 
ing on  the  raw  data.  The  low  frequencies  represent  velocities  that  varied  from 
12  to  39  ft/sec.  These  values  qualitatively  agreed  with  the  photopole  data. 
Thus,  they  represent  the  velocity  of  the  wall  behind  the  driver  section  In  DABS. 
The  obvious  clipping  of  the  data  was  caused  by  the  radar  unit  Itself— specifi- 
cally the  maximum  output  voltage  of  the  Integral  preamplifier.  Thus,  all  data 
above  or  below  approximately  ± 1.5  V were  lost  because  of  clipping  In  the  pre- 
amplifier. This  led  to  a reduction  In  the  number  of  high-frequency  data  points, 
which  are  presumed  to  represent  the  velocity  of  the  Interface. 

The  reduced  doppler  frequency  Is  shown  In  figure  5 along  with  the  difference 
between  It  and  the  shock  velocity.  This  latter  curve  can  be  Integrated  to  de- 
termine the  distance  between  the  clean  air  shock  front  and  the  Interface  be- 
tween the  detonation  products  and  clean  air  shock.  These  data  were  obtained 
by  measurlng^the  period  of  the  oscillations  and  converting  that  to  frequency. 
This  was  accomplished  with  50-usec/1n  plots. 

The  knowledge  that  the  10.5-GHz  beam  could  propagate  through  the  detonation 
products  was  useful  as  a design  criterion  for  the  second  microwave  system  In- 
vestigated. However,  the  low-frequency  data  superimposed  on  the  high-frequency 
data  severely  limited  the  amount  and  continuity  of  the  high-frequency  data, 
which  are  of  prime  Interest. 
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Figure  4.  Raw  Data  from  DABS  IC 
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Figure  5.  Radar  Velocity  Data  from  DABS  IC 


On  DABS  ID  and  IE  the  Input  to  the  preamplifier  was,  as  previously  described, 
high  passed.  As  can  be  seen  In  figures  6 and  7 the  quantity  of  high-frequency 
data  was  radically  Increased.  There  was  some  low-frequency  data  retained  be- 
cause of  the  extremely  high  gain  of  the  preamplifier  (=  10®).  Figures  8 and  9 
show  the  reduced  data  from  these  two  tests. 

The  velocity  data  from  DABS  IE  (fig.  9)  were  extremely  Interesting  In  that  there 
were  some  obvious  oscillations  on  top  of  the  data.  These  oscillations,  which 
varied  In  period  from  approximately  3 to  5 msec,  were  probably  due  to  the  exci- 
tation of  a transverse  mode  In  the  DABS  IE  tunnel  (the  tunnel  was  13  ft  wide). 
This  would  tend  to  Indicate  that  the  acoustic  velocities  back  In  the  flow  varied 
from  2600  to  4300  ft/sec.  Although  these  results  were  Interesting  and  believ- 
able In  magnitude,  the  precision  of  the  measurement  was  certainly  limited  by  the 
scatter  In  the  data.  That  this  appears  only  on  the  DABS  IE  data  Is  due  to  geo- 
metrical reasons.  The  DABS  IE  tunnel  was  the  largest  tunnel,  and  since  It  was 
the  only  one  with  a sloped  roof  design,  there  was  only  one  direction  (parallel 
to  the  floor)  that  a transverse  mode  could  be  excited.  The  oblate  DABS  lA  was 
certainly  not  condusive  to  this  effect,  nor  were  the  DABS  IC  and  ID,  since  they 
were  constructed  with  flat  roofs.  However,  In  these  latter  cases,  there  exis- 
ted the  possibility  of  exciting  two  different  transverse  modes  of  different 
frequency  and  phase. 

Based  on  the  velocities  Inferred  from  the  radar  data,  the  radar  beam  tracked 
the  Interface  between  the  detonation  products  and  the  clean  air  shock  or  the 
debris  right  behind  It.  Another  possibility  that  should  be  recognized  Is  that 
the  use  of  oxygen-compensated  explosives  such  as  ANFO  or  slurry  may  not  be  ad- 
vantageous to  the  radar  system.  The  absence  of  a burning  Interface  between 
the  detonation  products  and  clean  air  shock  may  allow  nearly  complete  transmis- 
sion of  the  beam.  This  will  have  to  be  tested. 


PARTICLE  VELOCIMETER 

The  secondary  result  of  the  DABS  IC  experiment,  viz.,  that  the  beam  could  see 
through  the  detonation  products,  was  of  prime  Importance  In  the  possible  devel- 
opment of  an  Eulerlan  particle  velocimeter  for  DABS.  Figure  10  shows  the  lay- 
out of  the  proposed  system.  It  Is  the  exact  analog  to  a LASER  velocimeter. 
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Figure  C.  Raw  Data  from  DABS  ID 
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Figure  7.  Riw  Dita  from  DABS  IE 
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Figure  8.  Radar  Velocity  Data  from  DABS  ID 


Shock  Velocity 
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Since  It  was  believed  that  the  extremely  dirty  atmosphere  in  DABS  would  not 
adequately  transmit  optical  radiation,  the  particle  velocimeter  was  chosen. 
However,  when  the  debris  particles  are  small  compared  to  the  wavelength,  as 
in  the  case  of  microwave  radiation,  the  beam  can  be  transmitted. 

The  particle  velocimeter  system  was  originally  proposed  by  the  author  to  the 
personnel  at  the  Radar  Target  Scattering  Facility,  Holloman  Air  Force  Base, 

New  Mexico.  Their  evaluation  of  the  proposal  was  that  it  might  work  with 
their  existing  100-GHz  unit.  The  reason  for  this  choice  of  the  millimeter 
wave  band  was  that  there  are  fewer  problems  at  the  higher  frequency  associated 
with  transmission  through  hot  gases.  The  reason  it  was  not  fielded  on  the  DABS 
series  was  that,  if  it  did  work,  the  cost  of  procurement  and  construction  of  a 
complete  unit  would  be  prohibitive  (more  than  $350,000).  However,  the  fact 
that  the  10.5-GHz  beam  propagated  through  the  plasma  opened  the  door  to  the  use 
of  frequencies  in  this  range  and  thus  the  costs  have  been  significantly  lowered. 
As  a result,  the  University  of  Texas  Applied  Research  Laboratory  was  funded  to 
do  a feasibility  and  cost  study  of  the  proposed  system.  The  results  of  this 
study  are  presented  in  appendix  B.  From  the  study  it  was  recommended  that  a 
16.2-6HZ  source  be  used  in  a geometry  similar  to  that  proposed.  Theoretically, 
this  system  is  sound  and  if  it  is  financially  feasible,  it  should  be  construc- 
ted and  tested  for  DABS. 


SECTION  3 
INFRARED  DETECTOR 


Some  photographic  experiments  were  performed  on  DABS  lA.  High-speed  motion- 
picture  cameras  photographically  recorded  the  motion  of  an  MX  type  pool  as  the 
shock  wave  passed  over  the  water.  Thick  glass  windows  were  used  beside  the 
pool  and  a light  was  shone  through  the  pool  for  Illumination. 

The  results  of  these  experiments  were  quite  confusing.  Shortly  after  the  shock 
arrival,  the  atmosphere  above  the  pool  turned  completely  black,  I.e.,  the  atmos 
phere  appeared  opaque  to  visible  radiation.  This  led  to  some  speculation  about 
the  possibility  of  treating  the  atmosphere  as  a black  body.  If  this  could  be 
done,  the  temperature  could  be  determined  from  the  Infrared  spectrum  emitted 
from  the  shocked  gases.  This  can  be  shown  by  Planck's  Black-Body  Distribution 
Function,  viz. , 


where  I(X,T)  Is  the  Intensity  as  a function  of  the  wavelength,  X,  and  absolute 
temperature  T;  C Is  the  speed  of  light;  h Is  Planck's  constant  (6.6252  x lO*^** 
Joule-sec)  and  k Is  oltzmann's  constant.  Figure  11  shows  a plot  of  this  func- 
tion for  different  temperatures.  Thus,  If  narrow  bandpass  filters  are  used  on 
a number  of  different  wavelengths.  It  can  be  determined  If  the  spectrum  Is  that 
of  a black  body,  and  If  It  Is  the  temperature  Is  uniquely  defined. 

The  first  experiments  were  directed  toward  developing  the  gage  package  for  DABS. 
Since  the  Interface  between  the  detonation  products  and  the  clean  air  shock  Is 
burning,  a plain  Infrared  detector  should  be  able  to  detect  Its  arrival  time. 

The  Information  gained  from  this  type  of  experiment  would  not  only  provide  a 
functional  check  of  the  Infrared  system  but  also  a cross  check  with  the  radar 
data  to  determine  whether  or  not  the  detector  was  really  tracking  the  Interface. 
This  approach  was  taken  first  primarily  because  of  the  Increased  complexity  and 
cost  Inherent  In  making  a five-channel  Instrument  with  narrow  bandpass  filters. 
This  latter  experiment  would  be  undertaken  only  If  the  first  was  successful. 
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Figure  11.  Spectral  Distribution  of  Black-Body  Radiation 
at  Various  Temperatures 


The  first  event  on  which  an  Infrared  detector  was  used  was  DABS  ID.  This  de- 
tector was  a lead  selenide  photoresistor  (Model  No.  4007)  manufactured  by 
Infrared  Industries.  It  has  a l-usec  time  constant  and  a response  as  shown 
In  figure  12.  This  detector  was  chosen  because  of  Its  sensitivity  In  the  ap- 
propriate range  of  wavelengths  (fig.  11),  and  since  the  van  recording  capabil- 
ity was  limited  to  20  kHz,  the  l-psec  time  constant  was  considered  adequate. 

The  conditioning  circuitry  Is  shown  In  figure  13.  Ac -coupling  Into  the  opera- 
tional amplifier  was  necessary  In  order  to  block  the  dc  offset  of  the  detector 
when  the  power  source  (500  Vdc)  was  turned  on. 

The  mechanical  design  for  the  mount  Is  shown  In  figure  14.  It  consisted  of  a 
2-1n-d1ameter  (Inside)  pipe,  4 ft  In  length.  The  bottom  of  the  pipe  was  covered 
and  It  had  a l/4-1n-w1de  slit  through  the  mounting  flange.  A second  silt  was 
used  to  limit  the  field  of  view  of  the  detector.  These  slits  were  oriented  In 
the  DABS  tunnel  so  that  they  would  be  parallel  to  the  shock  front.  The  detec- 
tor, suspended  by  three  threaded  rods  which  were  mechanically  coupled  to  the 
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Figure  12.  Infrared  Detector  Response 


lOO-kQ 


Figure  13.  frared  Detector  Circuitry  on  DABS  ID 
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Pencil  Lead  Mount 


Figure  14.  Mount  for  Infrared  Detector 


pipe  through  pencil  leads,  was  dropped  down  2 ft  Into  the  pipe.  The  purpose  of 
this  mount  was  to  ensure  that  when  the  shock  wave  passed  under  the  detector 
and  propagated  up  the  pipe  walls,  the  detector  would  free  fall  and  thus  be 
Isolated  from  the  shock.  A quartz  window  wus  added  for  further  protection 
of  the  detector.  The  calculated  fill  time  for  the  cavity  produced  In  front 
of  the  detector  was  20  msec  and  thus  the  detector  could  be  expected  to  func- 
tion correctly  for  this  length  of  time.  However,  when  the  oxygen-poor 
detonation  products  enter  this  cavity  they  will  begin  to  burn  and  cause  an 
Illumination  visible  to  the  detector.  This  problem  can  be  alleviated  simply 
by  locating  the  quartz  window  at  the  bottom  of  the  pipe.  Although  this  prob- 
lem was  recognized  It  was  not  acted  upon  simply  because  of  procurement  time 
and  not  knowing  when  the  quartz  window  would  arrive.  However,  because  of 
the  design  used  the  expected  signal  would  be  a montonically  Increasing  func- 
tion until  after  shock  arrival  because  of  secondary  scattering  of  the  mole- 
cules In  the  air  (primarily  CO^).  When  the  detonation  products  Initially 
arrive,  a saddle  point  or  peak  should  be  achieved.  After  the  burning  Inter- 
face between  the  detonation  products  and  the  clean  air  shock  passes,  the  un- 
burned detonation  products  will  begin  to  burn  In  the  cavity  beneath  the  de- 
tector. 

Before  Insertion  on  DABS  ID,  the  detectors  were  checked  with  a Sandia  Labora- 
tories black-body  source  In  conjunction  with  an  Iris  and  camera  shutter.  The 
Iris  was  set  so  that  the  detector  would  have  the  same  field  of  view  as  It 
would  In  DABS  and  the  camera  shutter  was  used  to  give  the  fast  rise-time  pulse 
necessary  because  of  the  ac-coupling  (15  Hz  high  pass).  The  signal,  observed 
on  an  oscilloscope,  was  approximately  10  V for  all  the  detectors  when  the 
black-body  temperature  was  2000”C.  This  agreed  with  detector  sensitivity  cal- 
culations and  thus  the  detectors  were  Installed. 

Four  detector  packages  were  Installed  on  DABS  ID.  These  packages  were  In- 
stalled In  the  roof  at  25.21,  48.25,  75.21,  and  99.33  ft  from  the  explosive- 
driver  end  of  the  DABS  facility.  Two  of  the  packages  failed  before  the  test. 

It  was  presumed  that  the  Input  was  overloaded.  This  caused  a failure  In  the 
operational  amplifier  power  supply,  which  resulted  In  a large  dc-offset  sig- 
nal. However,  since  this  occurred  during  the  countdown,  there  was  no  effort 
made  to  repair  the  detectors. 
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The  Infrared  data  are  shown  In  figure  15.  Figure  15a  shows  a peak  at  shock 
arrival  (determined  from  time-of-arrival  measurements).  This  was  probably 
caused  by  the  fact  that  at  close-in  stations  the  shock  front  Itself  Is  quite 
hot.  The  next  large  peak  (0.9  msec  later)  was  probably  the  arrival  of  the 
detonation  products.  However,  the  existence  of  a peak  and  saddle  point  be- 
tween 4.0  and  4.7  msec  makes  this  doubtful.  From  strictly  shock-propagation 
considerations,  It  seems  that  the  first  large  peak  (0.9  msec  after  shock  ar- 
rival) Is  the  arrival  time  of  the  detonation  products.  The  second  trace 
(fig.  15b)  Is  much  more  classic  In  nature.  At  shock  arrival  the  signal  be- 
gan to  Increase  In  magnitude  to  a saddle  point  and  then  Increased  further 
after  the  detonation  products  commenced  to  burn  In  the  detector  cavity.  The 
uncertainty  In  the  arrival-time  measurement  at  the  first  station  Indicates 
that  the  average  velocity  of  the  Interface  between  the  detonation  products 
and  the  clean  air  shock  was  In  the  band  of  4800  to  5100  ft/sec  when  It  was 
propagating  between  the  25-  and  99-ft  stations.  This  reason  for  the  afore- 
mentioned uncertainty  should  give  a heavier  weight  to  the  lower  velocity. 
These  results--the  measured  shock  velocity,  the  radar  velocity  output,  and 
the  Infrared  band  of  velocities— are  shown  In  figure  16.  It  can  be  noticed 
that  the  lower  limit  of  the  Infrared  velocity  band  agrees  quite  well  with 
that  expected  for  an  average  of  the  radar  velocity  measurements  during  the 
same  time  frame. 

Channel  limitations  made  It  necessary  that  the  number  of  Infrared  experiments 
be  limited  to  two  on  DABS  IE.  The  previously  mentioned  electrical  problems 
prompted  a change  In  the  conditioning  circuit  for  these  tests  (fig.  17). 

The  mechanical  configuration  was  Identical  with  the  exception  that  a 1 -In- 
diameter quartz  window  replaced  the  slot  at  the  bottom  of  the  pipe.  The  par- 
ticular problem  with  this  experiment  was  that  the  roof  on  DABS  IE  was  sloped 
at  26  deg.  That  presented  the  options  of  either  sloping  the  window  with  re- 
spect to  the  axis  of  the  mount  and  accepting  the  lower  transmission  or  allow- 
ing the  pencil  lead  mounted  detector  to  rest  against  the  side  of  the  pipe  and 
accepting  the  Increased  vulnerability  to  shock  loading.  The  latter  option 
was  chosen  and  the  detectors  failed  shortly  after  shock  arrival. 
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Figure  17.  Infrared  Detector  Circuit  for  DABS  IE 


The  Infrared  detector  needs  further  check  out  before  a multichannel  tempera- 
ture detector  Is  made.  The  problems  mentioned,  however,  are  trivial  enough 
that  there  Is  a good  prospect  for  success  with  the  detector  package.  There 
Is  some  question  as  to  whether  or  not  the  detonation  products  really  behave 
as  a black  body.  This  can  be  answered  only  by  using  a temperature  gage  and 
determining  whether  or  not  the  spectrum  Is  equivalent  to  that  of  a black  body 
This  will  be  tested  In  the  follow-on  effort  (DABS  Instrumentation  Development 
Phase  II). 


SECTION  4 
OPTICAL  PHOTODIODE 


Performing  optical  experiments  in  DABS  was  quite  intriguing.  The  possibility 
of  determining  temperature,  particle  velocity,  and  acoustic  velocity  certainly 
merited  investigation.  Specific  interests  were  directed  toward  making  tempera- 
ture measurements  by  the  Rotational  Raman  Scattering  Technique.  The  proposed 
idea  was  to  use  a pulsed  ruby  laser  system  as  a source  and  measure  the  fre- 
quency shifts  in  the  rotational  Raman  lines  as  a direct  inference  of  tempera- 
ture. However,  doing  this  was  questionable  because  of  the  complete  or  nearly 
complete  opacity  behind  the  shock  front.  (See  section  3.)  Therefore,  the 
purpose  of  this  experiment  was  to  determine  if  the  atmosphere  behind  the  shock 
front  was  really  opaque  to  visible  radiation. 

A 650-U  quartz  iodide  photographic  lamp  was  placed  on  one  side  of  the  DABS  ID 
tunnel  behind  a 2-in-thick  glass  window.  The  lamp,  used  in  the  spotlight  mode, 
was  aimed  across  the  tunnel  at  a photodiode  which  was  also  placed  behind  a 2- 
in-thick  glass  window.  Both  items  were  shock  mounted.  The  photodiode  was  an 
RCA  Model  C30822  N-Type  Silicon  Photodiode.  The  spectral  response  of  this  de- 
vice is  shown  in  figure  18.  The  response  in  the  infrared  range  is  extremely 
poor,  and  thus  the  signal  from  the  photodiode  was  dependent  only  upon  the  vis- 
ible light  present.  Since  the  cavity  in  which  the  lamp  was  placed  was  small, 
the  lamp  was  turned  on  by  a relay  switch  activated  from  the  timing  and  firing 
van  associated  with  DABS  ID.  The  lamp  was  turned  on  10  sec  before  detonation 
time.  This  prevented  overheating  and  subsequent  destruction  of  the  lamp  be- 
fore the  actual  test.  Power  for  the  lamp  was  supplied  by  a 5-kW  portable  gen- 
erator. Since  the  lamp,  specifically  the  filament,  was  quite  susceptable  to 
shock  damage,  the  current  to  the  lamp  was  monitored. 

The  electrical  circuit  for  the  photodiode  is  shown  in  figure  19.  The  rise  time 
of  this  circuit  was  tested  in  the  laboratory  with  a strobe  source  and  found  to 
be  less  than  100  nsec;  the  published  rise-time  is  7 nsec.  Since  a fast  enough 
light  source  was  not  available  to  verify  this,  the  published  value  was  assumed 
to  be  correct.  The  signal  cable,  an  RG331 , was  terminated  at  50  A Bell  & 
Howell  3700  B Tape  Recorder  operated  in  the  Wideband  II  FM  mode  (500-kHz 
response)  was  used  to  record  the  data. 
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Figure  19.  Photodiode  Circuit 
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The  light  source  for  the  potential  rotational  Raman  scattering  experiment  was 
a pulsed  ruby  system  that  emitted  a series  of  very  short  (20  nsec)  high-powered 
light  pulses.  The  argument  can  be  made  that  if  debris  is  the  principal  cause 
of  optical  absorption,  there  exists  a higher  probability  of  apparent  transmis- 
sion of  a short  pulse  with  a high-frequency  response  recording  system  than 
there  is  of  propagating  a continuous  beam  with  a comparatively  low-frequency 
response  system.  There  are  two  reasons  for  this,  first,  much  higher  light 
intensities  can  be  obtained  with  pulsed  systems  than  with  continuous  wave  sys- 
tems. Second,  the  short-length  pulse  propagating  at  the  speed  of  light  has  a 
much  lower  probability  of  interacting  with  the  discrete  scattering  centers 
(debris)  than  does  the  continuous  wave.  For  these  reasons  along  with  the  de- 
sire to  determine  if  a continuous  wave  could  propagate,  a high-frequency 
recording  system  was  used  in  addition  to  the  tape  recorder.  This  was  the  par- 
ticular reason  for  using  the  RG351  cable  (foamflex).  In  addition  a Tektronix 
555  delay  triggered  oscilloscope  swept  at  100  nsec/cm  was  used  to  record  the 
data. 

The  data  obtained  from  this  experiment  are  shown  in  figures  20  and  21.  The 
current  monitor  output  (fig.  20)  indicates  that  the  lamp  ceased  to  function 
12.2  msec  after  detonation.  The  arrival-time  data  for  DABS  ID  indicate  that 
the  shock  front  arrived  at  the  lamp  and  photodiode  locations  7.4  msec  after 
detonation  or  the  lamp  failed  4.8  msec  after  shock  arrival.  Figure  21  shows 
that  the  optical  photodiode  continued  to  function  even  after  the  lamp  failed. 
Thus  there  was  some  source  of  light  other  than  the  lamp.  In  fact  almost  im- 
mediately after  detonation,  a pulse  was  observed,  after  which  the  light  inten- 
sity increased  until  it  reached  tape  recorder  bandedge  (twice  the  intensity 
provided  by  the  lamp).  The  light  signal  remained  at  bandedge  until  12.6  msec 
after  detonation  and  then  proceeded  to  zero  or  near  zero  intensity  in  a period 
of  3 to  4 msec.  At  later  times  the  light  intensity  increased  again.  This 
could  have  been  due  to  venting  into  the  cavity  which  incapsulated  the  photo- 
diode. During  the  period  of  time  that  the  signal  was  bandedged,  the  oscillo- 
scope was  triggered  and  as  expected  the  trace  was  off  screen. 

A reasonable  hypothesis  as  to  what  occurred  is  that  right  after  detonation, 
a bright  flash  illuminated  the  tunnel.  (This  is  the  spike  on  the  photo- 
diode trace.)  After  the  burning  interface  between  the  explosive  products 
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Figure  20.  Current  Monitor  Output 


Figure  21 . Pfntodlode  Output  froM  DABS  ID 


and  the  clean  air  shock  passed  the  detector,  the  signal  magnitude  decreased  to 
a level  below  bandedge.  This  was  within  3 msec  after  the  burning  Interface 
passed  the  photodiode.  As  this  Interface  proceeded  to  propagate  down  the  tun- 
nel, the  light  Intensity  (Incident  upon  the  photodiode)  due  to  this  burning  In- 
terface decreased  rapidly.  Thus,  since  the  Intensity  decreased  over  a period 
of  3 to  7 msec,  the  detonation  products  did  transmit  light.  However,  the  lati?- 
tlme  transmission  of  this  light  was  somewhat  minimal. 

The  results  of  this  experiment  are  somewhat  conflicting  with  those  of  the  photo 
graphic  experiment  on  DABS  lA.  The  DABS  lA  results  Indicated  opacity;  the  DABS 
ID  photodiode  Indicated  transmission.  This  discrepancy  and  cost  considerations 
lead  to  the  conclusion  that  the  Rotational  Raman  experiment  would  not  be  at- 
tempted at  this  time. 
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SECTION  5 
DEBRIS  COLLECTOR 


The  stagnation  pressure  data  generated  in  DABS  show  an  anomalous  second  peak. 
One  of  the  postulated  reasons  for  this  peak  is  that  pressure  gages  were  af- 
fected by  debris  impact.  For  this  reason  a debris  collector  was  fielded  on 
DABS  ID.  It  consisted  of  a simple  pipe  mounted  on  top  of  the  second  post. 

The  inside  diameter  of  this  pipe  was  1.5  in  and  the  total  mass  of  debris  col- 
lected in  the  pipe  was  36.78  g.  This  yielded  an  areal  density  of  0.0458 
Ib/in^.  This  debris  was  subjected  to  a sieve  analysis  and  the  results  are 
shown  in  table  1.  The  prime  fault  with  the  experiment  was  that  there  was  no 
way  to  separate  the  debris  collected  from  the  flow  from  that  collected  as  a 
result  of  fallback. 

If  some  assumptions  about  the  average  velocity  at  which  the  debris  was  moving 
are  made,  the  impulse  can  be  determined  from  the  assumption  that  the  momentum 
of  the  particles  impacting  the  stagnation-pressure  gages  is  converted  to  an 
effective  pressure  reading  at  the  gages.  Thus  if  7 is  the  average  velocity 
and  m is  the  areal  density  of  the  debris  mass,  the  specific  debris  impulse, 
Iq,  is  mV. 


Table  1.  DABS  ID  Debris  Distribution 


Sieve 

No. 

Particle 

Size, 

in 

Mass, 

g 

> 20 

0.0331 

9.2 

40 

0.0165 

9.3 

60 

0.0098 

7.0 

140 

0.0041 

5.8 

200 

0.0029 

2.0 

< 200 

0.0029 

3.5 

43 


The  next  assumption  made  was  that  the  bulk  of  the  debris  was  picked  up  right 
after  the  detonation,  where  the  turbulence  Is  great.  Thus,  It  can  be  expected 
that  the  bul  if  the  debris  was  generated  In  the  area  close  to  the  detonation 
chamber.  If  this  was  the  case,  the  debris  was  mostly  located  near  and  behind 
the  detonation-products  Interface.  The  presumption  was  that  the  radar  unit 
tracked  the  Interface  Itself  or  the  debris  near  It  (fig.  16).  The  Infrared 
data  show  the  range  of  average  velocities  for  the  Interface  during  the  time 
period  noted  on  the  graph.  The  agreement  Is  fairly  good  and  at  the  times  of 
Interest  the  peak  velocity  of  the  Interface  varied  from  5100  to  3600  ft/sec. 
Since  the  particle  velocity  falls  off  rapidly  behind  the  Interface,  the  aver- 
age velocity  of  the  debris  In  this  region  Is  certainly  lower  than  the  peak  ve- 
locity. Thus,  If  the  presumption  that  the  bulk  of  the  debris  Is  near  to  or 
behind  the  Interface  Is  correct,  average  debris  velocities  In  the  range  of 
2000  to  4000  ft/ sec  could  be  expected. 


The  next  question  to  answer  Is,  "How  much  of  the  debris  came  from  the  fallback 
as  opposed  to  the  shock?"  It  Is  believed  that  the  large  particles  (greater 
than  0.0331  In  In  diameter)  came  from  the  fallback  or  came  so  late  and  so  slow 
that  they  did  not  appreciably  add  to  the  Impulses.  This  presumption  leads  to 
an  areal  density  of  0.0344  lb/1n^ 


The  last  assumption  made  was  that  the  second  peak  on  the  stagnation-pressure 
measurements  was  totally  due  to  the  debris  In  the  flow.  If  this  was  the  case, 
an  average  velocity  for  the  debris  can  be  calculated  from  the  Impulse  due  to 
the  second  peak.  This  was  accomplished  by  simply  smoothing  the  entire  curve 
and  taking  the  area  under  the  second  peak  (fig.  22).  In  both  these  cases  the 
specific  Impulse  was  approximately  3.25  psi-sec.  Thus 


or 


mV  = V = 3.25 

V - 3000  ft/ sec 


This  Is  a reasonable  value.  Thus,  as  much  as  one  Is  able  to  generalize  with 
one  datum  point,  the  second  peak  can  be  plausibly  explained  with  a debris- 
impacting model.  Obviously,  however,  the  second  pulse  on  the  stagnation  meas- 
urements was  not  totally  due  to  debris;  It  Is  more  likely  that  It  was  due  to  a 
combination  of  the  higher  density  detonation  products  and  the  debris. 
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Figure  22.  Second  Peak  In  DABS  ID  (1  of  2) 
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Figure  22.  Second  Peak  in  DABS  ID  (2  of  2) 
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The  prime  difficulty  with  the  experiment  was  the  Inability  to  correct  for  the 
fallback.  Thus,  the  next  obvious  step  In  the  development  was  to  separate  the 
debris  In  the  flow  from  the  fallback  debris.  Figure  23  shows  the  mechanism 
used  on  DABS  IE.  This  device  consisted  of  a pipe  (the  same  type  as  that  used 
on  DABS  ID)  mounted  on  a wing  with  a closing  mechanism  designed  such  that  the 
debris  In  the  rear  of  the  pipe  would  be  separated  from  the  fallback  debris  In 
the  front  of  the  collector.  The  closing  mechanism  was  a spring-driven  door. 

The  spring  provided  a force  of  approximately  20  lb  on  the  door.  The  trigger- 
ing mechanism  was  simply  a plate  with  a tongue  that  was  Inserted  Into  a slot 
In  the  door  when  the  device  was  In  the  cocked  mode.  Cocking  was  accomplished 
by  Inserting  a screw  through  a threaded  hole  In  the  top  of  the  pipe.  When 
the  shock  wave  reaches  the  device,  the  flow  Is  stagnated  by  the  trigger  plate 
and  the  resultant  force  on  the  plate  (between  15.000  and  22.000  1b)  simply 
blows  the  trigger  plate  downstream;  this  releases  the  door  and  allows  It  to 
close. 

The  triggering  mechanism,  however,  caused  some  concern  because  of  the  differ- 
ential pressures  created  along  the  door.  These  differential  pressures  resulted 
from  the  hole  In  the  side  of  the  wing  for  the  trigger  plate.  It  caused  the 
pressure  there  to  be  approximately  the  side-on  pressure  (nominally  600  to  700 
psi).  while  the  pressure  Inside  the  pipe  was  the  stagnation  pressure  (nominally 
a few  thousand  pounds  per  square  Inch).  Thus  there  was  the  strong  possibility 
that  the  debris  contained  In  the  flow  passing  the  door  could  Jam  the  door.  For 
this  reason,  other  concepts  were  considered.  These  Included  electrical-shutting, 
explosive-shutting,  and  pencil-lead  triggering  techniques.  However,  all  these 
were  rejected  because  of  their  complexity  or  poor  reliability.  With  this  ex- 
pressed reservation.  It  was  felt  that  the  chosen  option  was  the  simplest,  least 
expensive,  and  most  reliable  system. 

An  Important  factor  considered  In  the  design  was  the  time  It  would  take  for 
the  door  to  close.  The  door  would  not  start  to  close  until  after  the  shock 
front  reached  the  device  and  the  differential  loading  on  the  door,  resulting 
from  the  flow  mechanism,  became  less  than  20  lb.  Thus,  the  door  would  close 
the  fastest  with  no  differential  loading  on  It.  This  time  was  measured 
through  high-speed  photography;  the  results  are  shown  In  figure  24.  These 
data  show  that  the  debris  collector  was  halfway  closed  16  msec  after  the  shock 
arrival;  for  the  DABS  waveforms,  this  was  deemed  adequate. 
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The  system  on  DABS  IE  did  not  work  as  planned.  The  door  closed  only  halfway; 
It  was  jammed  by  a small  amount  of  debris.  There  was  no  noticeable  debris  In 
the  spring  cavity  below  the  door.  However,  this  problem  was  not  serious  since 
the  overburden  fallback  was  minimal;  on  DABS  ID  the  fallback  literally  covered 
the  debris  collector.  The  design  of  DABS  IE  with  Its  sloped  roof  minimized 
this  problem.  The  bulk  of  the  overburden  fell  back  outside  the  tunnel.  The 
collector  was  examined  Immediately  after  the  test  and  no  noticeable  debris  was 
found  In  the  front  of  the  collector;  all  of  the  debris  was  behind  the  half- 
closed  door.  The  total  debris  collected  was  15.66  g.  The  sieve  analysis  Is 
shown  In  table  2. 


An* areal  density  of  0.0195  Ib/ln*  can  be  Inferred  from  the  total  weight  of  the 
debris.  Figure  25  shows  the  output  from  one  of  the  stagnation  gages.  This  gage 
was  located  approximately  the  same  linear  distance  (75  ft)  from  the  driver  end 
as  the  debris  collector.  The  apparent  Impulse  due  to  the  second  peak  was  ap- 
proximately 3 psi-sec.  If  It  Is  presumed  that  this  was  entirely  due  to  debris. 


or 


mV  ■ 


0.0195  TT 


= 3 


V * 4923  ft/sec 


Table  2.  DABS  IE  Debris  Distribution 


Sieve 

No. 

Particle 

Size, 

In 

Mass, 

9 

> 20 

G.0331 

4.99 

40 

b.0165 

3.59 

60 

0.0098 

2.36 

140 

0.0041 

1.91 

200 

0.0029 

0.73 

< 200 

0.0029 

1.78 
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At  the  time  the  second  pulse  occurred  (13.5  msec  after  detonation),  the  velo- 
city of  the  Interface  between  the  detonation  products  and  the  clean  air  shock 
was  5500  ft/sec  (fig.  9).  Thus  the  average  debris  velocity  of  4923  ft/sec 
Is  probably  high  and  It  may  be  concluded  that  the  second  pulse  Is  not  entirely 
due  to  debris  Impact  but  a significant  portion  of  the  Impulse  In  this  case  Is 
most  likely  due  to  debris. 

Visual  Inspection  of  the  debris  showed  that  It  consisted  of  carbon,  concrete, 
sand,  and  wood.  The  bulk  of  the  carbon  probably  came  from  the  decomposition 
of  the  plastic  Jacket  surrounding  the  detonating  cord. 

It  seems  that  the  basic  concept  of  a debris  collector  Is  quite  valid  to  support 
the  analysis  of  the  data  from  DABS.  However,  the  design  needs  some  Improvement 
or  possibly  some  radical  change.  An  obvious  Improvement  would  be  to  modify 
the  door-closing  technique  used  to  separate  the  debris  from  overburden  fall- 
back. In  addition,  the  stagnation  of  the  flow  In  a plugged  pipe  Indicates 
that  there  Is  a possible  back  flux  of  debris  out  of  the  collector  as  the  out- 
side pressure  decreases.  Admittedly,  this  effect  should  be  small  since  all 

the  debris  goes  to  the  end  of  the  collector  and  Is  deposited  on  the  plug.  How- 
ever, since  there  Is  the  possibility  of  putting  the  collector  tube  Into  reson- 
ance during  unloading  and  sweeping  out  the  debris  Is  finite,  some  modification 

may  be  desired.  Additionally,  the  present  debris  collector  gives  no  Infortna- 

tlon  regarding  the  time  history  of  the  debris  flow.  Debris  collectors  to  do 
this  have  been  designed  but  little  success  has  been  attained.  Primarily  they 
operate  with  the  timed  detonation  of  shaped  charges  around  both  ends  of  a tube. 
It  Is  apparent  that  these  designs  can  be  Improved  and  this  should  be  undertaken. 
Another  device  that  may  be  of  use  because  of  the  high  cross  section  that  debris 
has  to  beta  particles  Is  the  beta  densitometer. 
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SECTION  6 

SHOCK-ON-SHOCK  INSTRUMENTATION 


One  of  the  flow  variables  of  primary  importance  in  DABS  is  the  mach  number  i" 
the  flow.  Mach  number  is  defined  as  the  vatio  of  flow  op  pavtiale  veto,  .j  to 
local  aaoustia  oelooity . The  Stagnation  pressure  is  inferred  from  pitot  tube 
type  measurements  from  DABS.  The  supersonic  pitot  type  measurement  technique 
is  diagrammed  in  figure  26.  With  the  notation  that  is  the  local  mach  number 
in  the  flow,  p^^  the  static  pressure  in  the  flow,  p^  the  stagnation  pressure,  p^ 
the  static  pressure  in  the  stagnated  region,  end  p^^  the  pitot  tube  pressure. 


and 


tt 


1 


(ref.  2) 


The  relationship  between  these  measured  quantities  is  greatly  dependent  upon  a 
knowledge  of  the  local  mach  number.  Thus,  the  experimental  determination  of 
local  mach  number  is  of  primary  importance  in  understanding  the  flow  conditions^ 
in  DABS.  Normally,  a shock-diffraction  experiment  is  performed  and  from  a meas- 
urement of  mach  angle,  the  mach  number  is  inferred.  Unfortunately  DABS  is  the 
worst  of  all  possible  worlds  to  do  this  type  of  experiment  in,  since  the  major 
part  of  the  flow  is  transonic.  Therefore,  an  experiment  suggested  by  R.  Shunk 
of  Electromechanical  Systems  and  B.  Sturtevant  of  the  California  Institute  of 
Technology  was  implemented.  The  theoretical  analysis  of  this  experiment  is 
contained  in  appendix  C. 


The  experimental  layout  is  shown  in  figure  27.  A 2-lb  charge  of  C4  explosive 
was  mounted  in  the  wall  of  DABS  IE.  Three  fast-response  (PCB  quartz)  gages 
were  placed  downstream  and  above  the  explosive  charge.  Piezoelectric  pins 
were  placed  further  downstream  (=  15  ft  past  the  gages).  These  pins  were  used 
to  trigger  the  detonation  system  for  the  explosive  charge.  Thus,  at  the  time 


Shapiro,  A.,  The  Dynarniae  and  Themodynarnias  of  Compressible  Fluid  Flow, 

Vol.  I,  The  Ronald  Press  Company,  N.Y.,  N.  Y.,  1953. 
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of  detonation  of  the  small  charge,  the  shock  wave  produced  would  be  added  to 
the  shock  wave  In  DABS  (shock  on  shock).  As  the  secondary  shock  wave  propa- 
gates In  the  flow  It  will  pass  over  the  three  pressure  transducers.  Since 
these  transducers  were  close  together  ( 4 In  apart)  and  far  from  the  source 
(4  ft),  the  secondary  shock  would  propagate  across  these  transducers  with  a 
constant  direction  with  respect  to  the  DABS  flow.  The  angle  between  this  sec- 
ondary shock  propagation  and  the  DABS  Is  defined  as  6.  (See  appendix  C.)  If  4> 
Is  the  angle  between  a line  connecting  the  explosive  charge  and  the  gages  and 
flow  direction  and  If  6 Is  the  distance  between  the  gages  (as  shown  In  appendix 
D). 


At* 

tane  - 

sine 

At* 


, tan(9  - (<>) 
tan0 


where  c Is  the  secondary  shock  velocity,  u Is  the  DABS  flow  velocity,  At^  Is 
the  time  between  secondary  shock  arrivals  at  gages  1 and  2,  and  At^  Is  the 
time  between  secondary  shock  arrivals  at  gages  2 and  3 (fig.  27).  If  the  mach 
number  of  this  secondary  shock  were  known,  the  experiment  would  simultaneously 
determine  both  the  particle  velocity  and  the  acoustic  velocity  In  the  flow, 
and  since  the  mach  number  Is  the  ratio  of  these  two  quantities.  It  too  could 
be  determined. 


The  particular  disadvantages  of  this  type  of  measurement  are  as  follows: 

(1)  only  the  mach  number  at  one  given  time  back  In  the  flow 
Is  determined, 

(2)  the  pressure  pulses  are  so  noisy  that  discriminating 
the  secondary  shock  from  the  primary  shock  or  other 
noise  Is  difficult,  and 

(3)  geometry  and  time  dictate  that  the  explosive  charge  be 
put  In  a steel  canister  mounted  In  the  same  wall  as  the 
gages.  (This  mounting  tends  to  direct  the  blast  from 
the  small  charge  away  from  the  wall  and  thus  limits  the 
secondary  shock  Impingment  upon  the  gage.) 


SS 


The  experiment  was  a qualitative,  but  not  a quantitative,  success.  Figure  28 
shows  the  data  trace  obtained  from  gage  1.  The  magnitude  of  the  secondary 
shock  Is  quite  small.  Although  this  shock  was  present  on  the  other  traces, 
the  one  on  this  trace  was  the  most  clearly  defined.  This  low  level,  combined 
with  the  noise,  precluded  the  possibility  of  making  the  time-of-arrival 
measurement  within  the  required  precision.  It  Is  expected  that  this  problem 
can  be  overcome  by  directing  the  shock  toward  the  gages  Instead  of  away  from 
them.  This  can  be  accomplished  by  placing  the  charge  In  the  floor  upstream 
from  the  gages.  In  this  manner  the  secondary  shock  would  propagate  more  In 
the  direction  of  the  gages.  The  Importance  of  this  type  of  measurement  man- 
dates the  further  development  of  this  technique  and  appropriate  Instrumenta- 
tion should  be  fielded  on  future  DABS  events. 
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Secondary  Shock  Arrival 
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Figure  28.  Shock-on-Shock  Experimental  Data  (Gage  1 Output) 


SECTION  7 

PRESSURE  AND  VELOCITY  GAGES 


STAGNATION  PRESSURE  GAGE 

Primary  development  of  the  bar  gage  has  already  been  reported  In  reference  3. 
The  only  modification  to  this  gage  was  the  addition  of  a 0.2-In- thick  shield 
of  RTV  silastic  filled  with  coarse  lead-powder  granules.  Although  this  addi- 
tion caused  some  ringing  and  slapping  effects  In  the  reflected  pressure  mode, 
no  effects  were  noted  In  the  dynamic  flow  encountered  In  DABS.  Apparently  this 
shield  decreased  the  posttest  dc-offset  of  the  gage.  This  offset  was  probably 
due  to  debris  Impact  which  caused  permanent  deformation  on  the  front  surface 
of  the  gage.  Obviously  the  shield  absorbed  most  of  the  momentum  of  the  debris 
particles  and  thus  lessened  the  effect. 


DRIVER-CHAMBER  PRESSURE  GAGES 

An  attempt  was  made  on  DABS  lA  to  measure  the  equilibrium  driver  pressure. 

This  pressure  Is  defined  as  the  Initial  driving  pressure  due  to  detonation  of 
the  detonating  cord  distributed  In  a known  volume.  On  DABS  lA  a bar  gage  was 
mounted  surface  flush  In  the  floor  In  the  center  of  the  driver  chamber.  No 
protection  was  used  on  the  front  surface  and  the  gage  failed  before  any  pres- 
sure could  be  determined  from  Its  output.  Another  attempt  was  made  on  DABS  IE. 

The  expected  pressure  can  be  calculated  by 

p • (y  - 1)  f 

where  P Is  the  driver  pressure,  E/V  Is  the  specific  explosive  energy  density 
In  the  driver  chamber,  and  y Is  the  ratio  of  Isobaric  to  Isovolumetric  specific 


3.  Simmons,  Kenneth  B.,  Development  of  Heaoreeiative  Bop  Gage^  AFWL-TR-76-65, 
Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico,  December 
1976. 
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heats.  With  a y of  1.2,  a charge  density  of  6.4  1b/ft’,  and  an  explosive  en- 
ergy of  2 X 10‘  ft-lb/lb,  the  resultant  pressure  Is  approximately  18,000  psi. 
Since  the  gages  were  located  near  the  detonating -cord  array,  peak  pressure 
could  be  much  higher  than  this  because  the  peak  pressure  of  the  shock  wave 
emanating  from  the  explosive  Is  higher  than  the  equilibrium  gas  pressure  In 
the  vicinity  of  the  explosive.  Since  these  detonation  wave  pressures  can  be 
many  times  greater  than  the  desired  measurement,  the  gages  must  either  be  able 
to  measure  them  and  survive  or  a mounting  scheme  must  be  devised  to  ensure  that 
their  effect  will  be  limited.  Two  different  schemes  were  fielded  on  DABS  IE— 
a 30,000-psl  Kullte  gage  with  a MOD  V debris  shield,  and  a ported  bar  gage. 

The  description  of  the  MOD  V debris  shield  experiment  as  given  by  J.  Quintana 
of  AFUL/DED-I  Is  as  follows: 

Tha  debris-shielded  Kulite  gage  sensing  sahme  is  shoum  in  figure  29. 

The  sensing  installation  to  measure  driver  ohamber  pressure  is  shorn 
in  figure  SO.  The  installation  consisted  of  the  regular  AFWL-developed 
hardware  oast  into  the  concrete  side  wall  of  the  driver  ohamber.  In  the 
preshot  configuration  the  Primaoord  bundle  was  approximately  8 am  from 
the  module  sensing  surface. 

The  ported  mount  for  the  bar  gage  Is  shown  In  figure  31.  The  nylon  Insert  was 
used  as  a shock  Isolator  for  the  gage.  The  port  geometry  was  modeled  after  the 
design  of  a B filter.  The  response  of  the  porting  scheme  Is  described  In  ref- 
erence 4.  The  10-to-90  rise  time  of  the  B filter  was  55  ysec.  In  this  case 
the  rise  time  will  be  slightly  longer  because  of  the  slightly  larger  cavity 
below  the  entrance  port  (0.36  In  diameter  as  opposed  to  0.306  In  diameter). 

This  difference  was  necessitated  since  the  original  design  was  for  a Kullte 
gage  with  a smaller  active  gage  area.  This  ported  mount  was  designed  to  act 
as  a low-pass  acoustic  filter  and  thus  protect  the  gage  from  the  high-pressure 
spikes  caused  by  the  detonation  waves  emanating  from  the  Individual  explosive 
charges.  Since  the  rise-time  requirement  for  the  measurement  was  In  the  ranop 
of  100  to  500  ysec,  the  design  was  Judged  adequate.  Additional  protection  iur 
the  gage  was  afforded  by  using  a lead  particle  load  RTV  compound  over  the  front 
surface  of  the  gage.  This  0.2-1n-th1ck  shield  protected  the  gage  from  debris 
Impinging  on  the  gage  through  the  0.125-1n-d1ameter  port. 


4.  Baum,  Neal,  Shook  Tube  Gage  Evaluation  Quick-Look  Data  Reportt  CERF  Letter 
Report  FI-3,  Eric  H.  Wang  Civil  Engineering  Research  Facility,  University 
of  New  Mexico,  Albuquerque,  New  Mexico,  August  1975. 
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Figure  30.  Sensing  Installation 
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7/8-20-UMEF-2A 


One  of  the  questionable  points  concerned  with  the  gage  was  the  weakness  of  the 
nylon  grommet  Isolator.  However,  the  gage  has  a finite  mass  (=  50  g)  and  the 
time  required  to  make  a meaningful  measurement  Is  relatively  short  (100  to  500 
usee).  Thus  with  an  18,000-psl  pressure  applied  to  Its  surface  and  the  strength 
of  the  nylon  threads  neglected,  the  gage  will  move  only  0.03  to  0.75  In  and  the 
space  behind  the  gage  Is  adequate  for  this. 


The  result  of  this  experiment  Is  shown  In  figure  32.  There  was  an  obvious 

catastrophic  failure,  which  was  probably  due  to  cable  failure  since  the  gage 

was  still  functioning  after  removal  from  DABS  IE.  The  atmosphere  was  obviously 

quite  harsh  (fig.  33).  The  front  surface  of  the  gage  mount  was  dimpled  Inward 

(2  to  3 nm).  According  to  J.  Quintana, 

Figure  34  ia  the  data  plot  from  the  debria-ahielded  Kulite  gage  aenaing 
inatallatioH.  Poatteat  inapeation  of  the  inatallation  revealed  the  ef- 
feata  of  the  severe  environment  (fig.  3S>),  One  web  of  the  debria  ahield 
waa  impaated  by  aompetent  debria ^ which  caused  severe  deformation  of  the 
web.  However^  the  debria  ahield  did  effectively  prevent  damage  to  the 
tratiaducer  frotn  the  direct  impact  of  debria  on  the  aenaitive  area  of  the 
unit.  The  geometry  of  the  steel  module  showed  very  little  deformation. 
However,  some  defotmation  occurred  in  the  connector  shell  in  that  the 
module  could  not  be  unplugged  normally  and  removed  from  the  body.  A 
post-mortem  diaacction  of  the  transducer  revealed  that  a portion  (shaped 
like  a truncated  cone)  of  the  silicon  disk  sensing  clement  was  dislodged 
from  the  backside,  and  this  deatiroyed  the  electrical  atrain-aenaing 
elements  located  in  the  active  area  of  the  transducer.  Microscopic  ex- 
amination of  the  cavity  created  in  the  backside  indicated  crystal  line 
shear  failure  from  maaaive  overatreaaing  of  the  element  due  to  over- 
preaaure.  Electrical  continuity  of  the  land-line  cable  to  the  sensing 
location  waa  checked,  and  all  conductors  were  intact  and  operatioml. 

Thus,  further  development  of  the  driver-pressure  measurement  technique  Is  ne- 
cessary. This  will  probably  Involve  the  use  of  much  more  rugged  gages  (ytter- 
bium or  carbon). 


MAGNETOHYDRODYNAMIC  VELOCITY  GAGE 

The  magnetohydrodynamic  technique  of  generating  power  has  long  been  known  and 
has  recently  come  Into  vogue  as  one  of  the  possible  alternatives  for  the  pro- 
duction of  electrical  power.  The  use  of  this  principle  to  design  a flow  or 
particle-velocity  gage  appears  attractive,  since  a self-generating  measurement 
Is  Involved. 
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Figure  33.  Bar  Gage  Mount  After  Test 


The  principle  Is  briefly  explained  In  figure  36.  If  two  parallel  conducting 
plates  (parallel  to  the  xz-p1ane)  are  oriented  so  that  the  flowing  gases  In  a 
shock  wave  can  pass  between  them,  the  flow  velocity  can  be  determined  through 
the  utilization  of  a magnetic  field.  If  F Is  the  force  on  a charged  particle 
In  the  flow,  V the  velocity  of  that  particle,  n^  the  number  of  particles  with 
charge,  B the  magnetic  field,  and  E the  electric  field,  with  Einstein  nota- 
tion the  Lorentz  force  equation  can  be  written 

F ■ n^q^(E  + V x B) 

Now  If  It  Is  presumed  that  there  Is  some  measurement  resistance  (gage  term- 
ination resistance)  between  the  plates  that  Is  much  greater  than  that  caused 
by  the  resistance  of  the  hot  gases  In  the  flow,  the  current  flow  between 
the  plates  will  be  essentially  zero  or 

or 


I 

I 

I 
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Figure  36.  Principle  of  Magnetohydrodynamlc  Velocity  Gage 
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If  V Is  the  voltage  between  the  two  plates  and  l Is  the  distance  between  the 
two  plates. 

Therefore,  If  the  distance  between  the  plates  and  the  magnetic  field  are  known, 
a reading  of  the  voltage  between  the  plates  becomes  a direct  measure  of  the 
flow  velocity.  This  measurement  Is  completely  Independent  of  the  charges  In 
the  flow  In  so  far  as  the  conductivity  of  the  flow  Is  high  with  respect  to  the 
measurement  resistance. 

The  promise  of  this  type  of  system  was  obvious  and  Artec  Associates,  Inc.  of 
Hayward,  California  was  funded  to  pursue  Its  development.  The  report  on  this 
effort  Is  contained  In  appendix  D.  In  conjunction  with  this  effort  an  experi- 
ment was  conducted  on  DABS  IB.  The  purpose  of  which  was  to  qualitatively  de- 
tennlne  If  the  conductivity  within  a DABS  flow  was  high  enough  for  a magneto- 
hydrodynamic experiment. 

The  experiment  consisted  of  placing  two  1-ft^  copper  plates  1 ft  apart  In  the 
floor  of  DABS.  One  surface  of  each  plate  was  exposed  and  the  other  was  Insul- 
ated with  REN  Plastics  RP-3269-1  Epoxy.  The  plates  were  Interconnected  through 
a 1-Msl  resistance.  The  signal  across  this  resistance  was  buffered  through  an 
operational  amplifier  and  a 75-n  output  Impedance  was  used  to  drive  the  signal 
cable. 

The  observed  signal  In  this  experiment  oscillated  between  peak  values  greater 
than  30  V.  It  was  believed  that  static  charge  deposition  on  the  plates  was 
the  reason  for  this  and  thr  some  seeding  would  be  needed.  The  chosen  seed 
material  was  potassium  chloride.  It  was  chosen  over  the  more  optimal  cesium 
chloride  because  of  cost. 

Briefly,  the  result  of  Artec's  experimental  effort  was  that  the  conductivity 
of  the  flow  for  DABS  conditions  Is  not  high  enough  to  obtain  a reasonably 
noise-free  signal  even  If  the  flow  Is  seeded.  Thus  It  was  decided  not  to 
field  these  gages  on  a DABS  experiment.  However,  In  spite  of  the  unfavorable 
signal -to-nolse  ratio.  It  was  obvious  that  the  velocity  signal  was  present. 

With  higher  pressures,  such  as  those  achieved  In  a nuclear  test  or  In  the 
trench  program,  the  performance  will  Improve  and  the  gage  should  be  consid- 
ered for  these  types  of  programs. 
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SECTION  8 

CONCLUSIONS  AND  RECOMMENDATIONS 


Phase  I of  the  DABS  Instrumentation  development  has  come  a long  way  in  the  ex> 
ploratlon  of  new  and  novel  techniques  for  blast  Instrumentation.  Although  suc- 
cesses have  been  significant,  there  remains  areas  In  which  further  development 
Is  needed.  This  will  probably  be  the  case  for  some  time,  since  DABS  measure- 
ments are  made  In  an  extremely  adverse  atmosphere. 

The  use  of  a traffic-control  radar  unit  to  obtain  useful  data  Is  one  of  the 
unique  techniques  that  has  come  out  of  this  program.  Although  one  has  yet  to 
be  built,  the  Eulerlan  microwave  particle  velocimeter  Is  theoretically  and  eco- 
nomically feasible.  The  need  for  particle  velocity  data  from  DABS  and  the  op- 
portunity to  develop  anemometers  that  would  be  useful  for  related  and  totally 
different  applications  make  the  velocimeter  one  of  the  most  promising  devices. 

Complementary  data  obtained  from  Infrared  measurements  of  the  time-of-arrival 
of  the  Interface  between  the  detonation  products  and  the  clean  air  shock  have 
been  useful.  Although  the  technique  needs  some  minor  refinement,  the  obvious 
extension  of  making  temperature  measurements  should  be  pursued.  This  type  of 
temperature  measurement  Involves  the  use  of  a number  of  Infrared  detectors 
with  narrow  band  Infrared  filters.  The  data  obtained  will  be  self-sufficient 
proof  or  disproof  of  the  technique.  If  the  data  show  that  the  spectrum  Is  a 
black-bod^'  spectrum,  the  temperatures  will  be  uniquely  defined;  If  the  data 
do  not  follow  a black-body  spectrum,  the  temperature  will  not  be  defined.  Con- 
tingent upon  the  failure  of  this  method,  some  of  the  more  exotic  Infrared  and 
opMcal  techniques  should  be  considered.  Some  of  these  techniques  are  described 
In  appendix  E.  Also,  additional  consideration  and  development  should  be  given 
to  AFWL's  heat-transfer  gage.  This  might  also  be  extended  to  Intrinsic  thermo- 
couples. 

The  results  from  the  optical  photodiode  were  somewhat  anomalous  In  that  they 
radically  conflicted  with  previous  photographic  data.  However,  these  data  re- 
opened the  possibility  of  doing  optical  experiments  such  as  holography  and  ro- 
tational Raman  measurements.  However,  since  these  experiments  are  extremely 
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expensive  to  field  both  In  equipment  and  labor,  they  should  be  undertaken  only 
after  a low-cost  feasibility  experiment  has  been  conducted.  If  this  feasibil- 
ity experiment  shows  results  that  Indicate  a high  probability  of  success  with 
these  experiments,  then  and  only  then  should  they  be  considered. 

The  debris -col lection  experiments  Indicated  that  debris  Impact  Is  a probable 
contributor  to  the  second  pulse  on  the  pitot  tube  measurements  In  DABS.  How- 
ever, more  work  needs  to  be  done  on  closure  mechanisms  for  the  debris  collec- 
tors. Additional  considerations  should  also  be  given  to  Improvement  of  the 
explosive-closing  debris  collector  that  traps  a segment  of  the  flow.  If  they 
can  be  made  to  work,  they  would  probably  give  a time-history  record  of  the 
size  and  density  of  the  debris.  However,  until  they  are  developed,  debris  col- 
lectors similar  to  those  presently  used  should  be  fielded.  The  upcoming  varia- 
tion of  explosive  drivers  necessitates  this,  since  the  absence  of  the  plastic 
Jacket  containment  of  the  explosive  will  certainly  affect  the  quantity  of  de- 
bris In  the  flow. 

Although  the  shock-on-shock  experiment  was  a qualified  success,  the  configura- 
tion of  the  experiment  should  be  changed  as  recommended  In  the  body  of  this 
report.  This  experiment  should  be  fielded  on  all  possible  DABS. 

Some  advance  must  be  made  If  driver-chamber  pressure  Is  to  be  measured  In  DABS. 
The  attempted  transduction  schemes  were  clearly  Inadequate.  The  most  promising 
devices  In  this  area  are  the  higher-stress  ground-shock  transducers  (carbon  and 
ytterbium).  The  magnetohydrodynamic  velocity  gage  clearly  cannot  be  used  on 
DABS  unless  higher  pressure  levels  are  achieved.  Based  on  experience  with  this 
gage.  It  clearly  functions  well  when  the  conductivity  Is  high  enough  and  may 
have  pertinent  application  on  other  types  of  tests  such  as  underground  nuclear 
tests.  The  bar  gage  Is  essentially  operational  and  can  be  used  In  lieu  of  or 
In  addition  to  normal  commercially  available  gages.  There  still  exists  a def- 
inite need  for  a debris- Independent  pressure  gage.  This  will  have  to  be  accom- 
plished If  the  gas  dynamic  properties  of  the  flow  are  to  be  ascertained.  AFUL 
has  been  working  on  this  problem  with  some  success  and  certainly  this  work 
should  be  continued.  If  this  effort  Is  ultimately  successful,  blunt  cylinder 
drag  measurements  should  be  undertaken.  These  measurements  can  be  used  to  de- 
termine the  mach  number  of  the  transonic  flow  (appendix  F). 
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INTRODUCTION 


The  purpose  of  the  Dynamic  Airblast  Simulator  (DABS)  Is  to  load  defense  struc- 
tures with  airblast  pressure  similar  to  that  expected  from  a nuclear  detona- 
tion. One  of  the  principal  quantities  to  be  measured  Is  the  particle  velocity 
In  the  shock  wave.  It  Is  proposed  that  a thin  aluminum-foil  target  In  the 
shock  wave  be  tracked  by  a doppler  radar  system.  The  evaluation  of  such  a sys- 
tem and  the  results  obtained  with  the  small-scale  prototype  DABS,  constructed 
by  the  University  of  New  Mexico  Civil  Engineering  Research  Facility  for  the 
Air  Force  Weapons  Laboratory,  are  presented.  A description  of  the  radar  unit, 
the  testing  of  Its  frequency  response,  a theoretical  look  at  the  possible  per- 
turbation of  the  signal  due  to  the  blast  wave  Itself,  and  the  experimental  de- 
sign and  results  of  the  prototype  DABS  radar  experiment  are  Included. 


RADAR  UNIT 

The  radar  unit  was  a Model  JF  100  Speed  Gun  manufactured  by  CMI,  Inc.,  Minturn, 
Colorado.  In  Its  normal  configuration  It  Is  a hand-held,  traffic-control,  dop- 
pler radar  device  which  gives  a digital  display  of  velocity  In  miles  per  hour. 
The  pertinent  specifications  as  given  by  the  manufacturer  are  as  follows: 

(1)  Frequency  — 10,525  ± 25  MHz 

(2)  Polarization  — Circular 

(3)  Beam  Width  --  8°  or  Less  (Side  lobes  suppressed  greater 

than  24  dB  down) 

(4)  Antenna  — 4-Inch  Aperture;  Circular  Horn  Type 

(5)  Microwave  Oscillator  — Solid-State  Gunn  Diode 

(6)  RF  Power  — 20  mW  Minimum;  100  mW  Maximum 

(7)  Receiver  Diode  --  Shottky  Barrier  Type  Rated  for  100-mW  Burnout 

The  unit  was  ordered  without  the  digital  display  and  Its  associated  electronics 
since  they  were  not  pertinent  to  the  experiment.  Thus,  the  total  unit  consist- 
ed of  an  antenna,  a preamplifier,  an  ampllfler/regulator,  a 9-V  regulator,  and 
the  housing.  The  system  was  powered  by  a 12-V  car  battery. 
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The  signal  was  taken  from  the  output  of  the  preamplifier  and  ac-coupled  (high 
pass  30  Hz,  3-dB  point)  because  of  a dc-offset.  The  only  other  modification 
made  was  the  removal  of  all  low-pass  filtering  from  the  preamplifier.  The  out- 
put Impedance  of  the  device  was  1000  Q.  This  recorded  signal  Is  the  amplified 
doppler  signal.  Its  relationship  to  the  velocity  of  the  target  Is 


(Al) 


where  f^  Is  the  doppler  frequency,  V Is  the  velocity  of  the  target,  f^  Is  the 
radar  frequency,  and  C Is  the  speed  of  light.  If  the  unit  of  velocity  Is  feet 
per  second  and  the  frequency  Is  In  Hertz,  eq.  (Al)  reduces  to  f^  = 21.40V. 


FREQUENCY  RESPONSE  TEST 

Since  the  system  is  normally  used  to  track  motor  vehicles  at  speeds  up  to  200 
mph,  it  was  deemed  appropriate  to  use  targets  with  velocities  approaching  those 
measured  in  a shock  wave.  The  fastest  target  readily  available  was  a rifle 
bullet  which  weighed  150  grains.  A Remington  Model  700,  7-mm  Magnum  was  used 
to  propel  the  target. 

The  rifle  was  aimed  so  that  the  bullet  would  pass  within  a few  inches  of  the 
radar  antenna,  and  the  doppler  signal  was  recorded  on  a Bell  & Howell  Model 
3700B  Tape  Recorder.  First,  the  bullet  was  fired  approximately  parallel  to 
the  radar  beam,  both  toward  and  away  from  the  radar  unit.  The  measured  velo- 
city was  in  both  cases  3198  ft/sec.  Since  the  listed  muzzle  velocity  of  the 
rifle  is  3260  ft/sec,  this  result  was  assumed  to  be  accurate.  However,  to 
further  check  the  system,  the  radar  unit  was  placed  so  that  the  angle  between 
the  beam  and  the  bullet  path  was  approximately  20  deg.  The  results  of  this 
test  are  shown  in  figure  Al;  the  measured  velocity  was  2803  ft/sec.  Since  the 
ratio  of  this  velocity  to  the  straight-on  velocity  equals  the  cosine  of  19  deg, 
it  was  verified  that  the  radar  unit  had  tracked  the  bullet. 

The  only  point  of  possible  concern  was  the  magnitude  of  the  voltage.  The  volt- 
age was  approximately  0.2  V peak-to-peak  (fig.  Al).  If  the  system  records  lower 
velocity  inputs  from  a moving  bullet  close  to  the  radar  unit,  the  preamplifier 
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0.2  msec/di V 

Figure  A1.  Doppler  Signal  from  Rifle  Bullet 


will  saturate  and  give  a signal  of  ± 5 V.  Thus,  the  signal  will  be  attenuated 
by  27  dB.  However,  It  was  felt  that  the  magnitude  was  more  than  sufficient  and 
signals  generated  from  targets  moving  at  velocities  of  4000  to  5000  ft/sec 
should  be  easily  discernible. 


BLAST  WAVE  PERTURBATION 

Since  the  radar  beam  tracks  a target  Inside  a shock  wave,  the  effects  of  the 
shock  wave  on  the  microwave  beam  must  be  considered.  Such  effects  can  possi- 
bly arise  from  the  Ionization  and  the  Increase  In  air  density  behind  the  shock 
front. 


The  effect  of  Ionization  on  the  Index  of  refraction,  n,  of  an  Ionized  medium 
can  be  expressed  to  the  first  approximation  by 


0)1 

n*  = 1 - 

0) 


(A2) 
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where  co  is  the  Incident  radiation  and  is  the  plasma  frequency  (ref.  1). 

This  relationship  shows  that  when  the  plasma  frequency  is  greater  than  the  in- 
cident frequency,  the  index  of  refraction  becomes  imaginary  and  the  medium  to- 
tally reflects  the  incident  radiation.  In  a cgs  system  of  units,  the  plasma 
frequency  may  be  expressed  by  the  relationship 


0) 


p 


4ffae^ 

m 


(A3) 


where  a is  the  number  of  free  electrons  per  cubic  centimeter,  e is  the  electron 
charge,  and  m is  the  mass  of  an  electron.  Thus,  in  this  radar  experiment,  in 
which  the  incident  radiation  was  10.525  OHz,  the  electron  density  would  have  to 
be  approximately  3.5  x 10*“  electrons/cm*  for  this  condition  of  total  reflec- 
tion to  occur.  Saha  (ref.  2)  and  Predvoditelev  (ref.  3)  indicate  that  electron 
densities  are  much  less  than  this  for  the  pressures  and  temperatures  expected 
inside  the  shock  wave  produced  by  a DABS  event.  However,  this  disagreement 
does  point  out  one  important  design  criterion  for  the  radar  system;  viz.,  it 
must  be  constructed  so  that  the  radar  target  intercepts  all  or  nearly  all  of 
the  microwave  beam  since  the  hot  gases  generated  by  the  detonation  used  to 
drive  the  shock  wave  will  also  reflect  the  beam  and  thus  generate  an  addition- 
al doppler  signal. 


The  second  perturbation  to  be  considered  is  the  effect  of  the  change  in  air 
density  in  the  shock  wave  on  the  doppler  signal.  The  origin  of  this  effect 
can  be  seen  if  one  rewrites  eq.  (Al)  as 

f j = - 2[d(nx)/dt]/X  (A4) 

where  x is  the  distance  between  the  radar  unit  and  the  target,  and  A is  the 
wavelength  of  the  radar  beam.  If  it  is  assumed  in  the  first  approximation 


1.  Jackson,  J.  D.,  Clasaiaal  Eleatvodynamiaa ^ John  Wiley  and  Sons,  Inc.,  New 
York,  N.Y.,  1962. 

2.  Saha,  M.  N.,  and  Srivastava,  B.  N.,  A Treatise  on  Heat,  3rd  Edition,  Indian 
Press,  Ltd.,  Calcutta,  India,  1950. 

3.  Predvoditelev,  A.  S.,  Phyaiaal  Gas  Dynamioa,  Pergamon  Press,  New  York,  N.Y., 
1961. 
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The  percentage  of  error  In  the  particle  velocity  measurement  because  of  this 
effect  Is 


n 


As  an  example.  If  the  Index  of  refraction  Increases  by  10  percent  and  the  ratio 
of  the  shock  velocity  to  the  particle  velocity  Is  about  5:4,  the  error  in  the 
particle  velocity  measurement  is  only  2.5  percent  with  n^  1. 


The  change  in  density  will  cause  a change  in  the  index  of  refraction  that  is 
much  smaller  than  that  estimated  in  the  example.  If,  as  a first  approximation, 
one  uses  the  Drude  model, 

^ _ dn 

p “ n - 1 

where  p is  the  density;  since  for  a ratio,  y*  of  the  isobaric  specific  heat  of 
air  to  the  isovolumetric  specific  heat  of  air  of  1.4  the  density  only  changes 
by  a maximum  factor  of  6 (ref.  4)  and  the  index  of  refraction  of  air  at  stand- 
ard temperature  and  pressure  is  only  slightly  larger  than  1,  the  change  in  the 
index  is  quite  small. 


The  last  point  to  consider  is  the  reflectivity  of  the  shock  front  itself  com- 
pared to  the  target.  But  since  the  reflectivity,  R,  of  the  front  may  be  shown 
to  be 


and  the  change  in  the  index  of  refraction  at  the  shock  front  is  quite  small, 
reflectivity  of  the  shock  front  is  negligible.  Thus  for  the  expected  condi- 
tions encountered  in  DABS,  the  blast-induced  perturbations  on  the  measurement 
of  the  particle  velocity  are  insignificant. 


4.  Zel'dovich,  Ya.  B.,  and  Raizer,  Yu.  P.,  Physios  of  Ohook  Wavcj  and  High- 
Tmperature  Hydrodyrumia  Phenomena^  Vol . I,  Academic  Press,  New  York, 
N.Y.,  1966. 
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PROTOTYPE  DABS  RADAR  EXPERIMENT 


Basically  a disposable  shock  tube  constructed  from  a corrugated  pipe  arch  was 
used  in  the  prototype  DABS  radar  experiment.  The  pipe  was  160  ft  long  and  its 
cross  section  was  approximately  oblate  {==  8 ft  wide  by  5 ft  high).  Figure  A3 
shows  the  basic  layout  for  the  radar  particle-velocity  experiment. 

The  target  consisted  of  three  layers  of  0.001-in-thick  aluminum  foil,  spaced 
0.5  in  apart.  The  target  was  placed  normal  to  the  direction  of  shock  propaga- 
tion and  completely  blocked  the  tube.  This  target  design  was  chosen  for  the 
following  reasons: 

(1)  Rise  Time.  The  round-trip  shock  transit  time  in  each 
0.001-in  foil  is  less  than  10  nsec.  Thus,  with  6 to  10 
transit  times  allowed  for  the  foil  to  come  up  to  the 
particle  velocity,  the  rise  time  of  the  target  becomes 
60  to  100  nsec.  Since  this  figure  is  much  better  than 
the  effect  of  possible  misalignment  and  lack  of  planar- 
ity with  respect  to  the  shock  front  and  since  the  ex- 
pected doppler  frequencies  are  in  the  range  of  90  to 
100  kHz,  the  target  should  be  adequate. 

(2)  Opacity.  Since  the  target  as  previously  explained  must 
shield  the  microwave  beam  from  the  detonation  products, 
a three-layer  spaced  system  was  chosen.  It  was  felt 
that  even  after  the  possible  turbulence  in  and  behind 
the  shock  front  had  shredded  the  target,  the  three-layer 
system  would  not  allow  the  microwave  beam  to  propagate 
through  it;  i.e.,  the  gaps  in  the  target  would,  for  the 
most  part,  be  smaller  than  the  half  wavelength  of  the 
microwave  beam  (~  1.4  cm). 

The  radar  unit  was  placed  20  ft  downstream  from  the  target  and  aimed  directly 
at  its  center  (fig.  A4).  Since  the  beam  width  was  only  8 deg,  the  reflected 
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Figure  A4.  Location  of  Radar  Unit  for  Test 


signal  came  from  the  target  and  not  from  the  walls  of  the  tube.  The  doppler 
signal  was  recorded  on  two  channels  of  a Bell  & Howell  Model  37006  Tape  Re- 
corder. One  of  these  channels  was  ac-coupled  (high  pass  10  kHz,  3-dB  point); 
this  channel  was  a Bell  & Howell  Wideband  II  FM  Channel  wHh  a 500-kHz  response. 
The  unfiltered  channel  was  a Bell  & Howell  Wideband  I FM  Channel  with  an  80-kHz 
response. 

The  results  of  the  experiment  are  shown  in  figure  A5.  The  initial  low-frequency 
signal  was  probably  caused  by  wall  vibrations  coupling  into  the  target  before 
the  arrival  of  the  shock  front.  The  period  of  the  data  was  measured  with  a 
traveling  microscope  and  the  data  were  divided  into  10  approximately  equal 
segments . 

The  average  period  measured  was  10.41  ysec  (96  kHz)  and  the  standard  deviation 
was  0.36  usee.  This  corresponds  to  a particle  velocity  of  4489  ft/sec. 


Figure  A5.  Doppler  Radar  Signal 


Figure  A6.  Particle  Velocity  Versus  Overpressure 


With  the  relationships 

/(y  - 1)  + (y  + 1)h^ 

“s  ’ 'oY S 

and  p 

(y  + 1)  + (y  - l)p^ 

Us  - Up  = u^ ^ (ref.  4) 

(y  - 1)  + (y  + 1)k^ 

0 

where  is  the  absolute  pressure  In  the  shock  wave  and  P^  Is  the  absolute 
pressure  of  the  unshocked  air  (12  psi  In  and  around  Albuquerque,  New  Mexico) 
and  with  a y of  1.4,  the  curve  of  oarticle  velocity  versus  overpressure  (fig. 
A6)  can  be  generated.  The  average  value  of  327  psi  agrees  quite  well  with  the 
Incident  pressure  measured  In  DABS.  At  the  closest  station  {-  20  ft  upstream), 
the  measured  pressures  varied  from  320  to  350  psi.  The  shock  velocity  approx- 
imately 10  ft  upstream  (5765  ft/sec)  indicates  a pressure  of  345  psi. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  prime  conclusion  from  this  study  Is  that  the  doppler  radar  system  can  meas- 
ure the  particle  velocity  In  the  airblast  wave  produced  by  DABS.  However,  It 
should  be  noted  that  the  scope  of  this  effort  only  permitted  the  measurement 
of  particle  velocity  at  one  point  (327  psi)  but  It  certainly  can  be  assumed 
that  the  system  will  function  at  lower  pressures  than  this.  The  system  should 
be  adequate  at  pressures  of  at  least  600  psi,  since  Ionization  should  not  pose 
a problem.  However,  this  conclusion  must  await  experimental  verification. 

A more  adequate  frequency  analysis  of  the  data  should  be  undertaken  In  order  to 
reduce  the  scatter  In  the  data  and  develop  a velocity  profile.  To  accomplish 
this,  high-frequency  digitization  (1  MHz)  of  the  tape  and  Fourier  analysis  of 
the  signal  are  presently  being  considered. 

Since  the  output  Impedance  of  the  device  Is  relatively  high.  It  should  be  buf- 
fered to  a much  lower  Impedance  so  that  the  line  may  be  terminated  to  obtain  a 
better  frequency  response. 

Experiments  In  which  metal  reflectors  are  used  for  the  microwave  beam  should  be 
undertaken  since  they  would  shield  the  radar  unit  from  the  flow  field  and  thus 
the  unit  could  be  reused. 
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INTRODUCTION 


When  a large  amount  of  explosive  is  detonated  at  one  end  of  a 
tunnel,  a shock  wave  with  velocity  up  to  Mach  6 or  7 is  formed  with 
corresponding  peak  pressure  in  the  atmosphere.  The  overpressure 
above  the  ambient  constitutes  the  shock  front,  along  with  rapid 
fluctuations  in  the  dynamic  pressure  occurring  behind  the  shock  front. 

The  phenomena  associated  with  such  a blast  consist  of  waves 
moving  much  faster  than  sound,  in  addition  to  "well  behaved"  acoustic 
waves,  large  but  finite  changes  in  pressure,  density  and  temperature, 
wave  fronts  causing  jumps  or  discontinuities  in  these  parameters  and 
finite  particle  velocities  associated  with  these  waves.  These  conditions 
are  not  consistent  with  acoustic  wave  theory  and  empirical  measure- 
ments and  calculations  are  required  to  determine  the  time  variation  of 
these  quantities. 

In  order  to  perform  such  calculations,  data  are  required  for  the 
variations  in  temperature  and  in  pressure  at  different  points  in  the 
tunnel.  What  is  wanted  is  the  time  history  of  the  pressure  and  its 
spatial  distribution,  a 4- vector  quantity.  The  response  time,  or  band- 
width, of  pressure  transducers,  and  their  susceptibility  to  destruction 
in  the  blast,  plus  other  limitations,  such  as  temperature  dependence 
of  the  pressure  response,  militate  against  their  use  in  the  main  stream 
of  the  tunnel  blast.  Interference  with  the  flow  behind  the  shock  front 
by  the  transducer  cannot  be  tolerated.  No  suitable  omnidirectional 
blast  transducers  exist  for  obtaining  free-field  time  histories  of  the 
pressure  in  the  blast  wave. 

Therefore,  an  approach  to  obtaining  the  desired  data  through 
velocity  measurement  and  temperature  measurement  has  been  proposed. 
The  measurement  of  velocity  time  histories  at  various  points  in  the  flow 
behind  the  shock  front  of  the  blast  is  the  subject  of  this  report.  It  is 
proposed  to  perform  such  measurements  using  a microwave  analog  of 
a laser  Doppler  velocimeter  (LDV). 
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MICROWAVE  DOpPLER  VELOCIMETER  THEORY 


A.  Background 

Two  types  of  velocity  measurement  approaches  using  Doppler 
velocimeters  are  covered  in  Refs.  3 and  4.  The  first  (Ref.  3)  optically 
homodynes  the  scattered  energy  from  a moving  particle  in  a fluid  with  a 
reference  beam  unscattered  by  the  medium.  The  Doppler  difference 
frequency  is  proportional  to  the  velocity  of  the  medium  at  the  scattering 
point,  which  is  known  from  the  scattering  angle. 

The  second  type  is  described  in  Ref.  4.  The  separate  coherent 
beams  are  focused  on  the  same  small  volume  in  the  fluid  from  different 
angles,  and  two  different  scattered  radiations  are  heterodyned  to  detect 
a Doppler  difference  frequency  which  is  proportional  to  the  velocity  in 
the  medium  for  the  volume  illuminated. 

B.  Theory 

' The  theory,  as  developed  in  Ref.  4,  is  readily  adapted  to  a microwave 
analog  of  the  laser  system  employing  two  beams.  This  theory  will  be 
expanded  to  provide  insight  to  the  design  of  a suitable  velocity  measure- 
ment system  for  application  to  the  tunnel  blast  wave  measurement 
problem.  Figure B1  illustrates  the  geometry  of  the  problem.  The  antenna 
apertures  are  assumed  to  be  circular  with  Gaussian  illumination  so  as 
to  produce  a Gaussian  pencil  beam  pattern  of  diameter  2bo  at  the  1/e^ 
level.  FigureB2a  illustrates  the  diffraction  region  for  each  beam  and  the 
intersecting  volume  V of  Fig.Bl.  The  two  radiations  are  mutually  coherent, 
being  derived  from  the  same  source  over  identical  path  lengths,  and  with 
the  same  linear  polarization.  The  planar  phase  fronts  from  the  two 
beams  interfere  over  the  volume  Ax  Ay  Az(  1 /e^) . A set  of  closely- 
spaced  planar  interference  fringes  in  the  beam  crossover  region  occur 
as  indicated  in  Fig.  B2b.  The  spacing  of  these  interference  fringes  is 
given  by  the  half- wavelength  A/ 2,  divided  by  the  sine  of  the  half  angle 
between  the  two  beams.  Figure B2c indicates  the  relative  signal  amplitude 
regions  for  the  1/e  and  1/e^  levels.  The  energy  from  either  beam  within 
this  ellipsoid  will  be  at  or  above  86.  5%  of  maximum.  At  the  parallel 
lines  of  phase  reinforcement,  the  two  waves  will  add  constructively, 
elsewhere  they  will  add  in  various  degrees  of  phase  opposition. 

Separation  between  the  maxima  is  y = ' ^ — , where  9 is  the  angle 

2sin(9/2) 

between  the  beams. 
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Figure  B1 . Conceptual  Microwave  Velocineter 
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Figure  B2.  Laser  Beam  Crossover  Region  (Ref.  4) 


1/e^  Relative  Signal 
Amplitude  Ellipsoid- 


Figure  B2  (concluded)  (Ref.  4) 


A particle  paaaing  through  this  region  (see  Fij.B3)  with  velocity 
components  vy  only,  will  cut  the  fringes  at  a rate  proportional  to  its 

velocity  and  thus  produce  an  alternating  scattered  signal  at  a frequency 
proportional  to  the  velocity.  The  only  component  of  velocity  detected  is 
that  oriented  normal  to  the  interference  fringes. 

Since  the  amount  of  energy  scattered  in  any  direction  is  solely  a 
function  of  the  fringe  pattern  intensity  at  the  point  of  scattering,  the 
frequency  is  independent  of  the  direction  of  scattering.  This  means 
that  a larger  scattering  volume  on  receive  will  increase  the  sensitivity 
of  the  detection  system. 


The  number  of  cycles  (N)  of  information  obtained  from  a scatterer 
passing  through  the  center  of  the  probe  volume  is  given  by 


^(l/e^)  = where  N is  the  number  of  cycles  between  the 

A/ 2 

l/e^  relative  signal  amplitude  points,  2bQ  is  the  1/e^  beamwidth,  and 


# is  the  angle  between  beams.  Obviously  more  cycles  (and  greater 
accuracy)  result  when  A is  small,  9 and  bg  are  large.  On  the  other 


hand,  bQ  large  reduces  the  positional  resolution  of  a measurement.  It 

is  reemphasized  that  this  development  applies  only  to  the  discrete  case 
of  a single  scatterer. 


C.  Differential  Doppler  Shifted  Scattering 

The  Doppler  shift  in  the  energy  scattered  into  any  particular  direction 
from  beam  1 (Fig.B4)i  is  given  by 


Af,  = y (f.  - *i)  • V 


(B1) 


where  f^  is  the  Doppler  shift  in  energy  scattered  in  the  direction  denoted 

by  the  unit  vector  a^,  Fj  is  a unit  vector  in  the  direction  of  beam  1,  and 

Y is  the  velocity  vector  of  the  scattering  center.  The  Doppler  shift  of 
energy  scattered  from  beam  2 in  the  direction  of  a^  is 

Af2  = (a^  -'i2)  • V (B2) 
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Figure  B3.  Cross  Section  of  Elllpsoldil  Interference  Region 
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S 


Figure  B4.  Doppler  Difference  Me-isurement  Geometry 
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The  output  of  a mixer  with  inputs  of  scattered  signals  from  B1 
and  B2  includes  a difference  term  between  these  Doppler  shifts, 
oscillating  at  a frequency  f(j: 

*cl  ""  ^^2  “ (ai  - ^2)  • V 


The  Doppler  difference  frequency  is  independent  of  ag,  and  as 


previously  noted,  the  detector  antenna  position  can  be  chosen  for 
maximum  collection  of  scattered  energy.  This  frequency  is  also 
proportional  to  the  velocity  component  in  the  direction  of  a~^  - 
which  is  perpendicular  to  the  bisector  of  the  two  incident  beams  and 
in  the  plane  of  the  two.  For  the  situation  of  Fig.Bl,  the  sensed  velocity 


component  Vg  is  given  by 


^1 


2 sin 


9 

T • 


The  differential  Doppler  signal  will  be  modulated  in  both  phase 
and  amplitude,  and  will  require  analysis  of  the  recorded  output  of  the 
homodyne  detector. 

The  frequency  spectrum  broadening  over  any  given  observation 
time  is  due  to  the  residence  time  of  a scattering  center  within  the 
defined  scattering  volume  or  transit  time,  spatial  variations  in  the 
velocity  within  the  sample  volume, and  phase  and  amplitude  variation 
differences  for  the  two  transmit  beams  over  their  separate  propaga- 
tion  paths. 

The  mixer  output  power  is  proportional  to  the  square  of  the  incident 
fields  scattered  from  the  two  beams; 

P(t)- [Ei(t)  + E2{t)j2  (B4) 

The  power  spectrum  P(Kpw)  of  the  mixer  output  is  the  Fourier  trans- 
form of  the  autocorrelation  function  for  a stationary  random  process; 

P(Kj,w)  = Rg/  <P(o)  P(r)>  e-Jw^’dr  (B5) 

— « 


where  the  brackets  denote  the  average  of  the  ensemble  for  the  auto- 
correlation function.  Then 


<P(o)  P(r)>-<[Ei(o)  + E2(o)]^[Ej(t)  + E2(t)]^> 


(B6) 
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The  transmitted  energy  at  the  scattering  volume  has  a time  varying 
waveform 


Ej(t)  = E2(t)  = A(t)  costw^t  + o(t)]  VB/; 

since  the  wave  is  affected  by  the  time  varying  character  of  the  propaga- 
tion medium.  Since  the  scattered  fields  at  the  receiving  antenna  are 
detected  as  the  sum  of  Ej(Q  and  E2(t).  it  will  have  time  varying 

amplitude  and  phase  depending  on  the  position  and  motion  of  the  scattering 
centers,  and  will  be  of  the  form: 


E(t)  = B(t)  cos  |[wq  + Aw(t)]t  4-  Aa(t)  + 


(B8) 


where  a(t)  is  the  varying  phase  angle  due  to  propagation  path  phase 
shift  of  the  two  transmit  beams  over  their  separate  propagation  paths, 
Aw(t)  is  the  varying  Doppler  shift  frequency  component  due  to  the 
composite  motion  effects  of  the  scattering  centers,  and  Aat(t)  is  a 
random  phase  shift  caused  by  the  random  positions  of  the  scattering 
centers. 

If  B6  is  expanded  we  obtain 

<P(o)P(r)>-<[Ei2(o)  + 2Ej(c)E2^(o)]  [Ei2(r)+2Ei(f)E2(t)+E2^(f)]>  (B9) 


Without  performing  the  expansion,  it  can  be  noted  that  certain 
terms  will  contribute  to  spectrum  spreading  and  others  will  not.  Of 
the  nine  terms,  four  containing  factors  of  the  following  form  integrate 
to  zero  value: 

<cos[Aa(o)]  > , < cos  [Aa(r)]  > , <cos  [Aa(o)  4- Aa(r)]> 

This  is  because  the  phase  shift  takes  on  all  values  and  the  average 
value  of  the  sinusoids  is  zero. 

The  terms  < E|^(o)Ej^(t)  > and  <E2^(o)E2^(r)  > will  result  in 

spectral  components  at  frequency  2(m^  -f  A<tf)  xnd  zero  (dc).  The  upper 

one  will  not  appear  in  the  mixer  output.  The  dc  terms  have  a varying 

amplitude  of  the  form  <'(l/4)[B2(o)B2(r)J^which  is  due  to  the  amplitude 
variation  in  the  transmit  beams  and  to  the  varying  positions  of  the 
scattering  centers.  It  is  not  due  to  the  Doppler  shift  effect,  or  velocity 
of  the  scattering  centers. 


98 


Two  other  nonzero  terms  at  zero  frequency  (dc)  result  from  the 
crosscorrelation  of  the  two  amplitudes  and  B2> 

The  final  term  contains  the  information  of  interest 

^E^(o)£2(r)  E2(o)E^(r)^ 

= ^ 2MB^(o) B2(r)B2(o)B|^(r)cos  I [^W2(r)  - Aw^(7)]r  (BIO) 

+ Aa2{r)  - Ao^(r)  + « 2(’')  " 


where 


M 


cos  [Aoj  (o)  + aj(o)]co8[Aot2(o)  + 


(BID 


Now  rewriting  BIO 
<EDo)E2(r)E2(o)Ej(r)> 

= < 4Bi(o)B2(r)B2(o)Bj(r)  |cos[Aaj{o)  + 0^(0)  + Ao^(o)  + O2(o)] 

+ co8[  Aafj(o)  + aj(o)- A«2(o)  - 02(0)]  I cos  | [Aw2(r)  -Awj(r)]r 

+ Ao2(t)  - Aaj(T)  + a^{r)  - aj(*’)|> 

and  it  becomes  apparent  that  there  are  numerous  cross-products  involving 
A a and  a . If  the  difference  phase  shift  represented  for  example  by 
the  argument 

[ Aaj(o)  + aj(o)  + Aa^{o)  + 02(0)  -Aa2(»^)  + Aa^  (r) 

- 02(t)  + «i(^)] 

is  uncorrelated  for  large  values  of  r,  then  the  variable  in  the  brackets 
is  a random  one,  and  its  cosine  averages  to  zero.  This  effectively 
would  destroy  the  term  in  [Aw2(^)  - AwD''')^  transit  time  of 

particles  across  the  sampling  volume  is  ty,  then  at  frequencies  less 

than  1/ty  the  random  phase  effects  will  mask  the  Doppler  shift 

frequency  difference  which  it  is  desired  to  detect.  Therefore,  the 
sampling  volume  should  be  kept  as  large  as  possible. 
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At  a velocity  of  300  fpa  or  91.4  m/aec,  a 10  cm  observation  zone 
results  in  a transit  time  ty  = 1.09  msec,  which  approximately  corresponds 

to  a frequency  of  1000  Hz.  It  is  unlikely  that  the  Doppler  difference 
frequency  would  be  observable  below  this  velocity. 

It  must  be  realized  that  Fourier  analysis  of  the  detector  output 
will  be  required  to  extract  the  velocity  information  from  the  broadened 
spectrum  due  to  the  random  phase  shifts  with  time.  The  receiver 
amplifiers  should  be  designed  to  accommodate  a spectrum  broadening 
of  5 percent  greater  than  the  maximum  Doppler  difference  frequency 
requirement. 

In  addition  to  the  variable  phase  shift  with  time,  the  amplitude 
varies  with  time  and  also  with  position  because  of  the  shape  of  the 
scattering  volume.  This  further  contributes  to  the  spectrum  broadening. 
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PLASMA  EFFECTS 


A critical  consideration  is  the  cut-off  frequency  for  propagation 
through  a medium  likely  to  be  highly  dispersive,  because  of  the  high 
temperatures  in  the  combustion  products,  and  the  resulting  ionization, 
as  well  as  the  solid  constituents  in  the  medium,  such  as  dust  and  unburned 
cordite  particles. 

A.  Shock  Wave  Phenomena 

1.  Simple  Shocks  - Relatively  high  gas  temperatures  and  pressures 
may  be  achieved  (simultaneously)  with  "simple  gas  shocks",  or  in  "blast" 
waves".  In  a shock,  the  motions  are  so  rapid  that  very  little  momentum 
or  energy  is  transported  across  the  flow  lines.  A shock  front  is  rather 
thin.  In  the  shock  the  speed  of  sound  is  greater  at  the  density  maximum 
than  at  the  plateaus  on  either  aide  of  the  shock  wave,  consequently  the 
shock  wave  gets  steeper  and  steeper  as  the  maximum  catches  up  with  its 
own  base.  This  rapid  increase  in  pressure  poses  a serious  problem  in 
designing  radiating  structures  (antennas  or  radomes)  for  a microwave 
velocity  measurement  system. 

From  considerations  of  conservation  of  mass,  conservation  of 
momentum,  and  conservation  of  energy,  the  physical  relations  connecting 
the  media  across  the  shock  can  be  deduced.  The  resulting  equations  are 
known  as  the  HUGONIOT  relations. 

2.  Shock  Waves  - The  motion  of  a disturbance  relative  to  a fluid 
is  known  as  "wave  motion",  and  the  speed  of  propagation  is  the  "wave 
speed".  Consider  one-dimensional  motion,  that  is,  the  case  of  plane 
waves.  A shock  wave  propagating  into  a stationary  fluid  sets  the  fluid 
into  motion  and  raises  its  pressure,  temperature  and  density.  The  "ideal" 
gas  "jump"  (shock  wave)  conditions  can  be  related  by: 

p Ui  0+1) 

— =■  = — = - Density- Velocity 

^1  U2  ( > - 1)  Mj2+2 

zy  -7  + 1 

— = — ■ ■ ■ ■■  Pressure 

Pi  7+  I 

T^  [27  Mi^  - 7 + 1][(7-  l)Mj^+2j 

— = — ■ - - - Temperature 

Ti  ( 7 + 1)^  Mj2 
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where  1.  2 denote  physical  quantities  in  front  of  and  behind  the  shock, 

Ui 

respectively;  Mi  = Mach  number  in  region  1,  equal  to  * and 

Csi 

c,^  = . y is  the  ratio  of  specific  heats:  Y = , the  ratio  of 

the  specific  heat  at  constant  pressure  to  the  specific  heat  at  constant  volume. 

f+2 

y is  also  equal  to  ly  " . where  f is  the  number  of  degrees  of  freedom  for  the 
fluid.  The  shock  wave  speed  is  denoted  by  c,. 

A "weak  shock"  is  defined  as  one  for  which  the  normalized 
pressure  jump  is  very  small,  that  is: 


AP 

Pi 


Pz-Pi 

Pi 


<<  1.0 


A "very  strong  shock"  is  defined  as  one  for  which  the  pressure 
ratio  P2/P1  !■  very  large. 


3.  Strong  Shock  (Rankine-Hugoniot)  Equations 
For  very  strong  shocks.  P2/P1  >>1.  then 


^2 

Pi  y-  1 
I2  ^ 7-1 

Ti  r+ 1 


El 

Pi 


and  Up 


If) 

(■sfr;  If) 


1/2 


where  Pj  and  P2  are  the  mass 


densities  in  front  of  and  behind  the  shock, respectively;  T2  *nd  Tj  are 
the  temperatures;  ai  is  the  speed  of  sound  in  front  of  the  shock,  and  c^ 
is  the  shock  wave  speed.  Thus  the  Mach  number  of  the  shock  is 
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M 


I = 


The  fluid  velocity  behind  the  ihock  is  Up.  Based  on  data  obtained  to 
date  from  DABS  experiments,  the  following  values  will  be  used: 


Cg  = 6000  ft/ sec 


Pi 


P2  600  psi 
T2  ^2500°K 


Consider  a gas  with  3 degrees  of  freedom,  i.e. , f = 3,  then 

y = i-i-2.  s — . From  this 
f 3 

fz  , 1±±  = 5£3ii  , 4, 

Pi  y.  I 5/3  .1 

Pz 

Similarly  for  a gas  with  4 degrees  of  freedom,  ■ = 5. 

Consider  now  the  shock  velocity  (Cg  = 6000  ft/sec),  the  sound 
speed  in  front  of  the  shock  ai,  and  the  associated  mach  number 


M 


I 


For  c,  = 6000  ft/sec,  ai  = 1135  ft/sec  (standard  temperature  and  pressure  in 
air). 


Ml  = 


6000 

1135 


5.3 
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I 


t. 


Next  consider  P2  = 600  psi. 

Pj  = 14.5 

psi,  then 

T^ 

>-  1 ^ 

600  \ 

Tl 

Y+  1 Pi 

>+l  \ 

14.  5/ 

with 

> = 3/2.  -iS- 

- 3/2-1 

1 

11 

: 8.28 

Tl 

3/2+1 

V14.5/ 

If  = 300°Ki  then  T2  % 2500°K.  Thua  we  have  an  eatlmate  of  the 

teniperature  from  which  to  compute  the  electron  density.  For  Y = 1.4 
(air)  T2  ^ ZlOO^Ki  and  for  conservatism  we  will  use  ZSOO^K  in  subsequent 

calculations. 


B.  Electron  Density  and  Plasma  Frequency 

The  plasma  frequency,  due  to  free  electron  density  in  an 

pe 

ionised  medium  of  electromagnetic  wave  propagation,  represents  a 
critical  frequency  below  which  incident  E-M  waves  are  reflected  from 
the  plasma.  This  frequency  depends  on  the  electron  density, 

The  relationship  is 

— . s 56.4  n^^^^  radians/sec  (MKS) 

< 

o 


The  plasma  density  results  from  the  ionization  of  gas  atoms  or 
molecules  as  a result  of  the  thermal  condition  (temperature)  of  the  gas. 
In  deriving  an  expression  for  the  total  number  of  ionized  particles,  the 
recombination  of  the  positive  ions  with  the  electrons  must  be  considered. 
When  thermodynamic  equilibrium  exists,  the  number  of  new  ions  equals 
the  number  of  recombinations.  Using  this  condition,  SAHA  derived  an 
equation  giving  the  degree  of  ionization  a in  terms  of  the  gas  pressure 
P and  the  temperature  T. 


2.4  X 10"^ 
P 


T^-^  exp(-Wi/kT)  where  P is  in  Torr 
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(mmHg),  ii  the  ionisation  energy  in  electron-volts,  k is  Doltsmann's 

constant  in  (ev)/°K.  Having  determined  a,  n^  may  be  determined.  The 

equation  is  extremely  sensitive  to  the  value  selected  for  W^,  which  depends 
on  the  medium. 

Let 

Wj  = 15  ev,  k = X 10* ^ g ^ 10-,5ev 

1. 602  X 10“*9 


then  for  T2  = ZSOQOK.  exp(-Wj/kT)  = 5.6  X lO'^l. 

If  Wj  = 12  ev,  exp(-Wj/kT)  = 6.  3 X 10-25;  or  if  Wi  = 13  ev,  exp(-Wi/kT)  = 
6. 1 X 10-27. 


It  is  apparent  that  the  value  of  a is  quite  sensitive  to  the 
ionisation  potential,  so  an  accurate  value  of  W|^  for  the  composite 

medium  of  air,  and  combustion  products  mixed  with  it,  is  needed. 

If  a value  for  of  11  ev  is  chosen,  and  a2  1,  the  equation 
simplifies  to 


^2  _ (2.4  X 10^)  (2500)^’^ 
3.  I X 10^ 


(6.55  X 10-23)  ^ i 5g  X 


The  density,  with  standard  temperature  and  pressure,  is  2.685  X 10^^ 


Pz  2.4 


- 6.  Thus  the 


.4 


. <1  *• 
molecules/m-’,  but  at  600  psi, » I and 

Pi  <*i 

density  is  about  6 times  as  great.  The  ionized  electron  density  is  therefore 
n^  = (2.685  X lo25)(6)(l.  25  X 10-7)  = 2 X I0^‘^electrons/m3 
For  this  value  of  n^,  obtained  as 

(4>p^  = 56.4  ne^^^rad/sec  = 8.976  n^^^^  Hs 


f pg  = 40.  3 GHz 
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For  other  values  of  we  obtain 


Wi(ev) 

o2 

10 

1.65  X 10-^2 

11 

1. 58  X 10-1^ 

12 

1. 52  X 10-^6 

13 

1.47  X 10-»8 

14 

1.41  X 10-20 

15 

1. 35  X 10-22 

ne(/m^) 

fpe(GHz) 

2.07  X 1020 

129 

2.02  X 10^9 

40.4 

1.99  X 10^8 

12.7 

1. 74  X 10l7 

3.74 

1.91  X 10‘^ 

1.24 

1.87  X 10^5 

0.38 

If  the  ionization  potentiai  and  the  assumed  temperature  and  pressure 
produce  an  electron  density  of  10^ ^/m^,  then  any  frequency 
in  X-band  or  above  should  propagate  satisfactorily  through  the  medium. 

1 . Plasma  Conductivity  (skin  depth)  - The  relationship  between 
the  conductivity  and  the  temperature  (in  the  absence  of  an  external 
magnetic  field)  is  given  by  Spitzer's  equation: 

t3/2 

9 = O.OIS  mho/m 

Z InA 

with  T the  electron  temperature  in  °K,  Z the  effective  ionic  charge, 
and  A is  given  by 

12  W («okT/e2)^'^ 

where  is  free  space  permittivity,  e is  electronic  charge,  n^  is  number 

density  of  electrons.  Thus  9 depends  on  temperature  and  electron  density. 
For  the  numbers  used  previously,  i.  e.  T = 2500°K,  and  selecting  n^  as 

lO^^/cm^,  we  obtain; 

T 3/2 

A = 4.9  X 10^^  ' — , with  T.  in  kev  and  n®/m^ 

n.>/2 


106 


T = 2500OK  = _i522-  ev  = .215  X 10"^  kev 
11.600 

A = (4.9X  10^^)  r 1544 

(10l2)l/2 

lnA>  7.342 

In  thes«  units  (kev  and  cm), 


9 


(4.8  X 10‘‘^)lnA 


(.215  X 10-^^^^ 
(4.8  X 10‘'^)(7.34) 


. 895  mho/cm 


= 89«  5 mho/m 

Recalling  conductivity  values  for  good  conductors: 


Medium 

mhos/m 

Sea  water 

3.5 

This  plasma 

89.5 

Aluminum 

3.54  X 10 

Copper 

5.8  X 10*^ 

It  is  apparent  that  the  plasma  is  a fairly  poor  conductor,  but  still 
dispersive. 

2,  Plasma  Skin  Depth  - The  skin  depth  6 is  calculated  from  the 
conductivlty>-assumed  to  bo  90  mhos/m,  and  is  based  on  the  following 
equation. 

4 = (irf 

where  f is  the  frequency  of  the  E-M  wave,  \i  is  the  permeability  (assumed 
same  as  air)  and  9 the  conductivity.  For  a frequency  of  10.  525  GHa 

« = [(ir)(l0.525  X l09)(4irX  I0‘7)(89)r^'^ 

s 5.  2 X 10"^  m or  approximately  1/2  mm. 

This  indicates  that  the  E-M  wave  is  absorbed/ reflected  by  the  relatively 
good  conductors  in  the  hottest  part  of  die  plasma. 
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EM  WAVE  SCATTERING 


The  various  types  of  scattering  to  be  considered  (since  scattered 
energy  must  be  detected  to  obtain  the  velocity  shifted  Doppler  frequency) 
are: 


Scattering  by  electrons 

Scattering  by  conducting  and  by  dielectric  spheres 
Scattering  by  "particles". 

A.  Scattering  by  Electrons  (Thomson  scattering;  Thomson  cross-section) 

If  a plane  wave  of  monochromatic  E-M  radiation  is  incident  on  a 
free  particle  of  charge  e and  mass  m^  (e.g. . an  electron),  the  particle 

will  be  accelerated  and  emit  radiation.  The  emitted  radiation  will  be  in 
directions  other  than  that  of  the  incident  plane  E-M  wave.  This  process 
may  be  described  as  "scattering"  of  the  incident  E-M  wave.  The  Thomson 
cross-section,  is  given  by^ 

d = (sin<^)^ 

4irm 

in  MKS  units  for  a single  electron, with  l^the  angle  between  the  direction 
of  polarisation  of  the  incident  field  and  the  direction  of  the  scatterer 
to  the  receiver.  This  reduces  to 

or  = 8 X lO’^O  n 

«T  e 

and  for  the  case  in  question,  assuming  a cubic  scattering  volume  of  10  cm 
on  a side, we  obtain  a maximum  value  ( ^ = 90°)  at  10.4  GHz. 

^ ^ 10*^0)(1.35  X 10l8)(10-3)  = 1.08  X lO’^^m^ 

= -110  dBm^,  an  undetectable  cross-section. 

B.  Scattering  by  Conducting  and  by  Dielectric  Spheres 

1.  Dielectric  Spheres  - The  medium  within  the  tunnel  after  the 
shock  front  passes  is  known  to  be  composed  of  the  (partially  ionized) 
ambient  air,  dust  from  the  tunnel  floor  and  roof,and  combustion 
products.  We  will  examine  the  backscattering  potential  of  the  dust 
particles  and  combustion  particles  in  terms  of  scattering  of  an  E-M 
wave.  The  instrumentation  situation  under  discussion  is  a bistatic  one. 
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with  biiUtic  scattering  angle  9/2.  A theorem  due  to  Crispin,  et  al, 
states  "In  the  limit  of  vanishing  wavelength  the  bistatic  cross-section 
for  the  transmitter  direction  and  receiver  direction  n^  is  equal  to 

the  monostatic  cross-section  for  the  transmitter- receiver  direction 
-f  ^ withl^^  ^ for  bodies  which  are  sufficiently  smooth".  'Ic  is  a unit 
vector  from  transmitter  to  target,  is  a unit  vector  from  receiver  to  tarpet. 

Thus,  in  our  case,  the  monostatic  cross-section  should  be  close  to  the 
actual  bistatic  cross-section  for  individual  particles  (eapecially  if  the 
particles  are  spherical),  since  the  bistatic  angles  will  be  less  than  90°, 
and  the  wavelengths  of  interest  are  of  centimeter  dimensions. 

In  the  Rayleigh  region,  a sphere  of  circumference  small  by 
comparison  to  the  wavelength  exhibits  a cross-section  varying  as 

1/A^  where  A is  the  wavelength.  Dust  particles  will  be  considered 
as  dielectric  spheres  of  various  diameters  uniformly  distributed 
throughout  the  scattering  volume  but  with  random  phase.  The  probability 
of  the  scattered  signal  amplitude  x lying  between  x and  x dx  from  the 
radar  resolution  cell  is^ 

p(x)  dx  = ip-  exp  j 

where  is  the  time  average  of  x^.  This  is  a Rayleigh  distribution, 
with  shape  independent  of  the  component  amplitudes. 

One  of  the  properties  of  this  distribution  is  that  the  rms  amplitude 
of  the  fluctuations  about  the  mean  is  proportional  to  the  mean,  i.  e. , the 
ratio  of  the  standard  deviation  to  the  mean  value  is  a constant.  The 
above  assumes  that  the  incident  wave  is  of  constant  power. 

A potential  difficulty  in  obtaining  adequate  performance  from  the 
systemi  lies  in  the  unwanted  fluctuations  in  four-dimensional  probability 
space.  These  are  the  marginal  density  functions  of  the  two  transmitted 
beams  and  that  of  the  receive  beam  because  of  fluctuations  in  the  medium 
within  the  experimental  space,  and  the  marginal  density  of  the  volumetric 
radar  cross-section,  discussed  above.  If  a Rayleigh  distribution  describes 
each  of  these  phenomena,  which  are  not  independent,  then  a four-dimensional 
probability  space  must  be  examined  in  detail  to  determine  the  constitutive 
parameters  of  the  medium  required  to  minimize  unwanted  signal  fluctuations. 
These  could  effect  .vely  mask  the  desired  fluctuations  resulting  from  a 
single  particle  passing  through  clearly  defined  interference  fringe  planes. 


109 


It  is  apparent  that  the  scale  of  the  distribution  of  cross-section  in 
the  measurement  space  must  exceed  the  spacing  between  the  fringes. 

If  scattering  centers  are  uniformly  distributed  across  the  measurement 
volume  there  will  be  no  discernible  signal  variation  proportional  to 
velocity.  Similarly,  if  the  phase  velocity  of  the  transmitted  beams 
varies,  the  fringe  positions  will  vary,  and  erroneous  outputs  will 
result.  Finally,  if  the  receiver  amplitude  fluctuates  because  of 
fluctuating  attenuation  over  the  propagatioi  path, an  erroneous  output 
will  result.  Essentially  the  standard  deviation  to  mean  ratio  must  be  as 
small  for  all  four  signal  modifiers. 


2.  Cross-Section  of  Dielectric  Spheres  - In  the  Rayleigh  region 
the  radar  cross-section  of  a small  dielectric  sphere  of  diameter  D is: 


X4 


2 


where  K depends  on  the  dielectric  constant. 


K s (m2  . l)/(m2  + 1),  and  m is  the  index  of  refraction.  If  we 
assume  sandy  dry  soil,  the  dielectric  constant  is  2.  55  at  X-band 
frequencies,  with  a very  small  dissipation  factor  which  will  be  ignored. 
Then 

|k|  -4m  ‘ -437  and,for  fs  10.525  GH» 

_ (.437)2y5p6  ^ (g  10^)D^ 

(.0285)4 


In  DABS  experiment  1-D,  a debris  collector  collected  36.78  g of 
material.  Table  b1  is  reproduced  from  Ref.  6. 

TABLE  B1 

DEBRIS  TOTAL  WEIGHT  36.78  (g) 


Diameter  (in.) 

D(m) 

Mass(E) 

Percent  of  Total 

Greater  than  .0331 

84074  X 10>8 

9.2 

25 

Greater  than  .0165 

41910  X 10-8 

9.3 

25.3 

Greater  than  .0098 

24892  X 10-8 

7.0 

19 

Greater  than  .0058 

14732  X 10-8 

5.8 

15.8 

Greater  than  .0029 

7366  X 10-8 

2.0 

5.4 

Less  than  .0029 

3.5 

9-5 

36.8 

100 
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Since  about  half  the  particles  are  greater  than  .0165  in.  in 
diameter  and  the  croia  section  varies  as  we  will  calculate  the 
effective  cross-section  on  the  basis  of  the  area  density  for  these 
particles.  The  pipe  collector  diameter  was  1. 5 in. , so  for  the  largest 
particles  the  area  density  was  8070  g per  m^.  For  the  next  largest 
it  was  8157  g/m^.  Multiplying  by  the  areas  of  the  respective  particle, 
sizes  we  get  4.478  X 10'^  g/particle  and  1. 125  X 10~3  g/particle 
respectively.  This  leads  to  an  estimate  of  2.05  X 10^  and  8.  27  X 10^ 
particles  of  the  largest  and  next  largest  size.  These  particles  were 
deposited  over  an  interval  of  some  20  msec  at  an  average  velocity  of 
5000  fps,  or  an  effective  dimension  of  100  ft,  or  30.  5 m.  This, 
considering  the  area  of  the  pipe,  is  a volume  of  .035  m^.  The  volume 
densities  are  then,  respectively,  5.9  X 10^/m^  and  23. 6 X 10^/m^. 

In  Ref.  6,  the  point  is  made  that  some  of  the  particles  are  due  to 
fallback.  The  authors  assumed  that  the  largest  particles  are  due  to 
fallback.  For  the  sake  of  conservatism,  this  same  assumption  will  be 
made  for  this  analysis,  and  a particle  density  of  236  X lO^/m^  will  be 
used,  with  a particle  D of  41910  X 10’^  m. 

dj  = (8.86  X 10‘^)(4.191  X 10"^)^  = 4.  8 X 10'^  W 

Assuming  a scattering  volume  defined  by  a 1.5°  circular  beam 
and  1 psec  pulsewidth  at  a range  of  2 km,  we  get  7r(26.  2)2(150)  = 

323477  m^.  These  were  the  parameters  of  the  measurement  performed 
in  Ref.  7,  which  obtained  a measured  cross-section  of  .8  m^  on  a dust 
cloud.  This  results  in  a dy  volume  reflectivity  of  247  X lO'B 

Using  the  particle  density  for  the  DABS  experiment  data  we  obtain 
a volume  retlecllvity  ffy  o£  x io3)(4.8  X l0-‘3)m2/m’,  or 

*V  = 11- X 10-8m2/m3. 

This  is  within  one  order  of  magnitude  of  the  X-band  measurement  of 
Ref.  7. 


3.  Cross-Section  of  Conducting  Spheres  - If  a system  cannot  be 
designed  with  adequate  sensitivity  for  detecting  scattering  from  dust 
particles,  it  may  be  possible  to  increase  the  scattering  cross-section 
by  "seeding"  the  medium  with  aluminum  dust,  constituting  conducting 
spheres.  The  amount  and  size  of  the  particles  required  will  be  deter- 
mined at  that  time,  in  a similar  manner  as  in  1 and  2 above. 


C.  Scattering  by  Particles 


Since  the  particulate  matter  collected  in  the  DABS  experiment  of 
Ref.  6 was  ail  dusti  it  is  presumed  that  high  velocity  products  of  the 
cordite  explosive  charge  in  particulate  form  are  carried  out  of  the 
tunnel  with  the  blast  front  or  behind  it.  These  particles  may  increase 
the  volume  reflectivity  just  behind  the  blast  front.  Calculations  of  the 
radar  cross-section  due  to  this  could  be  made  from  radar  data  available 
to  DABS  testing,  provided  the  Doppler  radar  was  first  calibrated  using 
a precision  sphere.  This  has  not  been  done  to  date.  For  the  following 
analysis  a value  of  dy  “ 10“^  m^/m^  will  be  used. 
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BISTATIC  RADAR  EQUATION 


Pj.  - received  signal  power  (W) 

- transmitted  power  (W) 

Gt  - transmitting  antenna  gain  in  target  volume  direction 
Gf  - receiving  antenna  gain  in  target  volume  direction 
V - scattering  volume  defined  by  antenna  beams  (m^) 

- transmitter  to  target  volume  range  (m) 

D|.  - receiver  to  target  volume  range  (m) 

Lp(t)-  propagation  loss  over  transmitter-to*target  path 
Lp(r)-  propagation  loss  over  receiver-to-target  path 
Lg  - system  losses 


For  the  system  described  in  II  above,  coherency  is  required 
between  the  two  transmitted  signals  and  this  will  be  achieved  by  splitting 
the  output  of  the  RF  generator  and  feeding  separate  antennas  over  identical 
wave  guide  paths.  Thus  the  P^  term  and  the  G^  term  will  include  both 

transmitting  beams.  Similarly  Lip(t)  will  apply  to  both,  on  the  assumption 

that  the  medium  will  be  statistically  homogeneous  and  isotropic  for  the 
geometry  selected.  These  loss  terms  Lp(t)  and  Lp(r)  will  include 

radome  losses.  The  selection  of  a radome  material  will  be  covered 
subsequently. 


Selection  of  range  to  the  scattering  volume  V will  be  determined  by 
its  dimensions,  antenna  apertures  and  far-field  requirements.  These  in 
turn  will  dictate  P^  to  achieve  a Py  at  least  6 dB  above  the  system  noise. 

Each  parameter  will  be  examined  in  turn,  for  each  of  two  systems,  one 
operating  in  the  X-band  at  a frequency  of  10.  525  GHz,  and  one  at  16  GHz. 
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A.  X-Band  System  Design  Parameters 


A required  scattering  volume  of  a 10  cm  cube  has  been  defined. 

This  will  affect  the  beamwidth  and  the  far  field  requirements  which  will 
be  addressed  initially.  The  tunnel  is  assumed  to  be  7.  2 ft  across,  or 
2.2  m.  If  the  intersection  angle  of  the  two  transmit  beams  is  #,  referring 
to  Fig.BI,  we  will  choose  $ = 120°.  Then  with  W = 2.  2 m and  beam 
crossing  at  Dj.  = 1.  1 m,  the  separation  between  the  two  transmit  antennas 
is  2s  = 3.81  m,  and  Dt=  2.2  m.  A half  power  beamwidth  of  ^ radians 

for  a uniformly  illuminated  aperture  D in  the  far  field  (at  = 

would  require  D = X = 22  X or  . 627  m.  But  this  would  mean 


that  the  far  field  criterion  could  not  be  met. 


20^ 

since 


for  D = 


. 627  is 


27.6  m.  We  are  thus  constrained  to  a near  field  analysis  of  the  situation, 
since  the  receive  range  to  the  scattering  volume  is  one  half  the  transmit 
range. 


For  a typical  antenna  with  gain  the  power  density  at  a 


distance  R in  the  far  field  is  given  by 

PtGt 


Wt  = 


The  half  power  beamwidth  in  degrees  is  given  by  69  when  a cosine 
illumination  function  is  employed.  Assuming  this  in  each  dimension, 
the  gain  is  given  by  27000/(69A/D)^> 

Therefore  G,  = X 103)(.  627)^  ^ 3,^ 

(69)2{. 0285)2 


and 

Pt(274) 

Wt  = -i 

(4ir)R2 

If  we  let  R = = 27.  6 m,  then  Wt  = . 0286  P^ 

From  Ref.  9,  curve  No.  34  (power  density  in  the  near  field  normalized 

to  unity  at  2D^/ X , when  the  normalized  range  R^  = ■'  = 2.:,^—  = 

2D2/A  27.6 

.0797,  the  power  density  increases  by  12  dB,  or  15.85  times.  Therefore 
Wt  = . 453  Pt  in  W/m^. 
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If  now  we  consider  the  scattering  cross-section  as  that  due  to  the 
scattering  volume  common  to  two  pencil  beams  of  10  cm  diameter, 
resulting  from  focusing  the  antenna  radiation  at  the  desired  range,  a 
common  volume  such  as  indicated  in  Fig.BS  results.  An  approximation 
which  may  be  used  for  the  volume  is  the  area  of  the  parallelogram  times 
the  width  of  the  half  power  beam,  which  is  a volume  of  .001  m^. 

Multiplying  this  times  the  particle  volume  density  and  the  individual 
particle  cross-section  we  obtain  a scattering  cross-section 

V a y = (.001)(11.34  X 10*®)  = 11.34  X 10*  ^^m^ 

The  power  which  is  isotropically  scattered  is  then 
WjV  ffy  = 5.  14  X lO-llP^  in  watte 


At  a range  R the  scattered  power  density  is 


W 


r 


5.  14  X lO-llp^ 

47rR2 


If  the  range  is  = 1.1m,  we  have 
5.  14  X 10-^lp. 

W = j = . 338  X 10*  ^Ip, 

4ir(i.  1) 


If  the  receiving  aperture  is  made  to  match  the  scattering  volume 
dimensions,  a sectoral  horn  of  width  20  cm  and  height  10  cm  is  indicated. 
Using  this  aperture  with  an  effectiveness  ratio  of  90%  we  obtain 


Pr  = W,A^  = (.338  X 10*^^P^)(.02)(.9) 
= 6.08  X lO-^^Pj 


The  loss  over  the  propagation  paths  is  given  by 


Lp(t)(dB/m)=  4.34  X 10“* 


[ 


ir2 

X 


( s,d’) 


Im(-K) 


where  the  summation  is  over  a Im^  volume,  D is  particle  diameter  in 
centimeters,  A is  in  centimeters,  I|^(-K)  is  the  imaginary  part  of  -K, 
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and  K depends  on  the  dielectric  constant  of  the  particles  as  described 
in  IVB2  above.  For  the  dust  we  can  assume  the  conductivity  = 

.0006  mho/m.  Then 


= 2.  55 


.0006(36ir  X 10^) 
(2t)(10.525  X 10^) 


= 2.  55  - j 103  X 10"® 


,,  m2  - 1 1.55  - j(103  X 10"^) 

m2  + 1 3.55  - j (103  X lO"^) 

K = .676  - j(2.53  X 10"‘*)Ij„(-K)  = 2.  53  X 10-4 


For  the  two  largest  sized  particles: 

Lp(t)  = (4.34  X 10-4)  — 2_  ^2.  53  x 10-4)[(5.9  X 104)(8.4  X 10"®) 

2.  85 

+ (23.  6 X 104)(4.2  X lO”®)] 

= 5.  66  X 10-9  dB/m 


which  can  be  neglected,  as  can  Lp(r).  Lg,  or  system  losses,  will  be 
estimated  at  2 dB. 


Pp  - A signal  to  noise  ratio  of  10  dB  as  a minimum  will  be  desired 
for  the  system.  This  means  that 

Pr  = 10  KToB„F„ 

We  will  assume  a noise  figure  F^  of  12  dB  (15.84);  Tq  = 3000°K  + 100°K  + 
4303°K,  and  B^j  - 120  kHz.  The  selection  of  Bjj  is  based  on 

- 2Vrsin(g/2)  _ 2(1828.  8)sin  60°  =nikHz. 

A .0285 

From  this  KTB  = (1.38  X 10-23)(7403.  6) ( 12  X 104)  = 1.22  X 10"^4w 

andP].  = 1.22  X 10”^^  W = -129  dBW  with  10  dB  signal  to  noise 
(S/N)  ratio. 
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- We  can  now  find  the  required  transmitter  power,  knowing  the 


required  receive  power.  Listing  the  values  of  the  parameters  as  obtained 
previously  we  have: 


r — ’ ” 

1. 22  X 10-^^W 

X-Band 

Ku-Band  (Sec.  C following) 

Pr 

-129  dBW 

-130  dBW 

Gt 

274 

24  dB 

31  dB 

Gr 

278 

24  dB 

31  dB 

A2 

(.0285)2m2 

-30.9  dBm' 

^ -34. 7dBm2 

11. 34  X 10-8m2/m3 

-69.5  dBm 

-61.9  dBm 

V 

.001  m^ 

-30  dBm^ 

-30  dBm^ 

Dt^ 

(2.2)  ^m^ 

6.  8 dBm^ 

2.  1 dBm2 

Dr^ 

(1.  I)^m2 

. 8 dBm^ 

. 8 dBm^ 

Ls 

1.  58 

2 dB 

2 dB 

(4W)^ 

1984.4 

33  dB 

33  dB 

To  find  the  required 

transmitter  power 

it  only  remains  to  set 

the  two  expressions  for  received  power  equal,  and  add  2 dB  for  system 
loss  to  the  result. 

6.08  X 10‘^^P^  = 1.22  X 10‘^3 
Pt  = 2 W 

After  allowing  for  2 dB  of  wave  guide  and  radome  loss,  3 W of  average 
power  appears  to  be  adequate  for  the  job. 

B.  Proposed  X-Band  Measurement  System 

A principal  consideration  of  the  system  is  its  mechanical  layout 
and  installation.  It  must  be  protected  from  the  effects  of  the  explosion, 
and  it  must  be  isolated  from  vibration  and  mechanical  shock  effects, 
while  still  being  rigidly  interconnected  by  wave  guide.  One  approach 
could  be  to  build  the  system  on  a wooden  pallet  which  can  be  installed 
in  an  accessible  cavity  beneath  the  tunnel  floor.  FigureBS  illustrates 
the  location  and  mechanical  layout.  The  microwave  equipment  could 
be  rigidly  mounted  to  a 3/4  in.  plywood  pallet,  which  could  be  installed 
on  a vibration  isolating  base  14  ft  X 3 ft  and  12  in.  thick.  An  excavation 
in  the  shock  tunnel  floor  would  be  prepared  into  which  the  pallet  could 
be  placed  with  provision  for  120  V 60  Hz  power  input  and  signal  cable 
output.  The  system  would  then  be  covered  with  a reinforced  and  braced 
1/2  in.  steel  plate  into  which  might  be  mounted  the  slip-cast, fused-silica 
radomes. 
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Figure  B6.  Microwave  Velocimeter,  Tunnel  Floor  Installation  Layout, 
Typical  for  120®  Bean  Intersection  at  10.525  GHz 
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C.  Recommended  Ku-Band  System 

1.  For  several  reasons  a Ku-band  (15-16  GHz)  system  is 
recommended.  These  are: 

a.  Scattering  cross-section  increases  as  1/A^  in  the 
Rayleigh  region. 

b.  Plasma  cut-off  frequency  and  skin  depth  effects 
will  be  minimized. 

c.  Off-the-shelf  antennas  are  available  at  reduced 
cost- -adequate  power  can  be  obtained  at 
reasonable  cost. 

d.  Semirigid  coaxial  cable  can  still  be  used  rather 
than  wave  guide  runs  to  minimize  vibration  effects. 

2.  The  following  list  of  system  components  and  specifications 
with  estimated  cost  figures  is  recommended  for  a feasibility  model  of 
a microwave  velocimeter.  The  frequency  of  operation  will  be  16.2  GHz 
or  A = 1.85  cm. 


Transmitter 

Solid  State  16  GHz  Signal  Source  VSU  9002  $ 410.00 

M-1035-00 

16.2^  100  MHz  Mech.  Tuneable 

150  mW  min.  output 

10  V P.S.  Heat  Sink  Cooling 

Antennas  (three) 

TRG  (Control  Data)  856-12  Collimating  Lens  Corrected 
Conical  Horns 
Aperture  Size  12  in. 

Gain  31  dB 

1st  side  lobe  E-Plane  -25  dB,  H-Plane  -17  dB 

HPBW4.70  $ 1,450.00 


Receiver 


Mixer 

Anaren  Model  7A129  (ortho-quad)  Balanced  Mixer 
VSWR  1.7 

Isolation  LO/RF  20  dB 
Conversion  Loss  9 dB 
Noise  Figure  9 dB 

IF  Bandwidth  dc  - 600  MHz  $ 355.00 
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Preamplifier  ARL  Designed  and  Built  $ 500.00 

Bandwidth  0-185  kHz 
Gain  40  dB  NF  10  dB 

High  Gain  Low  Noise  Amplifier  $ 600.00 

80  dB  Gain  NF  10  dB  300  MHz  B.W. 


An  alternative  transmitter  source  is  a Varian  two-cavity  klystron 
which  produces  5 W of  CW  power.  Price  is  $2420.  • This  would  require 
a klystron  power  supply,  at  an  estimated  cost  of  $4000. 


Wave  Guide  Connectors,  Wave  Guide 

Flanges,  Miscellaneous  Hardware 

Radomes  (3  each  slipcast  fused  silica)  Total 

3 dB  Power  Divider  (4  port  power  splitter)  Ku-600 

30  dB  directional  coupler  559-30 

Components  and  Parts  (using  solid  state  source) 


$ 950.00 
$ 4, iOO.OO 
$ 495.00 
$ 525.00 
$12, 360.00 


A block  diagram  of  such  a system  is  presented  in  Fig.B7.  The 
angle  between  the  transmitters  is  reduced  to  60°,  with  the  range  to  the 
scattering  volume  maintaitied  at  1.  I m. 

3.  Design  Analysis  - Based  on  a wavelength  of  A = 1.85  cm  the 
various  bistatic  radar  parameters  for  the  rccomntended  system  are 
determined  as  follows; 


ai  = 

4.  8 X = 

\1.85/ 

2.  7 X 10" ‘^m 

Ov  = 

(ll.  34  X l0-®)(5.632) 

= 6.  39  X lO""^ 

= -61.9  dB 

Dr 

1.  l m 

= . 8 dBm^ 

Dt 

1 . 27  m 

Dt^ 

= 2.  1 dBm^ 

At  a range  of  1.  27  m the  4.7°  HPBW  of  the  selected  antennas 
subtends  10  cm.  This  results  in  an  effective  scattering  volume  of 
approximately  I0"^m^,  or  -30  dB. 

Gt  31  dB 

G^  31  dB 
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Roof  of  Tunnel 


XMIT  Antenn 


Gunn  Oscillator  150  mW 


l^(t)  L|^r)  These  losses  ere  negligible  in  the  medium.  The  slip-cest, 

fused  siiice  redomes  ere  IS  in.  squere  milled  meteriel  (eveileble  from 
Brunswick  Corp. , Merion,  VA)  7 wevelengths  in  thickness,  or  ebout 
2.9  in.  in  thickness.  The  loss  tengent  for  this  meteriel  is  given  es 
.003  end  the  dielectric  constent  is  3.2.  The  thickness  is  given  by 

nA 

2(«,-sin^#)^^^ 

where  # is  meesured  from  the  normei  to  the  redome  surfece.  If  # is 
0®  end  d s 7.  24  cm  then 

„ = ILilUMJLli  = 14 

1.8S 


The  tensile  strength  is  given  es  4000  - SOOO  psi  with  flex  strength  of 
4S00  pel  end  modulus  of  elesticity  of  3.  S X 10^.  The  meteriel  porosity 
cen  be  reduced  by  impregneting  with  silicone  resin  which  reduces  its 
susceptibility  to  moisture,  greese.  end  other  conteminents.  From  Ref. 

10.  e formule  for  reflection  coefficient  0 is  given  es 

0 8 - exp(-24Lo  - 2{d)] 

I • r^  exp(-2#Lo  • 2jd) 

^ - sin^#y^^  end  # s incidence  engle  meesured  from  normei. 

If  y = 7,  end  # = 30®  ♦ = 14*(3.2  - .25)*^^  = 75.  54  rediens 
L is  the  ettenuetion  due  to  conductive  losses 


.-SLilLooJL  ^ 

.2i 


2(«.  sin^d)  2(3. 2 -.25) 


r is  the  interfece  reflection  coefficient  given  by 


«j.  - sin^# 

»nd  2 = 3.93  for  vertical  (perpendicular)  polarization. 

cos  # 

is  the  reactance  due  to  conductive  loss: 

= = 3.2  X 10-3 

2.974 

-.98  + j(3.2  X 10-3) 

1 + V3793-  j (3.2  X 10-3)  2.98- jP. 2 X 10-3) 


r 

2cos#[  y-sin^#J^^^ 

1 - /TTgVi  j(3.2X  10-3) 


r-i-.33 


Then 

-.33(1  - ^7)(16.2  X Kr^)^ 

^ I - (..33)2  X 10--*)  . ^.2j(75.54) 

|P(^  = .0069 


= -.054  - j .067 


The  transmission  at  this  angle  is  then  99%  and  losses  can  be 
considered  negligible  for  both  Lp(t)  and  Lp(r).  The  radome  size 
seiection  was  based  on  radome  installation  parallel  to  the  antenna 
aperture.  If  the  radomes  are  installed  normal  to  the  direction  of 
propagation  no  losses  would  result.  A loss  of  . 25  dB  will  be  estimated 
for  Lp(t)  and  Lp(r). 

Lg  will  be  estimated  at  2 dB 

Pt  = 5 W = 7 dBW 

A = 1.85  cm  = -34.  7 dBm^ 

Finally  we  obtain  in  dBW: 

P,  = 7 + 31  + 31  + (-34.7)  + (-61.9)  + (-30)  - 33  - 2.  1 - . 8 - . 25  -.  25  -2 
- -96  dBW 
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To  calculate  the  S/N  ratio  we  again  need  the  effective  bandwidth: 

B = fn  = 2(1828.  8)_sin30O  ^ ^ jqq 

n u 

.0185 

The  noise  power  then  is  kTB  = Nq 

Nq  = (1.  38  X 10-23)(7403.6)(105)  = 1.02  X lO-^^ 

Nq  = -140  dBW 

The  signal  to  noise  ratio  (S/N)  is  the  ratio  of  to  N^: 

S/N  = -96  + 140  = 44  dB.  In  other  words,  there  is  a 34  dB  advantage  in 
S/N  ratio,  over  the  X-band  design. 

4.  Mechanical  Layout  - The  system  will  be  installed  on  a pallet 
6 ft  long  by  2 ft  deep,  so  that  it  may  be  positioned  adjacent  to  the  tunnel 
wall  with  the  conical  horn  antennas  inserted  into  openings  preformed 

in  the  cement  wall  at  1/2  the  height  of  the  tunnel. 

5.  System  Engineering  - In  addition  to  the  required  components 
and  parts,  an  estimated  12  man-months  of  system  engineering,  at  a cost 
of  $42,000  would  be  required.  The  total  estimated  cost  of  the  system  is 
thus  about  $54,  360. 

In  the  interest  of  cost  effectiveness,  it  is  suggested  that  the 
transmitter  power  may  be  adequate  using  a 150  mW  solid  state  source. 
The  estimated  30  dB  excess  S/N  ratio  could  eliminate  the  requirement 
for  the  klystron  source  or  for  a TWT  amplifier.  The  system  could 
readily  be  modified  to  add  the  TWT  amplifier,  and  in  fact  could  be 
designed  with  the  add  in  capability,  should  additional  sensitivity  be 
required.  The  cost  of  a klystron  power  supply,  plus  the  klystron 
source,  which  is  heavy  and  bulky,  makes  it  a less  attractive  solution 
for  the  transmitter  requirement. 

6.  Millimeter  Wave  System  - At  96  GHz  the  wavelength  is 
approximately  one- sixth  that  at  the  selected  Ku-band  frequency.  The 
scattering  cross- section  would  be  improved,  therefore,  by  31  dB. 

This  could  make  the  system  design  at  this  wavelength  more  attractive, 
since  the  antenna  gains  might  be  12  to  15  dB  higher  also.  Unfortunately, 
a klystron  and  power  supply  would  be  required,  costing  about  $8500,  just 
to  get  250  mW  of  power.  In  addition,  millimeter  components  are 
substantially  more  expensive.  For  example,  a 30  dB  directional  coupler 
costs  $675,  and  a power  divider  costs  $625,  and  of  course  wave  guide 
and  other  components  would  be  somewhat  higher. 
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INTRODUCTION 


It  has  been  proposed  that  a "blast  gauge"  be  used  to  measure 
local  gas  velocity  and  speed  of  sound  during  teats  of  the  Dynamic  Air 
Blast  Simulator.  The  principle  of  the  device  is  to  generate  at  some 
point  in  the  simulator  a pressure  pulse  whose  amplitude  is  sufficiently 
large  as  to  be  distinguishable  in  the  noisy,  nonsteady  signal  generated 
by  the  primary  blast  wave,  and  to  measure  the  velocity  of  propagation 
of  this  pulse  between  two  relatively  closely  spaced  detectors. 

The  theory  for  the  propagation  of  a finite -amplitude  pulse  through 
a nonuniform,  nonsteady  flow  field  is  very  complicated.  In  order  to  make 
the  blast  gauge  a feasible  instrument,  significant  simplifications  in  the 
theory  of  its  operation  must  be  made.  In  this  report  we  show  what 
assumptions  have  to  be  made  and  point  out  some  of  the  features  of  the 
technique. 


THE  BLAST  GAUGE 

Consider  the  configuration  shown  in  figure  Cl  (for  simplicity,  a two 
dimensional  picture  is  shown;  the  formulas  can  easily  be  goneraliKed  to 
three  dimensions).  A charge  of  high  explosive  is  located  at  the  origin 
and  is  immersed  in  the  one-dimensional  flow  field  u(x,  t)  generated  by 
the  primary  blast  in  the  Air  Blast  Simulator.  Three  pressure  transdxicers 
are  located  in  an  array  at  points  (1),  (2)  and  (3),  point  (1)  being  distant 
r from  the  charge. 

At  some  time  (t=0)  after  the  passage  of  the  primary  shock  the  charge 
is  detonated.  The  blast  wave  from  this  explosion  propagates  relative  to 
the  nonuniform,  nonsteady  moving  fluid  with  velocity  c.  The  wave  inter- 
acts with  the  gradients  in  the  flviid,  and  its  strength  (e.  g. , its  Mach 
number  M~c/a,  where  a is  the  local  sound  speed)  changes,  in  the  same 
manner  as  when  a wave  interacts  with  an  area  change.  We  are  primarily 
interested  in  the  velocity  of  the  wave,  and,  if  the  wave  is  weak,  the 
variation  of  M due  to  the  interaction  with  nonuniformities  is  second  order, 
so  the  effect  on  c is  negligible.  However,  other  causes  of  variation  of 
the  velocity  c may  be  first  order.  For  example,  if  the  variation  of  wave 
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Charge 


Figure  Cl.  Geometry  of  Blast  Gage 
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strength  M along  the  path  length  R(t)  due  to  spherical  spreading  is  known 
a o r iori  . then  M (r)  can  be  accounted  for.  Furthermore,  in  general 
the  sound  speed  is  variable,  a (x,  t),  and  it  may  be  necessary  to  account 
for  this. 

The  geometry  of  shock  wave  propagation  in  multidimensional  flows 
can  most  easily  be  studied  according  to  the  theory  of  nonlinear  shock 
dynamics  of  Whitham.  For  weak  shocks  this  reduces  to  the  theory  of 
nonlinear  geometrical  acoustics  of  J.  B.Kellor,  which  has  been  worked 
out  in  much  greater  detail.  For  waves  of  infinitesimal  strength  it  is  the 
familiar  theory  of  geometrical  acoustics.  In  all  cases  the  wavefronts 
propagate  along  trajectories  known  as  rays.  In  our  problem  (figure C2) 
the  waves  propagate  relative  to  the  fluid  in  a direction  normal  to  the 
front  with  speed  c and  are  convected  along  with  the  fluid  velocity  u(x,  t). 
The  resultant  wave  speed  v is  the  vector  sxim  of  c_n  and  u,  where  is 
the  \mit  vector  normal  to  the  wavefront,  v is  tangent  to  the  ray  at  the 
wavefront.  Thus,  in  a moving  fluid  the  rays  are  not  normal  to  the  wave> 
fronts,  and  the  wave  speed  is  not  c. 

RAY  TRACING 

In  nonuniform,  nonsteady  flow  u(x,  t)  and  a(x,  t)  are  variable,  and, 
if  the  wave  strength  is  decaying,  M(r)  is  also  variable,  so  in  general  the 
rays  are  curved.  (Furthermore,  nonlinear  effects  cause  ray  curvature, 
but  we  are  neglecting  them  here.  ) Therefore,  in  order  to  determine  what 
part  of  the  wave  front  is  being  observed  at  any  point  in  space  and  time  it 
is  necessary  to  trace  rays  by  solving  the  vector  equation  for  the  ray 
trajectory 


dr 


= V, 


(Cla) 


or 


dx 

dt 


V 

X 


(Clb) 
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Figure  C2.  Propagation  of  Secondary  Shock 
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where,  from  figure  C2 


v^=  M (R)  a (x,  t)  cos  0 + u(x,  t),  (C2a) 

Vy  = M (R)  a (x,  t)  sin  0 (C2b) 

In  C2  we  have  used  c = Ma.  6 is  constant  along  a ray  and  is  the  label  denoting 
the  part  of  the  wavefront  that  is  being  observed.  Thus,  the  ray  emanating 
from  the  explosion  with  angle  0 can  be  constructed  by  integrating  the 
equations 


dx  = [M^R(t)j  a(x*t)  cos  0 + u(x,  t)  ]dt  (C3a) 

dy  = M(R(t))  a(x,t)  sin  0 dt  (C3b) 


to  give  the  parametric  representation  x(t;3),  y(t;0)  for  the  ray.  The 
path  length  R (t)  is  evaluated  from 

R(t;0)  =J  v{t;e)dt,  (C4) 


where 

v(t;0)  =|[M(R(t))  a(x(t),  t)J  + |^u(x  (t),  tjj]  + 2M^R(t)J  a^x(t),ty  u(x(t),  t^cosO^ 

^ (C5) 

2 2 2 

which  is  simply  v = c + u + 2uc  cos  0,  with  the  dependence  on  the  inde- 
pendent variables  written  out  explicitly.  It  is  clear  that,  in  the  general 
case,  C3  can  be  solved  for  the  ray  coordinates  together  with  the  integral 
equation  C4  ^or  the  path  length  only  by  an  approximate,  iterative  scheme, 
and  this  only  under  circiimstances  when  good  estimates  of  the  entire  flow 
field  u(x,  t) , a(x,  t)  are  available. 

Equations  C3  and  C4  would  in  general  be  used  with  the  blast  gauge 
as  follows:  Using  a preliminary  estimate  of  the  flow  field,  rays  are  traced 
from  the  explosion  to  the  locations  of  the  detectors  and  arrival  times  are 
thereby  predicted.  The  predictions  are  compared  with  measurements,  and 
any  differences  arc  used  to  correct  the  estimate  of  the  flow  field,  presum- 
ably by  some  kind  of  optimal  interpolation  procedure.  The  process  is  repeated 
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until  agreement  is  obtained.  This  is  clearly  a very  complex  and 
sophisticated  process,  particularly  if  several  different  detectors  and 
several  different  explosions  are  used  for  the  purpose  of  providing 
sufficient  data  to  adequately  define  the  flow  field.  Thus,  in  order  that 
the  blast  gauge  be  an  easily  applied  tool,  it  is  essential  that  further 
simplifications  be  incorporated. 

SIMPLIFIED  TECHNIQUE 

If  the  detectors  are  located  close  enough  together  that  Ax,  Ay  « L, 
where  L is  some  length  scale  of  the  flow  field  (given  by  the  gradients 
behind  the  primary  shock;  e.  g. , of  order  of  the  driver  length)  and  if 
Ax, Ay  « r,  then  locally,  in  the  triangle  of  the  detectors  the  waves  can 
be  treated  as  being  plane,  the  rays  as  straight  and  the  flow  as  uniform 
and  steady.  In  figure  C3  the  wav^  front  is  shown  at  three  different  times 
as  it  passes  the  three  detectors.  The  time  interval  measured  between 
the  first  pair  of  detectors  is  ^2  ' ^1  ^ ^1  between  the  second  pair 
is  t^  - t^  -^2*  figure  shows  the  distance  traveled  by  the  wave 

along  the  vector  v in  these  intervals.  From  the  geometry  shown  in 
the  figure  (cf.  dotted  line),  it  follows  that 


(c  + u cos  6)  At, 

■ — = cos  0 

(C6a) 

(c  + u cos  0)  At, 

, i . sin  a 

Uy 

(C6b) 

a - 0. 

(C6c) 

Therefore, 


(C7a) 

(C7b) 
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Figure  C3.  Propagation  Direction 


Equations  C7  are  the  fundamental  relations  for  the  blast  gauge.  They 
show  that: 

1)  It  is  necessary  to  have  an  array  of  three  detectors  as  shown 
in  order  to  determine  both  the  angle  6 and  the  velocities « and 

2)  It  is  necessary  to  observe  two  different  waves  (i.  e. , with 

2 different  angles  6^.  62)  in  order  to  resolve  separately  the  velocities 
u and  c.  That  is,  if  two  waves  are  observed,  closely  enough  spaced 
in  time  that  the  flow  field  does  not  change  significantly,  then  the  two 
measurements  yield,  say. 


so  that 


c -i-u  cos  0^  = Kj 
c + u cos  62  = ^2  » 


u = 


cos  0 j - cos  &2 


c = 


Kj  cos  02  - K2  cos  0j 
cos  @2  ■ cos  6 j 


(C8a) 

(C8b) 


(C9a) 

(C9b) 


On  the  other  hand,  if  it  is  known  "where  the  wave  came  from" 
(i.  e. , the  angle  Cp ),  then  another  relationship  between  c and  u is 
available.  <p  is  known  if  the  ray  r from  the  explosion  to  the  detector 
array  can  be  considered  to  be  straight,  i.  e.  , if  r « L.  Therefore, 
providing 


Ax,  Ay«r«L  (CIO) 

c and  u can  be  evaluated  from  one  wave  only.  For,  from  figure  C3 

V cos  (0  -Cp)  r:  c + U cos  0 (Cll  ) 

Expressing  v in  terms  of  c and  u,  and  using  C7b  to  evaluate  the 
right-hand- side. 
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(C12) 


c + u + 2uc 


c..9=[^^ 


COS  9 ’ 

cos  (6  - tp). 


Cl 2 can  be  used  together  with  equations  C7  to  evaluate  c and  u; 


(C7a) 

(C13a) 

(C13b) 


Note,  however,  that  for  6 = cp  r 0 the  additional  information  is  lost;  that 
is,  the  array  must  not  be  directly  downstream  from  the  explosive  source. 

For  "slight  offset"  (0«1  or  TT-0«I)  simplified  equations  result. 
First  consider  downstream-propagating  waves  (6«1).  Then  equation  C11 
becomes 


cp  ;= 


e 


and  C7  become 


0 •= 


Ax 

w 


“2 


. . .Ax 
c + u 


(C14) 


(C15a) 


(Cl  5b) 


Thus,  for  a downstream  wave  with  the  gauge  slightly  offset,  c and  u can 
be  evaluated  independently: 


(C16a) 
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(C16b) 


u - 


Ax 

At 


cp 


1 


AU 


On  the  other  hand,  for  upstream  propagating  waves  with  slight 
offset  (TT-0«  1),  Cn  becomes 


cp=20-n  ; c>u 


(Cl  7a) 


cp  = i{Q  -TT)  ; c<  u . 


(C17b) 


That  is,  for  c > u cp  must  be  almost  tt  while  for  c<  u cp  is  almost  zero, 
and  to  first  order  an  additional  relation  between  u and  c is  not  obtained.' 
In  this  case  C7  become. 


e=i- 

Ay 


^2 


(CliM) 


c 


(Cibb) 


In  summary,  with  a blast  gauge  satisfying  the  conditions  CIO  and 
with  the  array  set  downstream  from  the  explosive  source  with  slight  off- 
set, if  the  flow  locally  supersonic  (c<u),  first  the  downstream 
propagating  wave  passes  the  array,  yielding  c and  u according  to  CIS  , 
and  then  the  upstream  propagating  wave  passes, yielding  6 by  C17b  or 
C18a»  and  c-u  by  C18b  . In  subsonic  flow  (c>u)  the  same  redundant 
information  may  be  obtained  with  a blast  gauge  if  an  additional  array  of 
detectors  is  set  upstream  of  the  source. 

It  shoxild  be  noted  that  if  the  time  of  ignition  of  the  blast  (t=0) 
relative  to  the  time  of  arrival  at  point  (1)  (i.  e.  , if  the  time  t^)  is  accu- 
rately known,  if  r is  taken  as  Ax  and  if  M(R)  is  known  all  the  way  from 
the  soui'ce,  then  the  occurence  of  the  blast  can  be  taken  as  one  of  the 
times  of  arrival,  and  the  number  of  transducers  required  can  be  reduced 
by  one. 
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1 . BACKGROUND 


Measurement  of  alrblast  flow  velocity  In  a DABS  test  Is 

essential  for  adequate  definition  of  the  flow  environment, 

but  there  are  difficult  technical  problems  to  overcome.  The 

flow  velocity  Is  too  high  and  the  blast  environment  Is  too 

severe  for  conventional  anemometers  to  work.  Optical 

methods  for  measuring  flow  velocity  are  difficult  to  begin 

with,  and  are  Impractical  or  Impossible  when  account  Is 

taken  of  smoky  Impurities  In  the  explosive  gases  together 

with  problems  of  optical  alignment  In  an  explosive  field 

test.  Doppler  shift  radar  Is  a possibility,  but  depends  on 

adequately  launching  lightweight  reflecting  materials  In 

the  alrblast  flow  or  utilizing  Interface  conductance 

- * • • • 

phenomena.  Even  If  successful,  the  radar  method  only  deter- 
mines the  velocity  of  a particular  element  of  fluid,  and 
does  not  yield  directly  the  desired  result  of  gas  velocity 
history  at  a given  test  station. 

On  the  other  hand  measurement  of  flow  velocity  In 
explosively  driven  plasma  for  MHD  electric  power  applications 
Is  outstandingly  simple,  accurate,  reliable,  and  compatible 
with  field  test  conditions.  Standard  practice  at  Artec  In 
these  tests  Is  to  provide  a magnetic  field,  either  with  a 
permanent  magnet  or  an  electromagnet,  across  the  flowing 
gaa,  and  pick  up  the  Faraday  generator  voltage  with  simple 
electrodes  fabricated  from  brass  screws.  The  electrode 
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voltage,  as  measured  directly  across  the  electrodes  through 
a long  (100  ft)  terminated  cable.  Is  directly  proportional 
to  the  gas  flow  velocity.  Knowing  magnetic  field  strength 
the  relationship  between  velocity  and  voltage  Is  .iinown 
absolutely,  and  there  Is  no  need  for  gage  calibration. 

The  magnetic  velocity  gage  Is  an  outstanding  candidate 
for  DABS  tests  In  terms  of  sln^llclty  and  practicability, 
but  there  are  serious  problems  associated  with  low  gas  con- 
ductivity which  must  first  be  overcome.  In  explosive  MHD 
experiments  the  gas  conductivity  Is  high,  and  the  Internal 
resistance  of  the  Faraday  generator  Is  only  a few  mllllohms. 
By  contrast  the  theoretical  equilibrium  conductivity  of 
seeded  detonation  products  Is  so  low  that  Internal  resistance 
Is  many  megohms.  Nonequilibrium  effects  In  the  expanding  gas 
undoubtedly  Increase  the  conductivity,  but  the  generator 
resistance  Is  still  so  high  that  amplifiers  with  high  Input 
Impedance  are  required.  Measurement  of  generator  voltage 
at  high  Impedance  levels  Introduces  concerns  with  noise  and 
static  electricity  that  do  not  affect  explosive  MHD  measure-  . 
ments. 

The  present  program  was  Initiated  to  evaluate  the 
utility  of  the  magnetic  field  gage  In  measuring  dynamic 
flew  environment  In  DABS  tests.  The  program  Is  divided  Into 
two  phases:  Phase  I consists  of  gage  design  and  preliminary 
evaluation  tests  conducted  by  Artec  at  a local  test  site. 
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and  Phase  II  consists  of  evaluation  of  the  gage  In  scheduled 
large  scale  DABS  tests. 

This  report  describes  the  results  of  Phase  I tests  to 
date.  It  Is  concluded  that  problems  associated  with  high 
Inqpedance  measurements,  stray  capacitance,  and  charge 
accumulation  effects  preclude  advancing  to  Phase  II  at  this 
time.  The  data  are  not  conclusive  with  respect  to  the 
feasibility  of  the  concept,  and  there  Is  a reasonable  possi- 
bility that  the  measurement  problems  can  be  resolved  with 
further  Phase  I effort. 
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2.  INTRODUCTION  TO  MHD  VELOCITY  MEASUREMENT 
2 . 1 Faraday  Generator 

Michael  Faraday  originally  observed  that  a voltage  Is 
generated  when  a conductor  moves  through  a magnetic  field. 
Maximum  voltage  Is  obtained  when  the  conductor  motion,  the 
magnetic  field,  and  the  pickup  electrodes  are  perpendicular. 
.When  used  as  a velocity  gage  the  Faraday  generator  consists 
of  an  electromagnet  field  coll  (or  permanent  ma^iet)  and  a 
pair  of  electrodes  as  shown  In  Figure  D1. 

Open  circuit  voltage  across  the  electrodes  Is 

(01) 

ji 

where  B Is  the  magnetic  field  In  tesla  (1  tesla  ■ 10  gauss), 
u Is  the  gas  velocity  In  meters/sec,  A Is  the  Int ere le strode 
spacing  In  meters,  and  V Is  the  open  circuit  voltage  In 
volts.  This  equation  relates  the  open  circuit  output  voltage 
to  the  gas  velocity  In  terms  of  two  readily  determined 
parameters  B and  A . 

When  the  generator  Is  attached  to  a load  there  Is  a 
voltage  drop  due  to  resistive  losses  In  the  gas.  In  that 
case  the  generator  voltage  Is  given  by 

V'*  BuJA  - R1  (02) 
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The  equivalent  circuit  of  the  Faraday  generator  consists  of 
a voltage  source  In  series  with  the  Internal  gas  resistance 

R . 

Oas  resistance  Is  related  to  gas  conductivity  by  the 
equation 


(“) 

where  J Is  the  Interelectrode  spacing  In  meters,  6 Is  the 
gas  conductivity  In  mho/m,  A Is  the  electrode  area  In 
square  meters,  and  /(*  Is  a geometric  correction  factor  which 
Is  approximately  unity  for  large  electrodes.  The  correction 
factor  Is  usually  determined  by  direct  experiment  using  a 
convenient  conductor  (e.g.  copper  sulfate  solution). 

If  the  Input  Impedance  of  the  oscilloscope  or  Instru- 
mentation amplifier  Is  much  larger  than  the  Internal  resis- 
tance of  the  gage  the  resistive  voltage  drop  In  the  generator 
Is  negligible  and  open  circuit  voltage  Is  measured.  In  that 
case  the  measured  voltage  Is  proportional  to  Velocity,  and 
It  does  not  depend  on  the  value  of  the  gas  conductivity. 

A necessary  condition  for  the  velocity  gage  to  work  Is 
therefore  that  the  gas  resistance  be  much  smaller  than  the 
Input  Impedance  of  the  Instrumentation  amplifier.  In 
explosive  NHD  studies  the  Internal  gas  resistance  Is  about  a 
mllllohm,  and  In  combustion  MHD  for  large  power  Installations 
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the  Internal  resistance  Is  a fraction  of  an  ohm.  In  these 
cases  there  Is  no  problem  whatsoever  In  measuring  open 
circuit  voltages.  In  the  magnetic  velocity  gage  for  DABS, 
on  the  other  hand,  gas  resistance  Is  very  high,  and  there  Is 
considerable  difficulty  assuring  that  the  condition  Is  met. 
Furthermore,  as  will  be  discussed  later,  when  all  the  signal 
processing  occurs  at  very  high  Impedance  levels,  static  charge 
and  capacltatlve  coupling  become  circuit  design  concerns  not 
present  In  the  MHD  work. 

2.2  Qas  Conductivity 

A DABS  blast  pulse  Is  produced  by  detonating  an  explosive 
In  a confined  earthen  tube.  At  the  test  station  the  first 
portion  of  the  pulse  Is  shocked  air,  followed  by  expanding 
detonation  products.  Because  of  the  large  compression  of 
air,  together  with  Interface  mixing,  the  bulk  of  the  blast 
pulse  Is  composed  of  detonation  products. 

Conductivity  of  the  gas  Is  principally  determined  by 
temperature  and  the  Ionization  potentials  of  the  gas  con- 
stituents. At  the  low  temperatures  of  expanded  explosive 
gases  (less  than  1500  K)  It  Is  necessary  to  "seed"  the  gas 
with  metals  such  as  potassium  or  cesium.  As  long  as  the 
seed  material  Is  only  a few  percent  the  temperature  of  the 
gas  Is  essentially  unchanged  by  seeding.  Tenqperature  Is 
therefore  determined  by  the  properties  of  the  expanding 
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explosive  gases,  and  conductivity  will  be  determined  by  the 
seed  material  present. 

We  are  not  aware  of  any  tenqperature  measurements  In 
expanded  explosive  products,  but  there  have  been  computer 
calculations  based  on  equilibrium  thermochemistry  assuming 
Isentroplc  expansion  from  the  detonation  state.  The  latest 
and  best  of  these  explosives  codes  Is  the  TIGER  code 
(Reference  1).  Published  calculations  were  available  only 
for  ThT  and  RDX;  the  results  are  shown  In  Figure  D2. 

Temperature  ranges  from  about  1100 *K  to  about  800 *K. 
Higher  temperatures  are  obtained  from  a relatively  low 
density  explosive,  presumably  because  of  the  Increased 
entropy  caused  by  compression  at  the  detonation  front. 
Temperatures  are  lowest  at  low  pressures,  as  Is  expected. 

Actual  temperatures  In  an  experiment  may  differ  sub- 
stantially from  these  curves  for  several  reasons.  The 
explosive  in  a DABS  test  Is  an  array  of  PETN  detonating 
cords.  The  gases  expand  radially  out  from  the  cords  at 
first,  then  stagnate  and  expand  as  a whole  at  a lower 
pressure  to  drive  the  blast  wave.  The  process  of  stagnation 
and  reexpeuislon  Increases  the  entropy  of  the  gas,  which 
Increases  Its  tenqperature  at  a given  pressure.  Furthermore 
there  Is  burning  of  the  gases  with  air  In  the  blast  tube, 
which  will  Increase  the  temperature. 

Conductivity  of  a gas  depends  on  the  number  of  electrons 
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present  and  the  average  collision  cross-section  of  the  gas 
molecules.  Electron  number  density  may  be  estimated  from 
equilibrium  chemical  thermodynamics  using  data  on  expanded 
explosive  products  together  with  a selected  seed  material. 

These  calculations  are  quite  lengthy  and  tedious;  It  turns 
out  that  certain  electronegative  compounds  (e.g.  CN’and  CO^) 
play  an  Important  role  In  suppressing  electrons  and  must  be 
Included  In  the  calculations. 

The  results  of  these  equilibrium  conductivity  calcu- 
lations are  shown  In  Figure  D3.  The  conductivity  Is 
strongly  temperature  dependent,  varying  many  orders  of 
magnitude  In  the  temperature  range  of  expanding  explosives. 

The  actual  quantity  of  Ionized  seed  material  In  the  gas  Is 
minute;  at  these  ten^eratures  the  material  Is  barely  melted, 
much  less  vaporized,  and  the  partial  pressure  of  the  vapor 
state  Is  very  small.  Consequently  It  takes  very  little  seed 
material  to  achieve  saturation.  However,  considerations  of 
mixing  and  evaporation  rates  suggest  putting  several  percent 
of  seed  material  Into  the  explosive.  The  difference  In  con- 
ductivity between  cesium  and  potassium  Is  about  a factor 
of  30,  but  the  cost  difference  Is  many  orders  of  magnitude. 

Nonequilibrium  effects  may  be  dominant  In  establishing 
the  electrical  conductivity,  because  the  electron  concentration 
Is  very  small  and  the  expansion  time  is  very  short  ( a few 
milliseconds).  For  example  nonequilibrium  chemistry  Is 
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largely  responsible  for  automobile  smog  production;  equilib- 
rium concentration  calculations  show  vanishingly  small 
concentrations  of  many  noxious  species. 

If  seed  material  Is  Imbedded  In  the  DABS  explosive  to 
assure  exposure  to  the  extreme  detonation  conditions,  the 
likely  result  of  nonequilibrium  effects  Is  to  Increase  con- 
ductivity by  delaying  electron  recombination.  If  the  seed 
material  Is  outside  the  explosive,  relying  on  late-tlme 
mixing  effects,  non-equl librium  effects  would  result  In  lower 
conductivity  due  to  delayed  evaporation  of  the  seed  material. 
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3 . EXPERIMENTAL  RESULTS 


The  objective  of  this  proprisim  (Artec  Project  122)  was 
to  determine  the  feasibility  of  the  magnetic  velocity  gage 
principle  for  AP\7L  DABS  tests.  Preliminary  calculations 
Indicated  that  the  critical  parameter  is  Internal  resistance, 
determined  by  gaseous  conductivity.  If  Internal  resistance 
Is  less  than  the  input  resistance  of  practical  Instrumen- 
tation amplifiers,  gage  output  will  accurately  m.easure 
velocity;  If  Internal  resistance  Is  greater,  gage  output 
will  be  reduced  by  a factor  depending  on  gas  conductivity. 

Phase  I consisted  of  five  explosive  experiments  con- 
ducted at  the  PMC  test  site  in  Hollister,  California.  All 
shots  used  essentially  the  same  configuration,  shown  In 
Figures  04  and  05.  The  explosive  was  contained  Ih*  a 
10-foot  long  12-lnch  dlam.eter  cardboard  tube  tamped  by 
several  feet  of  sand.  Four  5-foot  long  bundles  of  400-graln 
detonating  cord  were  suspended  approximately  6 Inches  apart 
In  the  first  half  of  the  tube.  The  explosive  driver  tube 
was  connected  to  a 20-foot  long  steel  pipe  to  achieve  proper 
pulse  amplitude  and  duration.  The  transition  pipe  was 
flanged  to  a 3*5-foot  long  stainless  steel  test  section, 
which  contained  magnetic  field  coil  and  electrodes.  An 
additional  10-foot  length  of  cardboard  tube  connected  to  the 
test  section  helped  to  maintain  pulse  duration. 
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Figure  D4.  Experimental  Setup  of  Shot  122-2 


Exp«rimentol  Setup,  Shot  122-2 


Test  Section 

Figure  D5.  Pre-Shot  Photographs  of  Shot  122-2 
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Shot  122-1  suooeBSfully  verified  the  meohanloal  perfor- 
mance and  poet-shot  Integrity  of  the  test  components.  Shot 
122-2  was  the  first  eleotrloal  test. 

Seeding  of  the  explosive  In  122-2  was  accomplished  by 
mixing  cesium  chloride  with  an  approximately  equal  volume  of 
PETN  explosive  powder  In  paper  tubes  taped  to  the  explosive 
charges.  An  amount  of  seed  material  equal  to  IS  by  weight 
of  the  total  explosive  charge  was  used. 

Electrodes  were  fabricated  as  an  8 x 20  array  of 
separate  elements  connected  electrically  In  parallel » which 
was  shown  by  auxiliary  laboratory  tests  to  work  as  well  as 
solid  electrodes  (as  shown  in  Figure  D6).  Each  element 
consisted  of  a U-shaped  piece  of  stainless  steel  wire  passed 
through  holes  In  the  test  section  sidewall.  Plastic  sleeves 
provided  eleotrloal  Insulation,  and  epoxy  provided  meohanloal 
rigidity.  The  Inside  of  the  stainless  steel  test  section 
was  coated  with  Cotronlos  ceramic  Insulating  material.  The 
electrode  elements  protruded  approximately  3/8-lnoh  and  were 
canted  at  an  angle  of  30  degrees  from  the  flow  to  Improve 
eleotrloal  contact  with  the  gas,  as  It  was  felt  that  a 
cold  boundary  layer  would  add  a high  resistance  to  the 
eleotrloal  circuit  In  the  gas.  The  electrode  constant 
relating  resistance  to  conductivity  is  .4  wT^* 

The  objective  of  122-2  was  to  measure  the  average  con- 
ductivity of  the  gas  by  means  of  the  circuit  shown  In  Figure 
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Figure  06.  Test  Section  and  Electrode  Detail 
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D7.  The  capacitor  discharge  tine  Is  long  (minutes)  compared 
to  gas  flow  duration  (milliseconds),  and  so  acts  as  a con- 
stant voltage  source.  The  gas  resistance  Is  determined  by 
measuring  the  source  voltage  and  the  voltage  drop  across  a 
known  resistance.  Average  conductivity  Is  obtained  from 
measured  resistance  using  the  known  electrode  constant. 

The  explosive  driver  produced  250-400  pslg  gas  pressures 
In  the  test  section,  representing  acceptable  DABS  test  con- 
ditions. Conductivity  measurement  was  unsuccessful  because 
of  a short-circuit  which  occurred  between  an  electrode 
element  and  the  test  section  tube  wall  during  the  test.  It 
was  also  noted  that  epoxy  seal  material  was  broken  from 
several  other  electrode  elements,  and  some  gas  escaped  through 
the  mounting  holes  possibly  providing  additional  electrical 
shunt  paths. 

Shot  122-3  was  essentially  a repeat  of  122-2  with  Improve- 
ments In  electrode  mounting.  All  the  electrode  elements  were 
removed  and  critical  mounting  holes  were  beveled  to  reduce 
the  chances  of  a short  circuit.  New  sleeves  were  Installed 
with  generous  amounts  of  epoxy.  The  Cotronlos  Insulator 
coating,  which  had  evidenced  some  cracking  and  chipping,  was 
removed  and  replaced  with  a high  temperature  insulating  paint 
on  both  the  Inner  and  outer  diameters  of  the  stainless  steel 
test  section,  along  with  additional  layers  of  high  teiiq;>erature 
duct  tape.  The  test  section  was  Isolated  electrically  from 
the  steel  transition  pipe. 
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C - Constant 
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Vm  s Voltogt  Drop  Across  R 

0-  * Conductivity 

Figure  07.  Gas  Conductivity  Measurement 
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By  direction  from  CERP  the  seed  material  was  changed 
from  oeaium  chloride  to  potassium  ohlorlde  on  the  basis  of 
economy  for  large  soale  DABS  tests.  Because  of  its  low  oost, 
5%  by  weight  of  potassium  ohlorlde  seed  was  used  Instead  of 
IS.  No  other  aspects  of  the  shot  were  changed  from  122-2. 

The  driver  operated  about  the  same  as  the  previous  shot, 
but  again  conductivity  data  was  lost  by  an  electrode  short 
olroult.  Post-shot  inspection  revealed  that  gas  again  had 
escaped  through  a few  of  the  electrode  element  mounting  holes 
causing  a low  impedance  shunt  to  the  stainless  steel  test 
section. 

Because  of  time  and  budget  constraints  the  next  shot 
was  configured  to  provide  both  velocity  gage  data  and  con- 
ductivity data.  Magnetic  field  was  provided  by  an  80-tum 
saddle  coil  with  an  average  B/I  of  3 >6  gauss/amp  over  the 
electrode  volume.  Four  automobile  batteries  in  series  pro- 
vided enough  current  to  generate  a field  of  180  gauss.  The 
electrode  elements  were  reworked,  and  a solid  meohanloal 
backing  to  prevent  gas  leakage  was  provided  by  high  density 
Hydrostone  casting  plaster.  The  Hydrostone  was  thoroughly 
dried  to  maintain  high  impedance  levels  (in  excess  of  20 
megohms).  Separate  1-lnoh  dianwter  electrodes  near  the 
upstream  end  of  the  test  section  were  included  to  obtain 
conductivity  data.  The  stainless  steel  test  section,  the 
pressure  gage,  and  the  conductivity  gage  electrodes  were  all 
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carefully  Isolated  electrically  to  minimise  spurious  ground 
paths . 

The  conductivity  gage  apparently  worked  well  In  122-4. 

Raw  data  from  the  tapedeok  Is  shown  In  Figure  D8,  and 
reduced  data  Is  shown  In  Figure  D9.  Conductivity  ranged 
from  10  to  over  100  mlcromhos/m,  with  peaks  at  the  shook 
front  and  what  appears  to  be  the  burning  Interface  between 
air  and  explosive  products.  This  conductivity  corresponds 
to  a velocity  gage  resistance  from  .4  to  4 megohm,  which  Is 
high  but  may  be  acceptable  with  high  Input  Impedance 
amplifiers. 

Velocity  gage  data  In  this  shot  showed  large  oscillations. 
A major  factor  contributing  to  this  effect  was  determined 
after  tl  a shot  to  be  electrical  crosstalk  occurring  In  the 
gas  between  the  conductivity  gage  electrodes  and  the  velocity 
gage  electrodes.  The  conductivity  gage  had  nearly  300  volts 
across  It,  while  the  velocity  gage  had  about  5 volts.  Trickle 
currents  through  the  gas  produced  several  volts  at  the 
velocity  gage  electrodes  from  the  conductivity  gage  elec- 
trodes; this  effect  was  confirmed  by  lab  tests  using  copper 
sulphate  solution. 

Considerable  effort  was  put  into  the  electrical  design 
of  shot  122-5  to  Isolate  the  high  Impedance  velocity  gage 
electrodes  from  other  circuitry,  as  shown  In  Figure  DIO. 

A PCS  Model  402M41  voltage  follower  amplifier  was  used. 
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Pressure  Gage 
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Figure  08.  122>4  Raw  Data 


Tim*  (m«) 


Figure  D9.  122-4  Reduced  Oita 
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which  has  an  Input  Impedance  of  10^^  ohms.  The  conductivity 
measurement  was  eliminated*  and  In  addition  the  voltage 
level  of  the  velocity  gage  was  halved  to  eliminate  clipping 
by  using  only  two  batteries  to  power  the  field  cull  Instead 
of  four.  All  other  aspects  of  the  shot  were  Identical  to 
122-4. 

Raw  data  from  122-5  Is  shown  In  Figures  Dll  and  D12* 
and  reduced  data  Is  shown  In  Figure  D13.  In  general  the 
shot  performed  very  well  except  for  anomalous  pulses  In  the 
velocity  gage  data  at  shock  arrival  and  at  pulse  termination. 
Measured  velocities  appear  reasonable*  although  they  have 
not  been  compared  to  a detailed  theoretical  calculation  of 
the  experimental  setup.  The  velocity  data  appears  smoother 
than  the  pressure  data*  as  though  It  were  filtered  by  RC 
effects  In  the  generator. 

Post  shot  analysis  has  Indicated  that  the  likely  source 
of  the  anomalous  voltage  pulse  Is  a static  electricity 
effect  associated  with  the  shock  tube  boundary  layer.  Orosse 
(Reference  2)  has  Investigated  static  charge  phenomena  In 
shook  tubes*  and  found  coupling  effects  similar  to  what  was 
observed  In  122-5.  Grosse  concluded  that  the  mechanism  was 
a charged  boundary  layer  produced  by  shearing  and  compression 
effects  Inducing  a voltage  on  pickup  electrodes  by  capacitive 
coupling.  This  Interpretation  Is  supported  by  observations 
of  voltage  on  the  stainless  steel  tube*  which  Is  also  coupled 


Top:  Stainless  Steel  Voltage 
10  V/cm  5 ms/cm 


Bottom:  Velocity  Gage 
3.18  km/sec/cm 
5ms/cm 


Top:  Stainless  Steel  Voltage 
5 V/cm  2 ms/cm 


Bottom:  Pressure  Gage 

I06psi/cm  2 ms/cm 


Figure  Dll.  122-5  Raw  Data  From  Scope  Traces 
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Figure  012.  122-5  Raw  Data  from  Tape  Deck 


capacltlvely  to  the  surface  charge  but  with  capacitance 
history  different  from  the  electrodes. 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 
4.1  Conductivity  Data 

Shot  122-4  successfully  obtained  conductltlty  data  on 
the  test  gas.  It  ranged  from  10  to  100  mlcroirhos/m,  which 
Is  several  orders  of  magnitude  better  than  preliminary 
equilibrium  calculations,  but  consistent  with  expectations 
based  on  nonequilibrium  phenomena.  The  measured  conduc- 
tivities correspond  to  gas  resistances  In  the  velocity  gage 
ranging  from  .4  to  4 megohm.  Use  of  cesium  Instead  of 
potassium  would  theoretically  lower  this  resistance  by  a 
factor  of  about  50,  but  It  Is  a costly  material  for  large 
tests. 

V'lth  gas  resistance  In  the  megohm  range  it  Is  necessary 
to  use  a high  Input  Impedance  amplifier,  and  a commercial 
FET  an^llfler  with  Input  Impedance  of  10^^  ohms  was  used 
In  122-5. 

At  these  high  Impedance  levels  shunt  capacitance 
becomes  a problem.  For  example  consider  the  circuit: 


Scop* 
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This  circuit  acts  as  an  RC  filter  on  the  velocity  gage 
voltage.  To  achieve  reasonable  fidelity  in  a typical  DADS 
blast  pulse,  RC  < 1 msec.  For  a gas  resistance  of  one 
megohm,  this  condition  limits  allowable  shunt  capacitance 
to  1 nanofarad.  A higher  gas  resistance  would  require  a 
lower  shunt  capacitance. 

Sources  of  shunt  capacitance  are  connecting  cable, 
capacitance  of  electrodes  to  electromagnet  coll,  and 
capacitance  of  electrodes  to  metal  structures.  In  the 
particular  case  of  the  Artec  tests  at  FMC  test  site  the 
shunt  capacitance  was  3*0  nf.  This  was  unusually  high 
because  of  a stainless  steel  channel.  In  a DABS  test  where 
the  entire  gage  Is  Immersed  In  the  blast  It  can  be  made  of 
fibreglass,  resulting  In  lower  shunt  capacitance. 

Velocity  gage  data  was  successfully  obtained  In  Shot 
122-5.  It  was  smoother  than  the  corresponding  pressure  gage 
record,  so  some  degree  of  RC  filtering  was  probably  present. 
Velocity  magnitudes  appear  reasonable,  but  neither  computer 
calculations  of  the  experimental  setup  nor  alternate  velocity 
measurements  were  available  as  a cheek. 

A static  electricity  noise  problem  was  Identified  In 
122-5.  A sharp  negative-going  pulse  occurred  at  shock 
arrival,  and  a slower  oscillation  occurred  toward  the  end 
of  the  blast  wave. 

Static  eleetl*lel'^v  noise  was  considered  as  a potential 
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problem  for  high  log>edance  velocity  gage  meaeuremente  early 
In  the  program,  and  this  effect  seemed  to  be  borne  out  by 
an  experiment  at  AFWL  (Figure  D14  ) . In  that  experiment  30 
volts  of  electrostatic  noise  was  measured  between  two  metal 
plates  resting  on  a DABS  shook  tube  floor.  The  explosive 
was  not  seeded. 

Further  consideration  of  this  effect  indicated  that 
seeding,  if  at  all  effective^  should  eliminate  the  px^blem. 
Charge  concentrations  cannot  persist  in  a conducting  gas; 
they  equilibrate  on  a time  scale  of  the  order 


-•8 

Even  low  conductivities,  on  the  order  of  10’”  mho/m,  lead 
to  an  equilibration  time  of  less  than  a millisecond.  If 
the  gas  is  sufficiently  conductive  for  the  magnetic  velocity 
gage  to  work,  static  charges  cannot  be  carried  along  by  the 
flowing  gas,  and  this  type  of  noise  should  not  be  a problem. 

The  electrostatic  noise  measured  in  122-5  is  apparently 
of  a different  nature  from  the  noise  measured  by  AFWL. 

Noise  in  the  AFWL  test  (Figure  D14)  was  rapidly  oscillating 
with  oscillations  persisting  throughout  the  pulse  duration. 
Noise  in  Artec  Shot  122-5  consisted  primarily  of  a single 
oscillation  at  the  start  of  the  pulse,  and  a longer  single 
oscillation  toward  the  end  of  the  pulse. 


(04) 
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Figure  D14.  AFWL  Noise  Test  Raw  Data 
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This  type  of  noise  appears  to  be  related  to  surface 
electrostatic  phenomena  previously  observed  In  shook  tubes. 
Qrosse  (Reference  2)  has  made  extensive  observations  of 
these  electrostatic  effects,  and  concludes  that  they  can  be 
explained  by  a static  charge  In  the  boundary  layer.  The 
magnitude  of  this  charge  layer,  about  10  ^ coulomb/m  , Is 
sufficient  to  Induce  voltages  In  sensing  electrodes  by 
capacltatlve  coupling.  The  charge  Is  produced  by  compression 
or  shearing  effects  at  the  boundary.  Since  the  bouhdary 
layer  Is  col6,  the  conductivity  Is  low  and  the  charge 
persists  a long  time. 

Assuming  a surface  area  of  1 sq.  meter,  which  Is  about 
right  for  the  Artec  experiments,  the  10'^  coulomb  charge 
would  produce  a 1 volt  signal  across  a 1 nanofarad  capacitor. 
Thus  the  observed  voltage  Is  of  correct  order  of  magnitude 
to  be  explained  by  the  boundary  layer  charge  effect. 

Electrostatic  boundary  layer  effects  and  their  effects 
on  high  Impedance  measuring  circuits  are  poorly  understood 
at  present,  but  there  Is  no  reason  to  believe  that  they 
cannot  be  circumvented  by  an  appropriate  circuit  design. 

In  an  actual  DABS  test  the  magnetic  velocity  gage  would 
be  located  In  the  core  flow,  well  away  from  the  shook  tube 
boundary  layer.  In  that  case  the  only  boundary  layer  Is  the 
one  formed  over  the  electrodes  themselves,  and  any  resulting 
effects  are  likely  to  be  small.  Even  In  the  Artec  test 
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setup  the  electrodes  could  be  moved  further  away  from  the 
wall  to  minimize  capacitive  coupling  to  the  boundary  layer 
charge . 

4.2  Recommendations 

For  the  Immediate  DABS  test  program  our  recommendation 
Is  that  the  MHD  velocity  gage  not  be  fielded  In  a Phase  II 
test  effort.  It  Is  possible  to  field  a high  Impedance  gage, 
but  the  associated  practical  problems  of  static  noise,  stray 
capacitance,  and  spurious  pickup  are  not  likely  to  be  solved 
In  time  for  the  test  program.  It  Is  felt  that  with  further 
effort  the  high  Impedance  difficulties  could  be  overcome. 

For  tests  conducted  with  a higher  temperature  gas  (e.g. 
a tunnel  test  using  a high  energy  source  or  a high  pressure 
DABS)  the  magnetic  velocity  gage  Is  likely  to  prove  more 
effective.  Conductivity  Is  highly  sensitive  to  temperature, 
and  any  seeded  flow  over  about  1200  K would  have  Internal 
resistance  well  below  a megohm,  thereby  eliminating  high 
Impedance  problems. 
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INTRODUCTION 


The  Dynamic  Air  Blast  Simulator  (DABS)  is  an  explosively -driven 
shock  tube  being  developed  to  simulate  nuclear  blast  wave  i^enomena  on  a 
variety  of  structures.  As  a simulation  facility  it  should  reasonably  reproduce 
the  pressure  history  of  a nuclear  blast  wave.  Several  small  scale,  exploratory 
DABS  have  been  built  and  fired.  These  tests  have  generally  demonstrated 
reasonable  pressure  levels  and  pressure  time  histories.  However,  several 
questions  remain:  (i)  at  least  two  of  the  three  tests  to  date  have  shown 
"spurious"  double  spiked  pressure  traces;  the  origin  of  this  phenomenon  is 
unclear,  (ii)  Since  the  utility  of  this  facility  is  structure  loading,  a direct 
measure  of  the  dynamic  forces  is  important,  (iii)  Scaling  laws  need  to  be 
verified. 

To  provide  further  diagnostic  concepts  as  an  aid  in  understanding  the 
operating  performance  and  scalability  of  die  DABS  facility,  a meeting  was 
held  at  CERF,  May  18  and  19,  1976.  At  this  working  group,  the  DABS  test 
data,  diagnostics,  and  technical  issues  were  discussed.  Recommendations 
were  made  for  additional  diagnostics  for  the  next  test  scheduled  for  June  17th. 

This  letter  report  documents  this  consultant's  ideas  concerning  radia- 
tion measurements  that  might  be  attempted  on  the  next  and  future  DABS  tests. 

RADIATIVE  DIAGNOSTICS 


A.  Objectives 

Measurement  of  the  emitted  or  absorbed  radiation  is  an  important  and 
highly  developed  diagnostic  for  high  temperature  gases.  The  intensity, 
temporal  history,  and  spectral  content  of  radiation  can,  in  principle,  yield 
information  about  gas  density,  temperature,  state  of  excitation,  and  species 
identification  to  mention  a few  important  variables.  With  the  recent  availability 
of  lasers  covering  a wide  range  of  wavelengths,  many  sophisticated  optical 
measuring  techniques  have  been  developed. 
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However,  the  emphasia  on  this  study  has  been  to  consider  simple 
optical  systems  which  have  the  capability  of  providing  overall,  gross 
measurements.  The  rationale  for  the  above  emphasis  is  several  fold: 

(i)  The  very  hostile  environment  of  DABS  precludes  sophisticated  experi- 
ments and  expensive  equipment,  (ii)  The  time  scale  for  the  next  test 
(June  17th)  precluded  any  instrumentation  development,  (iii)  At  present  it 
is  more  cost  effective  to  proceed  with  simple  experiments  which  are  easy 
to  interpret.  More  costly  diagnostics  might  be  developed  later  if  required. 

The  recommendations  of  the  working  group  in  the  optical  area  were 
to  attempt  to  measure  infrared  emission  at  several  selected  wavelengths 
using  simple  IR  radiometers  and  transmission  in  the  visible  using  both  a 
pulsed  ruby  laser  (X  = 6943  and  a broadband  source.  The  following 
sections  briefly  describe  some  of  the  technical  issues  involved  in  attempting 
to  field  and  interpret  such  measurements. 

B.  Environment 

The  unique  environmental  conditions  of  DABS  make  it  difficult  to  design 
optical  diagnostics  which  have  unequivocal  interpretation.  For  example,  the 
construction  of  the  shock  tube  consists  of  a trench  in  the  desert  with  concrete 
walls,  dirt  floor,  and  a wood  roof.  This  construction  creates  unusual 
problems  with  radiation  diagnostics.  It  does  not  appear  possible  to  completely 
characterize  the  driven  gas  conditions  as  to  temperature,  density,  and  species 
concentrations.  There  is  evidence  for  a substantial  carbon  (soot)  content  in 
the  driver  gas.  This  may  result  from  the  fact  that  the  PETN  explosive  con- 
tains up  to  20%  of  "inert"  plastic  binder  material.  There  is  data  which 
indicate  that  the  driver  interface  is  a highly  turbulent,  combusting,  and 
radiating  gas.  Finally,  the  high  densities  of  operation  (^  atmospheric)  and 
long  path  lengths  (6-8  ft)  suggest  that  optical  opacity  may  be  serious.  There 
is  photographic  evidence  of  an  opaque  region  that  accompanies  the  blast  wave. 

From  previous  DABS  test  data^,  it  can  be  seen  that  the  average  velocity 
of  the  shock  is  ~6500  ft/sec  (2  mm/^sec).  Based  on  this  shock  velocity,  the 
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poat  shock  conditions  for  equilibrium  air  are: 

T2  ~ 2000°K 

Pj/p,  ~ 

P^/Pj  ~ 40 

where  T = temperature,  p = density,  and  P = pressure.  Subscripts  refer 

to  (1)  initial  unshocked,  and  (2)  post-shock  conditions.  A graph  of  the 

equilibrium  composition  of  air  for  2000  K is  shown  in  Fig. El.  For  the 

shocked  gas  density  of  DABS  (log  p/p^  ~ 0.7),  it  can  be  seen  that  the  gas 

would  contain  ~ 1%  NO,  0.04%  CO.  (ambient)  and  a small  amount  of 
-2  ^ 

NO^  (~  10  %).  These  are  the  only  infrared  active  species.  In  addition,  the 

electron  concentration  is  obviously  very  low  and  probably  dominated  by  easily 
ionized  impurities  such  as  Na.  The  assumption  of  equilibrium  conditions  is 
justified  because  of  the  high  gas  densities  involved.  For  example,  an 
estimate  of  the  time  to  form  NO  under  these  conditions  is  ~ 1 ^sec.  ^ 

Unfortunately,  the  explosive  driver  conditions  are  far  less  easy  to 

characterize.  Martens^  has  used  a hydrocode  to  calculate  blast  wave  strength 

for  explosive  loadings  and  geometries  comparable  to  DABS  tests.  However, 

thi4  code  does  not  determine  explosive  product  species  concentrations  or 

5 

temperature.  Cook  gives  "typical"  detonation  concentrations  for  PETN, 
but  the  details  and  validity  of  these  results  is  questionable.  Nevertheless, 
Cook's  results  would  imply  that  N^,  CO2,  H2O  and  CO  are  all  probably 
present  in  roughly  equal  concentrations.  The  temperature  of  the  detonation 
products  initially  is  about  2000  - 3000°K  and  decreases  as  the  gas  expands. 
There  is  evidence  of  a highly  luminous  contact  surface  which  propagates  down 
the  tube  just  behind  the  shock  front.  This  luminous  front  is  caused  by  further 
combustion  of  the  detonation  products  as  they  mix  with  the  hot  air.  A specifica- 
tion of  the  conditions  within  this  burning  interface  is  not  practical. 


180 


Mole  Percent 


Another  feature  of  this  facility  is  the  production  of  the  substantial 
amounts  of  carbon  (soot)  by  the  detonation/combustion  process.  The  origin 
of  this  free  carbon  is  not  clear,  but  one  reasonable  source  is  the  plastic 
binder  of  the  PETN  primer  cord  which  burns  and  vaporizes  in  the  detonation. 
For  example,  for  600  lbs  of  PETN  there  will  be  about  85  lbs  of  plastic.  If 

3 

10%  of  this  plastic  material  is  converted  into  carbon,  ~ 4 x 10  gms  of  C are 

produced.  If  this  carbon  is  uniformly  distributed  throughout  the  entire  tube, 

-5  3 

the  mass  loading  of  carbon  is  about  3x10  gms/cm  . The  particle  size 
distribution  is  not  known.  However,  soot  formed  by  hydrocarbon  combustion 
typically  shows  sub-micron  particle  size.  With  this  assumption  and  the  large 
mass  loading  of  C we  can  calculate  an  absorption  length,  L,  in  the  particle 
laden  gas  in  the  volume  absorption  limit. 


where  M is  the  nuiss  loading  of  particles  of  density,  p and  a i*  the  volume 
absorption  coefficient  for  C which  is  shown  in  Fig.E2,  For  the  mass  loading 
considered  L ~ 10  cm  and  the  gas  is  generally  optically  thick  in  the  infrared. 

The  situation  would  even  be  worse  (smaller  optical  depth)  in  the  visible.^ 

This  calculation  is  very  approximate,  but  it  would  appear  that  the  particle 
laden  detonation  products  would  make  it  difficult  if  not  impossible  to  perform 
absorption  or  molecular  scattering  (i.e.  Raman)  measurements.  However, 
if  the  gas  is  optically  black,  then  it  can  be  treated  as  a blackbody  (see  below). 

C.  Radiation  Techniques 

A number  of  infrared  optical  techniques  were  considered  to  determine 
the  temperature  of  the  shock  heated  gas  or  detonation  products  as  they  flow 
down  the  tube:  emission  from  neutral  Bremsstrahlung,  emission  from  CO2, 

NO,  etc.,  absorption  of  CO^  laser  lines,  and  two  color  measurement  of  blackbody 
radiation. 
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NORMALIZED  RAYLEIGH  ABSORPTION  CROSS  SECTION 


4 8 12  16  20  24  28 


WAVELENGTH.  Mm 

Figure  E2.  Volume  absorption  coefficient  vs,  wavelength  for  soot  and 
fused  silica,  T = 300®K, 


183 


1.  Neutral  Bremsatrahlung  is  continuum  radiation  produced  by 

inelastic  scattering  of  electrons  from  neutral  atoms  or  molecules.  The 

appropriate  cross  sections  have  been  experimentally  determined  for  the 

7 

major  constituents  of  air  N2  and  and  the  technique  has  been  used  for 

temperature  measurements.  However,  this  method  is  clearly  inappropriate 
for  this  application  due  to  the  very  low  and  undefined  electron  concentration. 

Z.  In  principle  the  measurement  of  the  absolute  emission  from 
CO^  (or  other  infrared  active  molecules)  behind  the  incident  shock  could  be 
used  to  determine  gas  temperature.  The  shock  heated  air  is  at  equilibrium 
with  the  CO2  at  its  ambient  concentration.  A simple  IR  radiometer  could 
be  constructed  with  appropriate  filter  and  calibrated  in  a small,  pressure 
driven  shock  tube.  An  example  of  this  type  of  measurement  is  shown  in  Fig.E3 
taken  from  one  of  the  author's  own  papers.  This  figure  shows  the  emission 
from  the  first  overtone  band  of  CO  at  2.  3p.  The  expected  temperature 
sensitivity  for  the  4.  3 band  of  CO2  would  be  similar  to  the  dashed  curve  on 
Fig.E3. 

The  difficulty  of  this  technique  is  that  the  interpretation  is  simple  only 
if  the  gas  is  optically  thin.  The  optical  opacity,  )i,  for  infrared  molecular 
radiation  may  be  defined  in  terms  of  the  product  of  path  length,  4,  and  density 
of  species,  P,  in  atmospheres 

yx  = P I ~ 0.  5 atm -cm  for  CO2  in  DABS 

This  is  a value  of  |a  for  which  the  fundamental  infrared  bands  of  species  such 
as  CO2*  CO,  NO  and  H2O  are  generally  black.  For  example,  Fig.E4  shows 
the  emissivity  (proportional  to  intensity)  of  the  4.  3^  fundamental  band  of  CO. 
for  various  values  of  Pi  for  1500  K.  In  the  region  of  the  band  head  (~2300  cm  ). 
it  can  be  seen  that  the  emissivity  curves  do  not  scale  with  Pi  for  Pi  > 0. 1. 

The  curves  do  scale  in  the  wing  of  the  band,  X < 2200  cm"^.  However,  the 
intensity  is  much  weaker  int\e  wings  and  obviously  more  difficult  to  detect 
and  me.asure  quantitatively.  In  principle,  optical  calibration  of  this  technique 
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Figure  E3.  Data  from  incident  shock  calibration  on  the  CO  first  over- 
tone emission.  The  parameter  B = S/A|y^n^Q  , defined  in 
the  text,  is  proportional  to  radiated  power  per  particle  per 
steradian.  The  dashed  curve  shows  the  temperature  depend- 
ence of  the  fundamental  band  emission. 
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Figure  E4.  CO|  4»3-ii  bend  emieelvlty  versus  wsvenumber  for  pure 
Doppler  line  shape  and  weak  line  approximation  for  T ■ 1500K. 
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could  be  carried  out  under  conditions  of  identical  Pi  even  for  Pi  2 0. 1. 

In  the  small  shock  tube  the  pressure  of  CO2  would  have  to  be  increased 
to  account  for  the  smaller  i.  However,  interpretation  becomes  less 
certain,  and,  as  the  gas  becomes  more  optically  thick,  the  technique 
samples  boundary  layer  gas  which  is  not  of  interest. 

3.  It  is  well  known  that  the  relative  population  of  rotational  states 
of  a specific  vibrational  level  of  a molecule  can  be  used  to  determine  the 
rotational  temperature  of  a gas.  Since,  under  DABS  conditions,  the  various 
internal  modes  of  the  gas  are  in  equilibrium,  a measurement  of  rotational 
temperature  would  yield  gas  temperature.  The  population  of  (he  rotational 
levels  of  CO2  can  be  measured  in  absorption  by  using  a probe  CO2  laser. 

This  measurement  has  several  inherent  advantages.  First,  the  CO2  laser 
radiation  is  not  absorbed  to  any  appreciable  extent  by  the  room  temperature 
CO2  in  the  atmosphere  (<  1%/km  for  most  CO2  laser  lines).  Hence,  the 
probe  laser  could  be  placed  a reasonable  distance  from  the  facility.  Second, 
behind  the  incident  shock  the  CO2  is  heated  to  ZOOO^K  and  this  temperature 
rise  provides  enough  increase  in  CO2  (100)  population  to  create  significant 
and  measurable  absorption  over  the  200  cm  DABS  path  length.  Finally,  deter* 
mination  of  only  two  rotational  state  populations  is  necessary.  It  is  possible 
to  operate  a CO2  laser  on  two  individual  transitions  simultaneously,  hence 
the  same  source  and  collection  optics  could  be  used  for  this  measurement. 

However,  the  technique  is  not  recommended  at  this  time  because  the 
absorption  and  scattering  of  the  laser  radiation  by  the  soot  particle  is  a 
possible  serious  degradation.  Furthermore,  in  the  driver  gases,  large 
turbulent  density  gradients  may  exist  and  prevent  propagation  of  a well 
collimated  beam  across  the  tube.  This  technique  could  be  implemented  at  a 
later  date  after  some  of  the  questions  concerning  particles  and  turbulence 
have  been  resolved. 

4.  Since  the  gas  is  probably  optically  black  because  of  a high  mass 
loading  of  carbon  particles,  one  very  simple  optical  technique  to  measure 
temperature  is  to  use  a two  color  radiometer.  Figure  ESshows  the  spectral 
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SPECTRAL  RADIANCE  (W  sr 


FI9U  C E5*  SpFTlral  radianrc  of  a blarkbody,  N,^,  al  ibr  Irmprraiurc  in  ’K  thiiwn  on  rarli  rune; 
beyond  100  m ibe  cur»r»  ronliniir  lo  lie  linear  «|dnlli-d  on  l••((Jrilllllli<'  nale*!.  The  di-itionjl  line, 
inlrruerlinK  llie  curtri  at  their  maaitna,  ahuwa  WienV  di*>|«laeenH'nt  law.  tSulidi«i*iims  of  the 
ordinate  icaie  arc  al  2 and  S.)  Figure  by  P.  K.  Caal. 
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radiance  of  a blackbody,  as  a function  of  wavelength  for  various 
temperatures.  It  can  be  seen  that  for  the  conditions  existing  in  DABS, 
i.  e. , T ~ 2000°K,  the  peak  of  the  blackbody  function  is  1-2 

Many  techniques  have  been  developed  to  determine  temperature  from 
measurements  of  blackbody  radiance.  For  the  present  application,  two 
points  must  be  considered:  (i)  the  radiance  from  the  detonation  products 
will  be  composed  of  contributions  from  both  molecular  bands  and  continuum 
from  the  particles,  (ii)  The  emissivity  of  carbon  particles  is  both  a 
function  of  size  and  wavelength. 

g 

The  first  point  is  the  more  important.  Ferriso  et  al  have  developed 
an  infrared  band -ratio  technique  for  temperature  determination  under  the 
condition  of  radiation  from  both  molecular  bands  and  particle  continua.  This 
technique  was  developed  for  measurements  in  rocket  exhausts  for  which  the 
gas  conditions  are  similar  to  those  being  considered  in  DABS.  Ferriso  et  al 
suggest  that  the  most  appropriate  wavelengths  for  observation  of  the  continuum 
radiation  are  at  2.3  and  4.0  ^.  FigureE6is  an  example  of  the  expected  emission 
from  the  high  temperature  molecular  radiators  that  can  be  expected  in  the 
DAB  environment.  Note  that  this  figure  is  for  3Q00°K.  As  is  evident,  the 
wavelength  region  of  2.  3 and  4.  0|a  are  natural  valleys  between  the  emissivity 
peaks  of  bands  from  CO2,  CO,  and  H2O. 

The  emissivity  of  carbon  particles  is  a somewhat  controversial  issue. 

9 

Existing  measurements  do  not  all  agree,  but  clearly  indicate  particle  size 
and  composition  dependence.  However,  it  appears  that  the  emissivity  of  small 
(sub-micron)  carbon  particles  is  probably  ~0.  9 and  reasonably  constant  in  the 
infrared  spectral  region.  Therefore,  in  the  present  application,  the  existing 
uncertainty  in  carbon  particle  emissivity  does  not  introduce  a serious  error. 

CONCLUSIONS 

The  following  recommendations  have  been  developed  in  this  consulting 

study: 

(i)  The  environmental  conditions  existing  in  the  DABS  facility  pre- 
clude developing  any  sophisticated  radiation  diagnostics  at  this  time.  In 
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Figure  E6.  Example  of  Expected  Emission  from  High  Temperature 
Molecular  Radiators  In  DAB  Environment  (1  of  2) 
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Figure  E6.  Example  of  Expected  Emission  from  High  Temperature 
Molecular  Radiators  In  DAB  Environment  (2  of  2) 
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particular,  the  possible  high  mass  loading  of  particles  and  the  large 
optical  depth  of  molecular  radiators  indicate  that  the  gas  (at  least  the 
detonation  products)  are  probably  opaque  and  volume  sampling  may  be 
impossible.  In  addition,  high  turbulence  levels  and  interface  combustion 
indicate  that  the  gas  is  inhomogeneous,  and  it  may  be  difficult  to  probe  the 
medium  with  any  external  source  (laser  absorption,  Raman  scattering,  etc. ). 

(ii)  The  environmental  conditions  described  above  suggest  that  simple 
radiation  diagnostics  should  be  developed  to  (a)  determine  qualitative  under- 
standing of  the  DABS  gas  conditions  and  (b)  used  in  conjunction  with  existing 
pressure  probes  to  provide  some  further  information  about  shock  tube  flow 
parameters.  In  particular,  it  is  important  to  ascertain  whether  the  gas  is 
optically  black  and/or  whether  a probe  beam  can  be  propagated  across  the 
tube.  Simple  transmission  measurements  using  a probe  laser  (Ruby  at  6943 
as  well  as  a quartz  lamp  would  provide  some  initial  information  in  this  regard. 
These  measurements  should  have  sufficient  spatial  and  time  resolutions  that 
the  signals  could  be  correlated  with  existing  pressure  probes.  This  correla. 
tion  might  yield  information  about  the  separation  of  the  shock  and  driver  gases 
and  help  understand  the  "two  spike"  pressure  records. 

(iii)  In  addition,  it  is  reccmmended  that  simple  infrared  radiometers 
be  developed  and  used  on  future  DABS  tests.  If  the  gas  is  optically  black,  it 
should  be  possible  to  determine  gas  temperature  by  IR  emission  measure - 
ments.  Either  simple  overall  band  measurements  («  caT  ) can  be  used  or  the 

g 

more  sophisticated  band  ratio  method  can  be  employed.  In  any  event,  it 
should  be  possible  to  provide  an  estimate  of  gas  temperature  which  can  be 
correlated  with  otiier  variables  of  shock  performance. 
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APPENDIX  F 

CERF  INTRUNENTATION  RECOMMENDATIONS 
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MEMORANDUM 


Due  to  the  apparent  nature  of  the  flow  field,  measurements  that  can 
reliably  be  accomplished  while  still  yielding  significant  Information  with 
regard  to  the  nature  of  the  flow  filed  are  as  follows: 

1.  Shock  Wave  Velocity,  which  Is  currently  being  measured  very  success- 
fully. 

2.  Force  Measurements. 

a.  Static  pressure,  being  measured  at  several  locations. 

b.  Total  pressure,  being  measured  at  several  locations. 

c.  Drag,  not  being  measured. 

3.  Radiative  heat  transfer,  not  being  measured. 

Local  Mach  number  can  be  successfully  determined  by  a simultaneous  measure- 
ment of  static  pressure,  total  pressure  and  drag.  The  drag  measurement  Is 
necessary  In  order  to  determine  the  effective  ratio  of  specific  heats.  If  this 
Is  assumed  known,  a drag  measurement  Is  necessary. 

The  following  governing  equations  relate  the  measured  quantities  to  Mach 
number,  M,  and  the  effective  ratio  of  specific  heats,  y. 


Total  pressure/static  pressure  ratio. 


Po2 


(1 


+ id  M^)^^ 


For 


M < 1 


Y 

■^2  + ^ ( 2y 

2 ^ 


-1 

Y.1)Y^  M>1 

"+T 


Drag  D_ 
PS 


(M) 


(Fla) 

(Fib) 


(F2) 
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This  drag  coefficient  Is  a single-valued  function  of  Mach  number  and  Is 
virtually  Independent  of  the  wake  and.  at  the  expected  conditions  generated 
by  the  MBS.  Independent  of  Reynolds  number.  This  drag  coefficient*  as 
a function  of  Nach  number  Is  shown  In  FloureFI  by  the  open  symbols.  The  closed 
symbols  represent  the  product  of  M^Cq  (H).  These  data  are  then  least  squared 
to  the  following  function: 

AM^ 

A and  b were  found  to  be: 

A > 1.048 
b - 2.487 

with  a goodness  of  fit  0.99 

0,  Po2>  and  P are  measured  and  S Is  the  cross-sectional  area  of  the  cylinder. 
Once  the  left-hand  side  of  equations  FI  and  F2  are  known,  a solution  for  y and 
N can  be  found  by  any  one  of  a number  of  numerical  techniques.  The  bisection 
method  or  a similar  technique  Is  recommended  rather  than  Newton -Raphson  due 
to  the  problems  encountered  near  Mach  one  In  the  sensitivity  of  the  result.  If 
Y Is  assumed  to  be  known  then  the  drag  measurement  Is  unnecessary.  It  should 
be  noted,  however,  the  relative  sensitivity  of  the  technique  for  an  error  In  y. 

At  M-1 . 


L A 0 . ^ in  (^)  ^ 


For  a 5X  error  In  y the  following  errors  result  In  Pp/p 


Y 

**o/p 

1.1 

31 X 

1.2 

18.5X 

1.3 

14. 5X 

1.4 

12. 5X 

It  Is  therefore  recommended  that  a measurement  of  drag  be  made  In  conjunc 
tion  with  pressure  to  determine  the  effective  y.  In  addition,  at  high  mach 
numbers,  M > 0.8.  Co  (M)  Is  approximately  constant  and  the  dynamic  pressure 
can  easily  be  determined  since: 


D 

~Cq5”  = 


9 


*Ref:  Hoernen,  Fluid  Dynamic  Study  Drag.  Figure  A,  Transonic  Flow  Chapter. 
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O DRAG  COEFFICIENT 
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Figure  FI.  Blunt  Cylinder  Drag  as  a Function  of  Mach  Number 


Also,  I may  be  approximated  between  M equalling  0.8  to  2.0.,  since  Cq 
Is  a regular  function  and  relatively  slowly  varying.  ^ 

Another  Interesting  observation  Is  the  effect  of  compressibility  on 
the  Beroulll  result. 


^ . (t  ♦ ■ 1 ♦ y t Jm’ ♦ +11^  N®  + 0(H®)  (F4) 

Po  - 1+q  1 m2  + ^ + 0 (m6)  (F5) 

r f r TT 


The  bracketed  term  to  order  M^  and  has  the  following  values; 


Y 

M 

1.2 

1.4 

0.8 

1.17 

1.17 

1 

1.28 

1.28 

1.2 

1.43 

1.41 

If  the  M > 1 value  Is  used: 


' 1 ♦ J (1.28)  (F6) 

then  q ' (F7) 

Which  gives  no  more  than  a 15%  error  In  g for  this  Mach  number  range.  Thus, 
a quick  approximation  for  g can  be  made  by  using  a reasonable  value  of  the  bracketed 
coefficient  at  a given  point  time  In  time. 

Radiative  Heat  Transfer 


Much  Information  can  be  gathered  through  a wide-wall  radiative  heat-transfer 
measurement.  H^O  and  CO2  are  expected  combustion  products,  and  therefore, 
monitoring  their  vibration-rotation  spectra  In  the  near-litfrared  should  yield 
Information  with  regard  to  the  nature  of  the  shock-heated  region  and  the  relative 
location  of  the  deflagration  and/or  Indicate  that  the  gas  Is  optically  thick  at 
this  wavelength,  the  resulting  radiation  can  be  directly  related  to  the  gas 
temperature  through  the  black  body  function. 
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for  A » 1 
gr  ■ (VF) 


Au 


(F8) 


VF  Is  a view  factor  and  Au  Is  the  bandwidth  of  the  filter.  Note  that  this 
approach  makes  the  measurement  Independent  of  the  gas  properties  and  relative 
concentrations. 

Even  If  this  Is  not  the  case,  this  experimental  measurement  will  yield 
Information  with  regard  to  the  relative  locations  of  the  shock  front  and 
flame  and/or  combustion  products  front. 

It  Is  suggested  that  two  measurements  be  made:  One  over  a narrow  region 
In  the  Infrared  and  one  broad  band  over  the  wavelength  region. 

Narrow  Region  2.5  - 2.8  y 

Broad  Band  .2 u 

These  can  be  accomplished  by  using  a synthetic  sapphire  window  (AA^O.)  and  a 
filter  for  the  narrow  region.  ^ ^ 


The  following  table  gives  the  locations  of  the  Infrared  band  heads. 


H^O 


X 

quality 

6.2  u 

very  strong 

2.66  u 

very  strong 

1.14  u 

strong 

.94  u 

strong 

CO2 

X 

gual Itv 

15  u 

very 

strong 

4.26  u 

very 

strong 

2.77  u 

very 

strong 

2.69  u 

very 

strong 

The  broad  band  radiation  will  give  evidence  of  solid  particle  radiation 
as  well  as  any  chemi -luminescence  occurring  from  the  combustion  process  If  Is 
occurring. 
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ABBREVIATIONS.  ACRONYMS.  AND  SYMBOLS 


B 

C 

E 

E/V 

F 

I 


6 

Y 

6 

X 

P 

♦ 


magnetic  field 
speed  of  light 
electric  field 

specific  explosive  energy  density 

force 

Intensity 

debris  Impulse 

mach  number 

driver  pressure 

dynamic  pressure 

absolute  temperature 

velocity 

secondary  shock  velocity 

doppler  frequency 

radar  frequency 

Planck's  constant 

Boltzmann's  constant 

areal  density  of  debris  mass 

number  of  particles 

stagnation  pressure 

pitot  tube  pressure 

static  pressure  In  flow 

static  pressure  In  stagnated  region 

charge  on  particle 

DABS  flow  velocity 

voltage 

time  Intervals 

distance  between  gages 

ratio  of  Isobaric  to  Isovolumetric  specific  heats 
angle  between  secondary  shock  propagation  and  DABS  flow 
wavelength 
density 

angle  between  line  connecting  explosive  charge  and  gages  and 
flow  direction 
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